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HAUPTAUFSÄTZE 


Bemerkungen zur Theorie der freien Turbulenz. 
Von L. Prandtl in Göttingen. 


Zur Einführung für den nachfolgenden Aufsatz von H. Görtler wird ein 
Überblick über den derzeitigen Stand der Grundformeln der turbulenten 
Ausbreitung gegeben. Es wird dabei auch zu den von H. Reichardt 
vorgeschlagenen Ansätzen Stellung genommen. 


Die neue Formulierung, die Herr H. Reichardt im Zusammenhang mit seinen ver- 
dienstvollen Messungen über die turbulenten Ausbreitungsvorgänge') zur rechnerischen Ver- 
folgung dieser Vorgänge aufgezeigt hat?), lädt zu einer Betrachtung der hierdurch gegebenen 
Gesamtlage ein. 

Vor nicht ganz 20 Jahren standen die beiden älteren Ansätze für die scheinbare Schub- 
spannung in einer turbulenten Strömung, der aus strengen mechanischen Gesetzen hergeleitete 
von Reynolds 


und der rein phänomenologische von Bousinesq 

(w, v’ = Geschwindigkeitsschwankungen, @ = zeitlicher Mittelwert der Geschwindigkeit, e ein 
Maß der Turbulenzstärke, das die Dimension einer kinematischen Zähigkeit hat und deshalb 
„scheinbare kinematische Zähigkeit der turbulenten Strömung“ genannt werden mag?)) zu- 
nächst ohne Vermittlung nebeneinander. Der Verfasser hat mit seinem Mischungswegansatz 
erstmalig eine Brücke zwischen den beiden Formeln geschlagen‘). Die primitivste Formu- 
lierung ergab für die Bousinesqsche scheinbare kinematische Zähigkeit e die Beziehung 

ou 


l war eine neu eingeführte Länge, der „Mischungsweg“, dessen Größe in einem Einzelfall 
mit Hilfe von (3) aus Messungen der Geschwindigkeitsverteilung und Ermittlungen über die 
Schubspannung r gewonnen werden konnte. Wesentlich dabei war, daß in dem Falle, wo 
die turbulente Strömung einen mit der Geschwindigkeit quadratischen Widerstand erzeugt, 
! unabhängig von dem Betrag der Geschwindigkeit, also eine reine Ortsfunktion wurde. 


1) Gesetzmäßigkeiten der freien Turbulenz. VD1-Forschungsheft 414. Berlin 1942. 

?2) Ebenda, vgl. auch seine frühere Mitteilung in Z. augew. Math. Mech. Bd. 21 (1941), S. 257—264 „Uber eine 
neue Theorie der Turbulenz‘“. 

3) Das Produkt 0 € ist identisch mit dem von Wilh. Schmidt in die meteorologischen Betrachtungen eingeführten 
Begriff „Austausch“, vgl. sein Buch „Der Austausch in freier Luft und verwandte Erscheinungen“, Hamburg 1925. 

4) Untersuchungen zur ausgebildeten Turbulenz. Z. angew. Math. Mech. Bd. 5 (1925), S. 136. 
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Formel (3) hat jedoch, abgesehen davon, daß diese Ortsfunktion vielfach keinen genügen«l 


einfachen Verlauf zeigt, noch das Unbefriedigende an sich, daß dort, wo = Null wird, also 


am Ort eines Maximums oder Minimums von “, bei endlichem ? die scheinbare Zähigkeit 
Null wird. Dies kommt auch in den an diesen Stellen unnatürlich zugespitzten Geschwin- 
digkeitsverteilungen, die bei regulärem Verlauf von r(y) und ! mit dieser Formel erhalten 
werden, zum Ausdruck. Der Verfasser hat deshalb schon in der erwähnten ersten Ver- 
öffentlichung einen Ausweg vorgeschlagen, der darauf hinauskommt, daß man die Formel (3) 
nicht auf den zeitlichen Mittelwert der Geschwindigkeit, sondern auf die augenblickliche Ge- 
schwindigkeit anwendet, deren Maximum wegen der Schwankungen nach Zufall seitlich ver- 
lagert wird, und daß man von diesen Augenblickswerten von e nun an jeder Stelle den qua- 
dratischen Mittelwert bildet. Führt man eine weitere Länge !’ ein, die positiv oder negativ 
sein kann und zeitlich nach Zufallsgesetz variiert, so kann man in erster Näherung 


9u_ du, 
dy 
setzen. Der erwähnte Mittelwert ist dann 


(4) 


/ 


Man kann sich leicht davon überzeugen, daß die Formel, die aus Formel (3) entsteht, wenn 
ou 
suchstatsachen weitgehend gerecht wird. Man kann allerdings dagegen auch einwenden, daß 
jetzt zwei empirische Größen zur Verfügung stehen, mit denen man sich natürlich einem 
gegebenen Versuchsmaterial besser anpassen kann als mit einer einzigen. Die verbesserte 
Formel ist freilich für die wirkliche Durchführung analytischer Rechnungen sehr unbequem, 
so daß der Verfasser in der Zwischenzeit öfters Versuche gemacht hat, für praktische Zwecke 
eine bequemere Formulierung zu finden. Das neue Versuchsmaterial, das inzwischen Herr 
H. Reichardt erarbeitet hat, ließ diese Gedankenreihe wieder aufleben und ermunterte zu 
einem neuen Anlauf in dieser Richtung. Als Ergebnis dieses Suchens ergab sich ein Ansatz, 
der, wie gleich betont werden mag, ziemlich roh ist, dessen Wert aber in seiner Ein- 
fachheit und Handlichkeit liegt, durch die es gelingt, auf einfacherem Wege als es mit dem 
Ansatz (3) möglich ist’), die wichtigsten formelmäßigen Beziehungen zur Theorie der freien 
Turbulenz zu gewinnen. Die Einfachheit hat natürlich einen Verzieht an Genauigkeit in den 
Einzelheiten im Gefolge. 

Bei dem in Rede stehenden Ansatz ist davon ausgegangen, daß die einzelnen Flüssig- 
keitsteile, die bei der turbulenten Vermischung quer zur Strömungsrichtung wandern, durch- 
aus nicht, wie es der Ansatz (3) annalım, sozusagen differentieller Natur sind, sondern daß 
sie vielmehr in ihrer Größenordnung an die Querabmessung der Vermischungszone heran- 
reichen. Es wird deshalb eine scheinbare kinematische Zähigkeit in der Weise gebildet, daß 
der maximale Geschwindigkeitsunterschied der gemittelten Strömung mit einer passenden 
aus der Breite der Vermischungszone entnommenen Länge multipliziert wird, wozu noch ein 
Zahlenfaktor kommen muß, dessen Wert aus den Versuchen bestimmt wird. Es wird also gesetzt 

Aus diesem Ansatz ergibt sich unmittelbar, daß e über den einzelnen Querschnitt konstant 
gesetzt ist, während es nach Formel (3) über den Querschnitt veränderlich wird, wenn, wie 
gewöhnlich, ! konstant angenommen wird. Besonders große Abweichungen gegenüber den 
älteren Rechnungen werden am Rande der Vermischungszone zu erwarten sein, da nach 
Formel (3) wie auch nach Formel (5) e dort gegen Nuil geht. Die Versuchsergebnisse weisen 
in die gleiche Richtung). 

Die Durchführung einiger Aufgaben mit Hilfe des Ansatzes (6) hat auf Wunsch des 
Verfassers Herr H. Görtler übernommen’). Seine Rechnungen, die auch einiges mathe- 
matisches Interesse bieten, und ein Vergleich ihrer Ergebnisse mit den Reichardtschen 
Messungen finden sich in seinem hier anschließenden Aufsatz. 


an Stelle von. der Ausdruck auf der rechten Seite von (5) eingeführt wird, den Ver- 


5) Vgl. etwa die zum Teil ziemlich verwickelten Rechnungen von Herrn W. Tollmien, Z. angew. Math. Mech. 
Bd. 6 (1926), S. 468478. 

6) Vgl. auch die Reichardtschen Ausführungen über den Verlauf der Austauschfunktion der freien Turbulenz 
in dem demnächst in dieser Zeitschrift erscheinenden Aufsatz „Uber Impuls- und Wärmeaustausch in freier Turbulenz“. 

?, Hier ist anzumerken, daß Herr Reichardt, ohne daß es dem Verfasser bekannt war, den Ansatz (6) zur 
gleichen Zeit auch in Erwägung gezogen hat und zwei einfache Beispiele dafür gerechnet hat, ihn dann aber zugunsten 
seines veröffentlichten Ansatzes wieder verworfen hat, da dieser die Versuchsergebnisse mit höherer Genauigkeit 
wiedergibt. (Die erwähnten Beispiele betrafen die ebene Strahlausbreitung und den ebenen Nachlauf, vgl. das zweite 
und dritte Beispiel von Herrn Görtler.) 
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Zu den Reichardtschen Betrachtungen sei zur Orientierung des Lesers das Folgende 
ausgeführt: 

In seinen „Gesetzmäßigkeiten der freien Turbulenz“ werden die bei freier Turbulenz 
bisher beobachteten Geschwindigkeitsprofile kritisch gesichtet und besonders sorgfältige 
_ eigene Messungen hinzugefügt. Andererseits wird unter Ausnützung aller theoretisch-hydro- 
dynamischen Beziehungen die Aufgabe gelöst, die praktisch nicht meßbaren Größen, wie die 
Quergeschwindigkeit und die in der scheinbaren Schubspannung nach (1) auftretende Größe 
wv’, aus den allein gemessenen Gesamtdruckprofilen #? (,) rechnerisch zu ermitteln, wobei 
einerseits „ausgebildeter Strömungszustand“, d. h. Affinität der Geschwindigkeitsprofile in 
allen Abständen & des Bereiches, vorausgesetzt und in erster Näherung Druckkonstanz an- 
genommen wird. Aus den Versuchen ergibt sich, daß die Gesamtdruckprofile mit recht guter 
Annäherung durch Fehlerfunktionen (bzw. Fehlerintegralfunktionen) dargestellt werden können. 
Damit können die in dem erwähnten Formelsystem auftretenden Differentiationen und In- 
tegrationen bequem ausgeführt und so die nicht meßbaren Größen ermittelt werden. 

In einem weiteren Schritt hat nun Herr Reichardt ein einfaches System von Diffe- 
rentialformeln aufgestellt, durch das alle bisher untersuchten Fälle von freier Turbulenz zu- 
sammengefaßt werden konnten. Die erste dieser Gleichungen ist bekannt. Sie wird aus der 
Eulerschen Gleichung für eine volumbeständige Flüssigkeit durch Bildung der zeitlichen 
Mittelwerte erhalten und bildete schon die Grundlage der alten Reynoldsschen Be- 
trachtungen. Sie lautet für ebene Strömungen 


/p = 7 


Bei der freien Turbulenz kann in (7), wie auch sonst üblich, in eıster Näherung Po ge- 


setzt werden, wodurch die Gleichung zu einer „Kontinuitätsgleichung des «-Impulses“ wird. 
Hierzu fügt Reichardt eine zweite Gleichung hinzu, die phänomenologischen Charakter 
trägt und die lautet: 


- 


Hierin spielt offenbar die „Übertragungsgröße“ A (eine Länge, die aber nicht mit dem 
Mischungsweg identisch ist) eine ähnliche Rolle wie die scheinbare kinematische Zähigkeit & 
in Gl. (2). Aus Gl. (7) und (8) läßt sich «® fortschaffen und es bleibt dann 


was bei A=const. mit der partiellen Differentialgleichung der eindimensionalen nicht 
stationären Wärmeleitung übereinstimmt, wenn an Stelle von » die Zeit # tritt und y die 
räumliche Entfernung bedeutet. Die Fehlerfunktion und das Fehlerintegral gehören be- 
kanntermaßen zu den Lösungsfunktionen dieser Gleichung. Bei der Turbulenzaufgabe ist A 
im allgemeinen nicht konstant, sondern eine Funktion von x, deren Bau sich durch einfache 
Impulsbetrachtungen od. dgl. aus der Forderung der Affinität der Geschwindigkeitsprofile 
herleiten läßt. 

Der Verfasser fand es seinerzeit bedenklich, daß in Gl. (9) nicht die Beschleunigung 

2 

en, sondern nur a stand, was die «-Achse unerwünschterweise auszeichnete. 
Er empfand dies als eine Verletzung der Newtonschen Relativität, nach der die Kräfte in 
einem mechanischen System von der Zufügung einer konstanten Geschwindigkeit (hier einer 
solchen in der y-Richtung) unabhängig sein müssen. In (8) hätte man deshalb eher den von 
einer konstanten Zusatzgeschwindigkeit unabhängigen Ausdruck «v’ erwarten sollen als 
ue=uv-+wv. Auch Herr Reichardtgibt an, daß er das Hinzutreten von « v als merkwürdig 
empfand. Der bei Reichardt angegebene Funktionsverlauf von @v und wv’ zeigt Jedoch, 
daß man das Glied «v nicht entbehren kann, wenn man die Fehlerfunktion für die Gesamt- 
druckverteilung erhalten will. In den Reichardtschen Beispielen ist allerdings jedesmal 
eine bestimmte Strömungsrichtung durch die Problemstellung ausgezeichnet, die dann folge- 
richtig als «-Richtung definiert wird, so daß die nachträgliche Zufügung einer konstanten 
v-Geschwindigkeit nicht mehr in Betracht gezogen werden darf. Herr Reichardt beruft 
sich auch im Hinblick auf den Erfolg seiner Rechnungen mit einem gewissen Recht hierauf; 
wenn man aber an allgemeinere Beispiele denkt, z. B. an Nachlauf oder auch Strahl- 
ausbreitung in einer gekrümmten Strömung, dann wird man andere Wege gehen müssen. 
Die Ansätze, die von der scheinbaren Schubspannung ausgelren, sind von diesem Einwand frei 
(erfordern aber freilich eine Integration der hydrodynamischen Differentialgleichungen!). 419 
. 16* 
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Berechnung von Aufgaben der freien Turbulenz auf Grund 
eines neuen Näherungsansatzes. 


Von H. Görtler in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Es werden einige Probleme der freien Turbulenz auf Grund des in den 
vorangehenden Darlegungen ron L. Prandtl erörterten neuen Näherungs- 
ansatzes: „Austausch über den einzelnen Querschnitt konstant und propor- 
tional dem Produkt aus (reschwindigkeitsunterschied und Vermischungs- 
breite“ durchgerechnet. In allen behandelten Fällen ergibt sich gute 
Übereinstimmung der berechneten (reschwindigkeitsverteilungen mit den 
Beobachtungen, abgesehen natürlich von den unmittelbaren Randgebieten 
der Vermischungszonen, wo ja der Austausch in Wirklichkeit erlischt. 


Auf Anregung von Professor Prandtl werden im folgenden einige Probleme der freien 
Turbulenz auf Grund eines neuen Näherungsansatzes für die scheinbare Schubspannung in 
einer turbulenten Strömung: 

y 
mit 
durchgerechnet, (@ zeitlicher Mittelwert der Geschwindigkeit, e scheinbare kinematische Zähig- 
keit, b Maß für die Breite der Vermischungszone, # dimensionsloser Proportionalitätsfaktor), 
d. h. die kinematische Zähigkeit wird über den einzelnen Querschnitt konstant angenommen. 
Wie Professor Prandtl in seinen diesem Aufsatz vorangehenden Darlegungen ausgeführt 
hat, ist zu erwarten, daß gegenüber dem alten Ansatz 


der in den Extremstellen von @ zu unnatürlich zugespitzten Profilen geführt hat, weil dort 
nach (3) & verschwindet, der Ansatz (2) im Inneren der Vermischungszone zu guter Überein- 
stimmung mit dem wirklichen Verlauf führen wird, daß aber am Rande der Vermischungs- 
zone Abweichungen zu erwarten sind, weil dort der Austausch in Wirklichkeit erlischt. Dies 
wird durch die folgenden Rechnungen bestätigt. Um der Rechnung Grenzen zu setzen, habe 
ich mich auf die interessierenden ebenen Probleme der freien Turbulenz beschränkt. 


) 
1. Die ebene Vermischungszone zweier Strahlen. 


Zwei Parallelstrahlen mit konstanter mittlerer M 
Geschwindigkeit U, bzw. U, mögen durch eine 2 
Wand (negative «-Achse, siehe Bild 1) getrennt 
sein, die in x=0 abbricht, so daß sich von hier 
an eine turbulente Vermischungszone der beiden 
Strahlen ausbildet. Gefragt ist nach der Verteilung 
der mittleren Geschwindigkeit für x >0. 


Wie in allen folgenden Beispielen soll der 
Druck in erster Näherung als konstant angesehen 
werden. Die Pewegungsgleichung und die Konti- [7A 
nuitätsbedingung lauten dann allgemein: Bild 1. Ebene Vermischungszone zweier Strahlen. 

du du 


ip’ 


da), 


wobei wir der Einfachheit halber den Querstrich bei @, ® hier und im folgenden fortlassen. 
u bedeutet also den zeitlichen Mittelwert an der betreffenden Stelle. 
Im vorliegenden Fall ist aus dimensionsanalytischen Gründen b=c.r, wo ce eine Kon- 
stante bezeichnet, und somit 


Die Geschwindigkeitskomponenten « und v sind als reine Funktionen von ) = yj.x anzusetzen. 


| Pa 
4 
dar y 
XUM 
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Führen wir zur Befriedigung der Kontinuitätsgleichung die Stromfunktion y vermöge u — 2 


tr ein und setzen wir N 
so hat y die Gestalt 
Es ist dann 
und durch Einsetzen in die nun allein noch zu befriedigende erste Gl. (1.1) ergibt sich als 
Differentialgleichung') für F: 
+20? F F’—=0 
wo 


1 
mit 
An Randbedingungen ist zu fordern: 
1+4 für y=o, 


F'n)= 
E; für „= 


damit « die vorgegebenen Werte außerhalb der Vermischungszone annimmt. Gleichzeitig wird 
v—=const für „= +». 
Für unsere Ausführungen erweist es sich als zweckmäßig, an Stelle von 7 als unab- 


hängige Veränderliche endgültig 


einzuführen. Wir machen dann für F den Reihenansatz’) 


Setzt man diese Reihe in (1.6) ein, so ergibt der Koeffizientenvergleich gleicher 4-Potenzen 
das folgende rekursiv zu lösende System von Differentialgleichungen für die Beiwertfunktionen 
F,(&) e#=1,2,3,...: F,=2£), in welchem Striche Ableitungen nach £ anzeigen: 
Fi 
FÜ +2EFY= 
Fir = 


mit den Randbedingungen 


F,(&)=0 für 
Es bleibt dann in den Lösungen F,, F,,.... eine additive Konstante willkürlich. Die Tat- 
sache, daß hier eine Unbestimmtheit in der Lösung des Problems bestehen bleibt, ist nur 
Ausdruck einer grundsätzlichen Unbestimmbarkeit der Lösungen der Grenzschichtdifferential- 
gleichung. Wir kommen später darauf zurück. 
Zunächst berechnet sich F, unter Erfüllung der Randbedingungen (1.10a) zu 


. (1.10a). 


0 


@ 


P(a)= \ 


0 
das Fehlerintegral und d, eine Integrationskonstante bezeichnen. 


1) Diese Differentialgleichung liegt auch der Berechnung der laminaren Grenzschicht an der Platte zugrunde, 
natürlich bei anderen Randbedingungen, siehe H. Blasius: Grenzsehichten in Flüssigkeiten mit kleiner Reibung. 
Z. Math. u. Phys. Bd. 56 (1908), S. 1 bis 37. Ich schlage zur Lösung jedoch einen anderen Weg ein, der hier wesentlich 
einfacher zum Ziele führt. 

2) Dieser Reihenansatz entspricht folgender Iteration: Man gehe aus von der mittleren Strömungsgeschwindig- 
keit u=(U, +U3) und also von F=F®=n(=o-1F,) als nullter Näherung. Für eine erste Verbesserung 
F= F«0 + Fcı> bestimmt sich die Korrektur F‘'’ unter Vernachlässigung quadratischer Glieder in F<!> und in dessen 
Ableitungen aus (rn). Dabei muß das Verhalten von u für n=+» in Ordnung bringen 
(u=U, bzw. =U3). Durch Fortführung dieses Iterationsprozesses gelangt man zur Reihe (1.9). 


ch. 245 
y 
d 
en 
[S- 
(te 
en 
en 
en 
in 
1) 
>) 
g- 
r), 
rt 
3), 
rt 
n- 
;S- 
es 
be 
für 
—1 für E=—o; 
1), 
n. 
n- 
n. 
XUM 


2. . Math. Mech. 
246 Görtler, Berechnung von Aufgaben der freien Turbulenz Be Nr.5 Okt. 108 


Für. die weiteren Gleichungen des Systems (1.10) wählen wir zweckmäßig die Schreib- 


weise 
indem wir den Faktor 
(Fehlerfunktion) rechts herausziehen. Dann ist das Integral Fy’ von (1.12) formal 


und somit durch Quadraturen hieraus 


3 
F,=\®,(ß) \ 
0 0 


=| | ®,(P)\ 
0 0 


C,, D, und d,(=2) sind Integrationskonstanten, die mit dem ebenfalls noch unbestimmten 
d, festgelegt werden müssen. Dies und die rekursive Ausrechnung der Hilfsfunktionen R, (£) 
für den Formelapparat (1.13) ist die allein noch zu bewältigende Aufgabe. Es sind ® (£) und 
D,(£) bekannte tabulierte Funktionen’). Zur Berechnung der R, hat man dann, wie sich 
zeigen wird, nur noch Quadraturen auszuführen, und dasselbe gilt nach (1.13) dann auch 
für die interessierenden Größen des Strömungsfeldes. Dies bedeutet für das numerische 
Rechnen einen erheblichen Vorteil. 

Die Integrationskonstanten C,, D, und d, sind vollkommen festgelegt, wenn wir neben 
den Randbedingungen (1.10a) noch etwa fordern, daß «(0)= U, d. h. F’(0)=0, also 


Daß diese Forderung keine unerwünschte Einschränkung unserer Betrachtungen bedeutet, 
wird noch ersichtlich werden. Für r=1 ist bereits von selbst F/ (0)=0 erfüllt. 
Es wird nun zunächst 


wobei F\, nach (1.11) bis auf eine additive Konstante d, bekannt ist. Setzen wir diesen Aus- 
druck in (1.13) ein, so erhalten wir 


00 


( wobei die Beziehung \®. ()ada=—— > Io, (&) — NZ \ benutzt wurde), Da die beiden 
ersten Ausdrücke gerade Funktionen in £ darstellen, folgt aus F}(+ ®)=0, daß 
? 


Wegen F;(0)=0 muß ferner 

sein, also 


00 
(Es gilt allgemein, daß d, erst mit C,,,, D,;+, durch die Randbedingungen für F}; ,(£) be 
stimmt wird.) 
Für v=3 ergibt sich zunächst wegen 
{F,F} + F 


0 


3) E. Jahnke und F. Emde: Funktionentafeln, 2. Aufl. Leipzig 193, S. 97 ff. 
4) Hier und im folgenden wurden die Integrationen graphisch ausgeführt. 


. (1.13a) 
_ 
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Wegen 


D,(P) \ R,(a) da 2a, | . . . (1.20) 
0 


mit R,—=R, 2d,®,(£), ergibt die Forderung F; (0)=0 zunächst 


und F}/(+x)=0 kann nur befriedigt werden, wenn der 2. Term als gerade Funktion ver- 
schwindet, also 

d,—V x 
da ja die übrigen Terme ungerade Funktionen in & sind, und wenn 


3 


. 2.2... (LM. 
0 0 


In dieser Weise kann man bequem fortfahren. Es möge noch das Glied mit v»—=4 be- 
rechnet werden. Es wird sich zeigen, daß die Konvergenz der Entwicklung auch im un- 
günstigsten Falle so gut ist, daß dieses Glied und im allgemeinen bereits F, praktisch keine 
Beiträge mehr liefern. 

Da ist, wird 


0 


wobei F, die noch unbekannte Konstante d, enthält. Wegen 


| R(a)dadß— 24, | . . . (1.3) 
0 0 v0 


mit R, —2d,®,(£), schließt man zunächst aus (0)=0, daß 

sein muß, und aus F}(+x)=0, daß die ungerade Funktion ©, P(£) verschwinden, also 


sein muß, da die übrigen Glieder gerade Funktionen in & sind, und daß 
3 
0 0 


Die noch unbestimmte Integrationskonstante d, in F, würde sich erst bei Ermittlung 
F, ergeben. Man überlegt sich jedoch leicht, daß d,=0 ist, 


Zahlentafel 1. 


Fy F, F, 


0,007 0.007 — 0.001 0,008 0,000 0,003 
0,025 0010 0,001 0,013 0,001 0,004 
0.044 0,009 0,004 0,013 0.002 0,002 
0.058 0005 0,006 0,008 0.002 0,001 
0,062 — 0,001 0,007 — 0,001 0,002 0,008 
0,056 — 0,005 0,006 — 0,009 0.008 0,007 
0,044 — 0,007 0,004 — 0,015 0.005 0,010 
0.030 — 0,006 0,000 — 0,017 0.007 0.008 
0,019 — 0,005 — 0,002 — 0,015 0.008 0.002 
0.010 — 0,008 — 0.005 — 0,011 0.008 — 0,005 
0.005 — 0,002 — 0,007 — 0,007 0,006 — 0,009 
0,002 — 0,001 — 0,008 — 0,004 0,005 — 0,009 
0,001 0,000 — 0,009 — 0,002 “0,008 — 0,008 
0,000 0,000 — 0,009 — 0,001 0.002 — 0,005 
0.000 0,000 — 0,009 0,000 0,001 — 0,003 


| 
247 
| 
| | 
| 
n 
d 
h 
e 
). 
), 
1st. 
von 
& F, Fi F, 
| -0,161 0 1.128 0 
0,2 | — 0,139 0,223 1,084 0,000 
0.4 | — 0,073 0,428 0,962 0,003 
0,6 0,031 0,604 0,787 0,010 
| 0,8 0,166 0,742 0,595 0,021 
1.0 0,328 0,843 0,415 0,033 
1.2 0.503 0,910 0,267 0,045 
14 0,690 0,952 0,159 0,055 
1,6 0,883 0,976 0,087 0,062 
1.8 1.080 0,989 0,044 0,067 
2,0 1,278 0,995 0,021 0,070 
2,2 1,477 0,998 0,009 0,071 
2,4 1,677 0,999 0,004 0,072 
2,6 1,877 1,000 0,001 0,072 
2,8 2,077 1,000 0,000 0,072 
3,0 2,277 1.000 0,000 0,072 
XUM 
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Die numerische Auswertung unserer 
Ergebnisse liefert die Zahlentafel 1 für 
F,, F, und P=1,2,3,4; für ist 
F,=&). Es ist dabei zu beachten, daß F,, 

+1, Fy gerade oder ungerade Funktionen 

von & sind, je nachdem » ungerade oder ge- 

rade ist. Auf Grund dieser Tafel berechnet 
| | man «(£) vermöge 


05 


u 
Bild 2. Ebene Vermischungszone zweier Strahlen: Verlauf der er => ”F,(). . (1.29), 
Funktionen für v»=1,2,3,4. Es ist ferner 1. 

wobei O</=1. Im Falle U,—=0 
(Vermischungszone zwischen Strahl 
und ursprünglich ruhender Luft) 
ist 4—=1. Selbst in diesem für die 
Konvergenzgüte ungünstigsten Fall 
liefern die Glieder F} und F} keine 
wesentlichen Beiträge mehr (siehe 
Bild 2). Um so mehr trifft dies zu, 
wenn /<1 ist, da die Potenzen 
4” die Güte der Konvergenz ver- 
bessern gegenüber dem Fall /=1. 
Entsprechend ist 


7,0), 


y=alo'Zu—v). 


| | 
-91 | 17777] Wirschreiben 77 2 ’v,(d), 


Bild 3. Ebene Vermischungszone zweier Strahlen: Verlauf der Funktionen und zn ve Bild 3 zuer Funk- 
für v—1,2,3,4. tionen v,(&) für v—=1,2,3, 4. 


Nun noch ein Wort zur Unbestimmtheit der Lösung, welche oben durch die zusätzliche 
Forderung u (0) = U willkürlich behoben wurde. Es ist eine wohlbekannte Eigenschaft der 

“ar? 
du dv 


daß mit einer Lösung « (x, y), » (x, y) auch 
und das zugehörige »* (,y) auf Grund der Kontinuitätsgleichung: 
y + Flo) y + 


eine Lösung darstellen. Das gilt nicht nur im laminaren Fall, wo e=»=const ist, sondern 
auch im vorliegenden Falle des turbulenten Austausches, wo e=e(x) im allgemeinen va- 
riabel ist. 

In Fällen wie dem vorliegenden, in welchen eine bestimmte Stromlinie nicht vorgegeben 
ist, bleibt mithin die Lösung bis zu einem gewissen Grade unbestimmt. Die einleitend ge- 
stellte Forderung, daß v/u = const für y= +» sein soll, ergibt f(x) =ax (a eine Konstante). 
Wegen =oy].x ist somit neben der von uns hergestellten Lösung «(£), »(&) auch 


eine Lösung. Man erkennt leicht, daß die in unserer obigen Rechnung zutage getretene Un- 
bestimmtheit gerade diese Willkür zum Ausdruck bringt. Mit der zusätzlichen Forderung: 
u = U haben wir somit lediglich das Geschwindigkeitsprofil «(£) gegenüber noch freien 


5) Siehe L. Prandtl: Zur Berechnung der Grenzsehichten. Z. angew. Math. Mech. Bd. 18 (1038), S. 77 bis 9. 
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kleinen Verschiebungen in £, also gegenüber kleinen Drehungen der Vermischungszone um 
den Ursprung <= y=V0, festgelegt®). Auf die Gestalt des Profiles «(£) hat dies natürlich 
keinen Einfluß. 

Die Forderung « (0) = !/ ist lediglich gewählt worden, weil sie die Konvergenz unserer 
Entwicklungen besonders begünstigt. Denn es ist a v= und 
da von selbst F/(0)=0 ausfällt, müßten die höheren Näherungen F}, F}, ... die Korrektur 
für jede Verschiebung des Profils «(£) aus der durch «(0) = U gekennzeichneten Lage zu: 
sätzlich aufbringen. So aber, indem wir u (0) = fordern, bringen die höheren Näherungen 
nur noch kleine Korrekturen am Rande der Vermischungszone (siehe Bild 2), während der 
mittlere Teil schon durch die ersten Glieder der Reihe gut wiedergegeben wird. 


Man hätte unter Verzicht auf die gute Konvergenz statt «(0)=TU auch z. B. » (0) =0 
fordern können, d.h. die Gerade y=0 zur Stromlinie machen können. Hieraus wird ersicht- 
lich, daß in jeder Vermischungszone zweier paralleler Strahlen die von 2=V, y=0 aus- 
gebende Stromlinie eine Gerade ist, und zwar in der obigen Schreibweise die Gerade y=— ax '). 


Zahlentafel 
erste Nüherung 


vlU, uU, vlU, 


0,081 
0,093 0,388 0,092 
Bild 4. Ebener Strahlrand (U;=0): Geschwindig- 0,128 0,292 0,121 
keitsverteilung. ). 0.176 0.214 0.159 

0,227 0,155 0.201 
— 0,269 0.111 0,239 
o Messungen von 
Meichardt 0,300 2 | 0,079 0,274 
0,314 0.056 0,304 
0,323 0,039 0.329 
0,324 0,026 0,350 
0,017 0,368 
0,010 0,382 

Bild 5. Ebener Strahlrand: Vergleich mit Druckmessun- _ 0,003 0,399 
gen von H. Reiehardt. Auch die Lage des Profil- 


wendepunkts stimmt mit dem Ergebnis Reiehardts . 0 0,408 
(nach besonderen Druckdifferenzenmessungen) überein. 


T 


Im 


Bild 4 und Bild 5 beziehen sich auf den Fall =1 (d.h. U,=0). In Bild 4 ist das 
Geschwindigkeitsprofil « (£) für A=1 gezeichnet (siehe auch Zahlentafel 2). Der Wendepunkt 
des Profils liegt nicht in &=0 (wo u—=T/,/2), sondern schwach nach der Seite höherer Ge- 
schwindigkeiten verlagert bei £=0,15 (als Nullstelle von F(£) ermittelt). Zum Vergleich ist 


gestrichelt die erste Näherung u/TU, -4.(d+ F}) eingezeichnet. Sie läßt einmal die Klein- 


heit der höheren Glieder der Entwicklung in Erscheinung treten, zum anderen zeigt sie das 
gesetzmäßige Abweichen der Funktion «/U,—"/s von einer in & ungeraden Funktion, indem 
die ausgezogene Kurve links und rechts von £=0 oberhalb der gestrichelten bleibt. In Bild 4 
ist auch der Verlauf von v (£) eingetragen so, wie er sich für unsere spezielle Lösung ergeben 
hat. In Bild 5 ist w’(£)/U} aufgetragen. Auch hier liegt der Wendepunkt nicht dort, wo 
='l2(£=0,35), sondern etwa bei £=0,4 (als Nullstelle von ermittelt. Dies 
steht nicht nur qualitativ sondern auch quantitativ in völliger Übereinstimmung mit dem auf 
Grund von Druckdifferenzenmessungen von H. Reichardt*) gewonnenen Ergebnis. Aber 


6) W. Tollmien fordert bei seinen entsprechenden Reehnungen nach dem alten Ansatz (3) für den Falli=]1 
statt dessen: v—0 für y„>+»; siehe W. Tollmien: Berechnung turbulenter Ausbreitungsvorgänge. Z. angew. Math. 
Mech. Bd. 6 (1926), S. 468 bis 478. 


x) v(x,0) ve(X,—ax) _ 
ur(X,y) ulX,y+ax) 


—a, also mit ——=0 auch =—a. 
1X, 0) ur (x,— ax) 


“ s) H. Reichardt: Gesetzmäßigkeiten der freien Turbulenz. VDI-Forsch.-Heft 414 (1942). Vgl. dort Bild 2 auf 
ite 9. 


h. 
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st 
n 
st 
v0 
ıl 
t) 
e 
1 
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Ur 
Wendepunkr. 
1, | | | 
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n 
P- 
). 
1° 
n 
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auch der Gesamtverlauf der in Bild 5 wiedergegebenen Kurve stimmt gut überein mit den ein- 
getragenen Meßpunkten nach Beobachtungsmaterial von H. Reichardt (abgesehen natürlich 
von den Rändern der Vermischungszone). Nach Reichardt a.a.O. ist die Breite der Ver- 
mischungszone —- gemessen als Entfernung jener Stellen &, wo (u/U,) gleich 0,1 (£= — 0,345) 
und 0,9 (£=0,975) ist — durch 0,098 x gegeben. Daraus ergibt sich (für diesen Fall U, = 0) 
der in unserer Theorie freibleibende Streckungsfaktor zu o=13,5. 


2. Das ebene Problem der Strahlausbreitung. 


Eine Flüssigkeit ströme aus einem engen Spalt (für die Rechnung eine singuläre Linie) 
und vermische sich mit der umgebenden ruhenden Flüssigkeit. Sieht man den Druck wieder 
als konstant an, so ist der Impuls in «-Richtung (Hauptrichtung des Strahles) 


| u?dy= const. 


Macht man für « den Ansatz «= so ist daher — alse mit „= f(y)=F’(n): 


y=yzxF(n), 

Ist b, die gewählte charakteristische Breite des Strahls in der Entferung «=s vom 

Austrittsspalt, so ist ferner ax — = und also mit 


Zu erfüllen ist demnach mit dem Ansatz (2.1) die Differentialgleichung 
V 


die auf 


führt. Zweimalige Integration führt auf 
(2.4), 


wenn man die Zählung der Stromlinien y = const durch die Forderung F(0) = festlegt und 
ferner aus Symmetriegründen verlangt, daß F’ (n) eine gerade Funktion in n sein soll. U, be- 
zeichnet den Wert von win “=s, y=0. Die Differentialgleichung (2.4) besitzt dann das 
durch F(0) = festgelegte Integral 


yet, 
o= . . . . . . . . . . . . (2.6), 


welches auch den Forderungen des Problems: F’(+x)=0 von selbst gerecht wird. Dem- 
nach ist?) 


mit 


zangon, 


> 6of-: 
Es 


v s 1 
T, = V: Coj?on— 35 Zangen) 
o ist den Experimenten zu entnehmen. Nach den Messungen von H. Reichardt (a.a. 0. 


S. 13) ist der Abstand jener beiden Stellen des Profils v: 7 ° mit dem Wert 0,5 durch 
s 


0,155 gegeben. Da nach (2.7) diese Stellen durch 07 = + 0,606 gegeben sind, ermittelt man 
hieraus: o = 7,67. 


9) Diesele Lösung hat, unabhängig von diesen Rechnungen, H. Reichardt angegeben (a.a.0., S.20, Fußn. 2). 


| 
XUM 
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Es sei bemerkt, daß das Problem der ebenen laminaren Strahlausbreitung auf denselben 
Differentialgleichungstypus führt'°), und daß demnach das nach der vorliegenden Theorie 
resultierende Geschwindigkeitsprofil im turbulenten Strahl zu jenem des laminaren Strahls 
affın ist. (Dort wächst die Strahlbreite allerdings nicht proportional zu x sondern zu x "z!) 

In Bild 6 bzw. 7 haben wir das Geschwindigkeitsprofil « (x, y) (siehe auch Zahlentafel 3) 
und das Stromlinienfeld in dimensionsloser Auftragung dargestellt. Zum Vergleich ist im 
Bild 6 auch das von Tollmien'') nach dem Ansatz (3) berechnete Profil gestrichelt ein- 


| 


Bild 6. Ebener Strahl: Geschwindigkeitsverteilung. 
Gestaltlicher Vergleich mit dem u-Verlauf nach 
Tollmien. 


0 


Bild 7 (rechts). Ebener Strahl: Stromlinienverlauf 
(affin verzerrtes Bild!). Zahlentafel 3. 


1 
0,786 
0,420 
0,181 
0,071 
0,027 
0.010 
0,004 


0,001 
Bild 8. Ebener Strahl: Vergleich mit Druckmessungen ; 
von H. Reichardt. 0 


0 
2 03 2 


gezeichnet. Es muß dabei beachtet werden, daß die Streckungszahlen o und (2c?)”"Is (es ist c 
die im Tollmienschen Ergebnis freibleibende Integrationskonstante) in diesem Bilde völlig 
willkürlich einander gleichgesetzt wurden. Es kommt uns hier nur auf einen allgemeinen 
gestaltlichen Vergleich an, wobei die in der Einleitung genannten Merkmale der Ergebnisse 
beider Theorien deutlich in Erscheinung treten. 


2 
In Bild 8 geben wir schließlich einen Vergleich unserer Kurve für y: Br mit Be- 


obachtungsmaterial von H. Reichardt. Die Übereinstimmung ist wiederum abgesehen 
von den Randgebieten gut'?). 


3. Nachträgliche Berechnung der Druckdifferenzen in erster Näherung. 


War bisher in den behandelten Beispielen zur Vereinfachung Konstanz des Druckes 
angenommen, so wollen wir nun nach dem Vorgang von Tollmien'*) nachträglich die Druck- 
unterschiede in erster grober Näherung ermitteln. Die Kleinheit der sich ergebenden Druck- 
differenzen wird die Annahme konstanten Druckes in unseren bisherigen Entwicklungen als 
gute Näherung bestätigen. 

Der von Tollmien eingeschlagene Weg besteht darin, die Druckdifferenzen, welche 
durch Impulsänderung auf Grund der oben berechneten Geschwindigkeiten resultieren, zu er- 
mitteln, wozu die bisher gestrichene zweite Bewegungsgleichung 


” 10) Siehe S. Goldstein: Modern developments in fluid dynamies. Bd. ], S. 145 ff. bzw. die dort zitierte Original- 
iteratur. 

11) W. Tollmien:a.a.0., S. arıfl. 

12) Im Gegensatz zu einer Bemerkung von H. Reichardt a.a.O., S. 20, Fußnote 25. 

13), W. Tollmien: a.a.O., S. 476 ft. 


t. 1942 
ein- 
rlich 
Ver- 
345) 
()) 
4 
inie) | \als | | 
| | N | 
| | af 2 1 
| | 
| | 
vom 
13 / 
| Theorie | | 
0 0 
2.2), | | | | | 1,0 0,039 
2.3) | | | | | | | 2,0 — 0,341 
25 -— 0,427 
Ber | 3,0 — 0,468 
2.4), | 3,5 — 0,486 
— 0,494 
und —05 
be- 
das 
2,5) 
2.6), 
)em- 
27). 
. 
1. 25). 
XUM 
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dient. (Diese Rechnung ist als erster Schritt eines Iterationsverfahrens aufzufassen, wobei 
nunmehr unter Berücksichtigung des so gewonnenen Druckverlaufs die Geschwindigkeiten 
aus der ersten Bewegungsgleichung erneut — in zweiter Näherung also — berechnet werden, 
hieraus dann mittels (3.1) den Druckverlauf in zweiter Näherung usw.) 

1 or 


Im Falle des ebenen Strahlrandes erhält man er PLFF" und somit 


"2 


was wegen 20° FF" = — F’" (siehe (1.6)) auf 


führt. Striche bedeuten dabei überall Ableitungen nach „ Mit „=—x und ,=+x 
erhält man 


2 


2,2, = (3.3). 


Für den Fall A=1(U,=0) ergab der Vergleich mit dem Experiment (siehe oben) o = 135 
und also 


Der Überdruck im Strahl ist also sehr klein. Die entsprechende Rechnung, die Tollmien 
nach dem Ansatz (3) durchführte, ergab mit den damaligen Messungen den Faktor 0,0058 
an Stelle von 0,0055. (Es ergibt sich also eine bessere Übereinstimmung, als man er- 
warten darf.) 


Im Falle des ebenen Strahls erhält man entsprechend zunächst enge (nFF'y 
und damit 
1 1 1 "a 
und für 7,=0, 7,=» somit als Überdruck in Strahlmitte gegenüber dem Strahlrand mit 
1,67 
1 oU?(« „oU?(x 


wo U(2)=u(2,0) die Geschwindigkeit auf der Mittellinie y=0 des Strahles bedeutet 
(YzU(«)=ysT,). 

Die Druckdifferenz gegenüber der ruhenden Luft ergibt sich dadurch, daß man den 
Unterdruck am Strahlrand in Größe des Staudrucks der seitlichen Zuströmung in Rechnung 
setzt. Man erhält im Falle des Strahlrandes (4=1) wegen (-&)/[U, = — (20)! FI F,(— x) 


(nach Tollmien 0,0048 statt 0,0046). Im Falle des ebenen Strahles ist r(#,%) 


= — (20)' (x), und man erhält 


Die Druckdifferenzen sind in der Tat so klein, daß ihre Vernachlässigung in den voran- 
gegangenen Rechnungen gerechtfertigt ist. 


4. Probleme des ebenen Nachlaufs. 


Es möge zum Abschluß nach dem neuen Ansatz (6) das Geschwindigkeitsfeld im Nach- 
lauf hinter einem unendlich langen Zylinder und hinter einem aus solchen parallelen und 
äquidistanten Zylindern bestehenden ebenen Gitter in erster Näherung für weit stromabwärts 
gelegene Stellen berechnet werden. Der Druck werde dort wieder als konstant angesehen. 


| 1 
Q 
XUM 
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Ist U, die ungestörte Anströmgeschwindigkeit weit vor dem Zylinder, so setzen wir für 
hinreichend große Entfernungen stromabwärts in erster Näherung 


wo b=b(x) die Nachlaufbreite darstellt. Die Lage der Geraden z<=0 senkrecht zur An- 
strömungsrichtung ist nicht näher festgelegt. Nach dem Impulssatz ist mit Rücksicht darauf, 
daß kein Druckgefälle vorhanden ist und außerhalb des Nachlaufs keine Reibungskräfte 
bestehen 


mit 


b2 
\u,dy=const, 
also nach (4.1), (4.2) i 


Es ist somit für hinreichend große «: 


und man erhält für «,, ®, formal dieselben Differentialgleichungen wie im Falle des lami- 
naren Nachlaufs'*), es tritt nur e, an Stelle der kinematischen Zähigkeit v. Wir können uns 
daher kurzfassen. Da sich aus der Bewegungsgleichung (1.1) n='/» ergibt, berechnet man 
leicht aus den bezüglich «, und v®, linearisierten Differentialgleichungen 

ou, _ ou, 


D, dx . . . . . . . . . . (4.4) 


unter Beachtung der Randbedingungen: u, = für y= +» 


Dabei ist 
(4.6) 


und die Integrationskonstante W stellt den Zylinderwiderstand pro Längeneinheit dar: 


r . 


U; 


(d Durchmesser und e„ Widerstandszahl des Körpers). 


In Bild 9 ist die Geschwindigkeitsverteilung 
dargestellt. Gestrichelt eingetragen ist ferner das 
uü,-Profil nach den Rechnungen von H. Schlich- 

' ting'”) auf Grund des Ansatzes (3), ebenfalls in 
erster Näherung für hinreichend große x. (Die 
beiden «,/|t, \na,-Kurven sind für den Vergleich so 
| gezeichnet, daß sie denselben Flächeninhalt mit der 
3 


| 


waagerechten Achse einschließen, also denselben 
Widerstand liefern.) 

Den Vergleich der Fehlerkurve für «, mit den 

Sachlauf hinter einem Experimenten hat_ bereits H, Reichardt'‘) ge- 

mit dem ”, - Verlauf nach H. Schlichting. geben. Es ist wiederum gute Übereinstimmung bis 

auf die Randgebiete festzustellen. 
Ist y,,, der Abstand von der Mittelachse, in welchem «, /|%, \nax = — 0,5, so ist mit guter 
Näherung nach den Messungen von Reichardt und von Schliehting”) y,/(c„dx)':='. 


| 
N 
2 


14) Zum laminaren Nachlauf siehe S. Goldstein, Proe. Roy. Soc. (A) 142 (1933), S. 545—573. 

15) H. Schliechting: Über das ebene Windsehattenproblem. Ing.-Arch. Bd. 1 (1930), S. 533 bis 571. 
1#) H. Reichardt:a.a.0.S. 14, Bild 1. 

17), Reiehardt: a.a.0. S. 14, Bild 11. 
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Wegen e*=0,5 für £=0,833 ergibt sich somit der Parameter «, zu 


& = 0,0450 W/o U, = 0,0225 c„dU,. 


Im Falle der laminaren Strömung '*) ergeben sich bei der Berechnung höherer Nähe- 
rungen Schwierigkeiten, da man mit einer Entwicklung der Gestalt 


über das Glied f, hinaus nicht zum Ziele kommt. Entsprechende Schwierigkeiten würden 
sich für das von uns behandelte Problem ergeben, in welchem darüber hinaus « nur hin- 
reichend weit draußen konstant ist. 


Im Falle des Nachlaufs hinter einem aus parallelen und äquidistanten zylindrischen 
Stäben bestehenden ebenen Gitter ist 5 hinreichend weit stromabwärts gleich der Gitter- 
weite B zu setzen. Mit 


ergibt sich hier aus der ersten (linearisierten) Gl. (1.1): n=1, und als Differentialgleichung 
für f(n) (mit „= y/B) geschrieben lautet sie 


Somit ist 


Ohne Einschränkung kann y=0 gesetzt werden, und indem wir das Vorzeichen von A positiv an- 
nehmen, haben wir insgesamt nur festgelegt, daß y=0 die Mittellinie zwischen zwei be- 
nachbarten Gitterstäben bilden soll. Die Periode der Lösung in y ist gleich B, also in n 
gleich 1 anzunehmen, somit ist /—=1, d.h. 


Mit Rücksicht darauf, daß fax —fnin =24A ist, besteht der Zusammenhang: 
) 


A=1/8R°%. 


Als Lösung ergibt sich somit: 


U, B_ 2ny U, . 


Nach den Messuugen von R. Gran Olsson’”) und G. Cordes?°) für noch mäßige 
Werte von x/B werden diese Ergebnisse recht gut bestätigt. Der Vergleich ergibt (unter 
Benutzung der Tabelle 5, S. 271 in dem Aufsatz von Gran Ölsson) 


x = (0,103)? 7 — 0,0333. 


Andere Meßergebnisse zeigen systematische Abweichungen des Abklingungsgesetzes von dem 
Potenzgesetz n—=1 und zwar gibt die Auswertung von Messungen Drydens sogar n=4,5°"). 
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18) S, Goldstein: a.a.0O., S. 555. 

19) R. Gran Olsson: Geschwindigkeits- und Temperaturverteilung hinter einem Gitter bei turbulenter 
Strömung. Z. angew. Math. Mech. Bd. 16 (1936), S. 257 bis 274, 
20) G. Cordes: Untersuchungen zur statischen Druck g in turbulenter Strömung. Ing.-Arch. Bd. 8 (197), 
S. 245 bis 270, 


21) Siehe hierüber L. Prandtl: Beitrag zum Turbulenz-Symposium. Vortrag auf dem Internat. Mechanik- 
kongress 1938 in Cambridge (USA). In diesem Zusammenhang verweise ich noch auf die Messungen von J. G. Edler 
von Bohl: Das Verhalten paralleler Luftstrahlen. Ing.-Arch. Bd. 11 (1940), S. 295 bis 314. 
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Der schwimmende Balken. 
Von Friedrich Schiel in Dresden. 


Die Behandlung des Balkens, der auf dem gleichmäßig federnden (= elastisch 
isotropen) Halbraum aufliegt, ist mit großen mathematischen Schwierig- 
keiten verbunden. Wieghardt hat einen Näherungsansatz vorgeschlagen, 
um die Schwierigkeiten zu verringern. Im folgenden wird gezeigt, daß das 
mechanische Modell dieses Ansatzes eine schwere, nichtbenetzende Flüssig- 
keit ist. Die Berechnung eines auf einer solchen Flüssigkeit schwimmenden 
Balkens ist als Näherungsverfahren zur Behandlung ron Bauwerken auf 
nachgiebigem Baugrund geeignet. 


Belastet man einen gewichtslosen, auf einer Flüssigkeit schwimmenden Balken, so wird 
er im allgemeinen auf Biegung beansprucht. Mit diesen Beanspruchungen befassen sich die 
folgenden Darlegungen im Hinblick auf die wertvollen Analogieschlüsse, welche sich auf das 
Verhalten von Balken auf nachgiebiger Unterlage ergeben. Es wird dabei keine 
Vollständigkeit angestrebt, sondern jede Gelegenheit benützt, um durch Näherungsannahmen 
zu rechnerischen Vereinfachungen zu gelangen. 


1. Flüssigkeit ohne Oberflächenspannung. 


Ein elastischer Balken mit rechteckigem Querschnitt von konstanter Breite b soll auf 
einer Flüssigkeit vom Einheitsgewicht y schwimmen. Unter Belastungen, welche in der 
lotrechten Ebene durch die Balkenachse erfolgen sollen, wird er in die Flüssigkeit einsinken 
und sich durchbiegen. Der Auftrieb, mit dem die Flüssigkeit auf die Flächeneinheit der 
Balkenunterseite drückt, ist 


wenn w die Einsenkung gegenüber dem waagerechten Flüssigkeitsspiegel ist. Man sieht, 
daß hier die bekannte Bettungsziffertheorie gilt, welche zur Berechnung von Balken auf nach- 
giebiger Unterlage verwendet wird‘), Dem Einheitsgewicht y der Flüssigkeit entspricht die 
Bettungsziffer, wodurch die Dimension kg/em’ bzw. t/m? dieser Größe ohne weiteres ver- 
ständlich wird. Allerdings gibt es so schwere Flüssigkeiten nicht, deren Auftrieb der Nach- 
giebigkeit z. B. von Baugrund entspricht, zu dessen Kennzeichnung Bettungsziffern von der 
Größenordnung 500 bis 50000 t/m? verwendet werden. 


Die Einsenkungen w gehorchen der Differentialgleichung 


d’ w 
worin E und J die bekannte Bedeutung (Elastizitätsmodul und Trägheitsmoment des Balkens) 
haben und g die Belastung des Balkens bezogen auf seine Längeneinheit ist. 


2. Ansatz mit Oberflächenspannung. 


Wir beschränken uns auf die eindimensionale Aufgabe. Der Balken sei also z. B. ein 
Streifen von der Breite 5, herausgeschnitten aus einem unendlich breiten Balken. Alle Be- 
lastungen sollen längs der Breite gleichmäßig verteilt sein. Die Flüssigkeit soll den Balken 
nicht benetzen, so daß die Flüssigkeitsoberfläche an der Balkenunterseite entlangläuft. 

Der Druck auf die Unterseite des, Balkens besteht dann nicht nur aus dem hydro- 
statischen Anteil ybw, sondern auch aus der Wirkung der Öberflächenspannung o der 
Flüssigkeit. Wäre die Flüssigkeitsoberfläche frei, so bestünde zwischen ihrer Krümmung 
und dem Flüssigkeitsdruck yw die durch den ersten von Laplace aufgestellten Satz der 
Kapillaritätstheorie gegebene Beziehung: 


r, 


worin r, und r, die Hauptkrümmungsradien der Flüssigkeitsoberfläche sind, von denen hier 

wegen der zylindrischen Krümmung der eine unendlich ist. Beschränken wir uns auf kleine Nei- 

gungen, so kann die rechte Seite von (3) o Ts geschrieben werden. Auf die Balkenunterseite 


? 1) Die Übereinstimmung des Verhaltens einer durch die Bettungsziffer gekennzeichneten Unterlage mit dem 
einer Flüssigkeit ist von Schleicher ausgesprochen worden. Beton und Eisen Bd. 26 (1927), S. 434. 
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wirkt der Unterschied zwischen dem hydrostatischen Druck yw und dem Einfluß der Öber- 

flächenspannung, und man erhält somit die Differentialgleichung, welche die Einsenkungen w 

des Balkens beherrscht, wenn man in Gl. (2) das Glied bw durch yb wog ersetzt, 


also 
(4. 


Die Besonderheiten an den Balkenenden werden später besprochen. 

Wir wollen einen anderen Ausgangspunkt für die Behandlung unserer Aufgabe be. 
schreiben, nämlich das Gesetz, nach welchem die Flüssigkeitsoberfläche bei Belastung einer 
Stelle einsinkt. Auf eine unendlich ausgedehnte Flüssigkeitsoberfläche wirke eine gleichmäßig 
verteilte Linienlast. Wir betrachten 


einen Streifen von der Breite b 

\ NY yim nach Bild 1. Auf das Stück b der 

\ Lastlinie entfalle die Last P. Wir 

Bild 1. Einzellast auf einer Flüssigkeitsoberfläche. bezeichnen die Abstände von P 


mit &£ Da wir es mit einer freien Flüssigkeitsoberfläche zu tun haben, gilt nach Laplace 


unmittelbar 
d’w 
| 
Die Randbedingungen für die Lösung bestehen darin, daß w im Unendlichen verschwindet 
und daß im Angriffspunkt die Kraft P mit den beiden Öberflächenkräften ob im Gleich- 
gewicht stehen muß, was sich unter der oben gemachten Voraussetzung kleiner Neigungen 
folgendermaßen ausdrücken läßt: 
de 


20b =P für E= 


Damit ergibt sich schließlich 


Zur Erhöhung der Anschaulichkeit in der weiteren Behandlung führen wir die konstante 
Subtangente s der Exponentialfunktion ein und erhalten 


“ 


) 2bsy 


Die Eigenschaften der Flüssigkeit werden durch die beiden Konstanten y (=Kraft pro 
Volumen) und s (= Länge) gekennzeichnet, welche mit der Oberflächenspannung o in folgendem 
Zusammenhang stehen: 

o=y®#. 

Wir wenden uns jetzt wieder dem Balken zu. Den Druck auf seine Unterseite be- 
zeichnen wir mit p. Wäre p als Funktion von £ gegeben, so könnte man die Einsenkung w 
nach Gl. (6) berechnen, indem man darin P=pd? setzt und die Anteile superponiert. Von 
der Balkenmitte aus zählen wir die Abstände x. Der Balken reicht von —a bis +a. Man 
würde somit‘ w als Funktion von «x erhalten: 


3. Ansatz von Wieghardt. 


Im 2. Jahrgang (1922) dieser Zeitschrift hat Wieghardt (Seite 165 bis 185) den Balken 
auf nachgiebiger Unterlage behandelt. Dieser Aufsatz hat eine gewisse Berühmtheit erlangt, 
weil er lange Zeit?) die einzige genaue Berechnung eines Balkens auf nachgiebiger Unterlage 
enthielt, welche nicht auf der einfachen Bettungsziffertheorie fußt, sondern auch die „Druck- 
ausbreitung* berücksichtigt. 


2) Den unendlich langen Balken auf dem gleichmäßig federnden Halbraum hat Marguerre behandelt, Z. angew. 
Math. Mech. Bd. 17 (1937), S. 224 u. f., den endlichen Balken bzw. den Plattenstreifen Borowicka, Ing.-Arch. Bd. I 
(1939), S. 113 u. f. 
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Die Ergebnisse des Zahlenbeispiels, welches Wieghardt behandelt, wurden wiederholt 
angeführt®), um die Richtigkeit von Näherungsrechnungen für den Balken auf dem gleichmäßig 
federnden Halbraum darzutun. Dessen Eigenschaften bilden aber nicht die Grundlage der 
Untersuchungen von Wieghardt, sondern er ersetzt das von diesem Halbraum bekannte 
Gesetz, wonach die Einsenkungen bei Belastung eines Punktes umgekehrt proportional mit 
der Entfernung abnehmen durch eine Exponentialfunktion, um die Aufgabe einer mathematischen 
Behandlung zugänglich zu machen. 

Wieghardt erhält als Ausgangspunkt seiner Betrachtungen eine Integralgleichung, 
welche mit unserer Gl. (7) völlig übereinstimmt! Man kann also sagen, der Wieghardtsche 
Halbraum ist eine schwere, nicht benetzende Flüssigkeit. Voraussetzung ist dabei, daß die 
Neigungen der Oberfläche klein bleiben, was bei Balken auf nachgiebigem Baugrund immer 
mit genügender Genauigkeit zutrifft. Genauer müßte man bei Wieghardt nicht vom Halb- 
raum, sondern von der Halbebene sprechen, da er nur den Balken, nicht die Platte betrachtet. 

Die Benützung der schweren, nichtbenetzenden Flüssigkeit als Gedankenmodell erlaubt 
nun, die Ergebnisse von Wieghardt zum größten Teil ohne den von ihm benützten mathe- 
matischen Apparat zu erhalten. Die Differentialgleichung (4), welche Wieghardt aus der 
Integralgleichung (7) herleitet, ergab sich z. B. unmittelbar aus den Eigenschaften der Flüssigkeit. 


4. Starrer Balken. 


Die Behandlung des starren Balkens ist mit unseren Hilfsmitteln durchaus kein 
„schwieriges Problem“. Wenn der Balken symmetrisch belastet wird, erfolgt das Ein- 
sinken in die Flüssigkeit nach Bild 2. Wir 
unterscheiden die Einsenkungen innerhalb und 
außerhalb des Balkens durch die Bezeich- 
nungen w; und w„. Die Abstände x zählen 
Bi wir vom Balkenende nach außen. Dann ist 
ild 2. Starrer Balken. 


x 
w;—= konstant = 


Auf die Balkenunterseite wirkt der konstante Druck ybw für die Längeneinheit, und an den 
Enden je eine Kraft D, welche dem Knickwinkel der Flüssigkeitsoberfläche verhältnisgleich ist. 
Bezeichnen wir mit P die gesamte symmetrische Belastung (in Bild 2 durch eine Einzellast 
angedeutet) des Balkens von der Länge I, so erhalten wir: 


m 
D= =. b oder mit 


D=ysbw 
=2D+ymwbl 
1427 


Mit s=0 erhält man aus (8) die Einsenkung nach der Bettungsziffertheorie. 


w= ybl 


5. Vergleich der Flüssigkeit mit dem Halbraum. 


Wenn man die Theorie des schwimmenden Balkens als Ersatz für die Theorie des 
Balkens auf dem gleichmäßig federnden Halbraum benützen will, liegt die Frage nach der 
Größe der Federungszahl (Elastizitätsmodul) des Halbraums nahe, der einer Flüssigkeit 
mit gegebenen Konstanten y und s annähernd entspricht. Wir vergleichen Gl. (8) mit den 
Ergebnissen, die Schleicher‘) bei der Berechnung der mittleren Einsenkung von gleich- 
mäßig belasteten rechteckigen Bereichen des Halbraums erzielt hat. 

Richtiger wäre es natürlich, nicht den starren Balken nach (8) zum Vergleich heranzu- 
ziehen, sondern entsprechend dem Vorgang in der Arbeit von Schleicher einen Balken 
ohne Biegungssteifigkeit, bei dem die Einsenkung der Flüssigkeitsoberfläche bei gleichmäßig 
verteilter Belastung die Form einer Kettenlinie annimmt. Die mittlere Einsenkung ist aber 
nur wenig größer als beim starren Balken, so daß mit Rücksicht auf den Näherungscharakter 
der Rechnung die einfache Gl. (8) verwendet werden soll. 


3) Habel, Näherungsbereehnung des auf dem elastisch-isotropen Halbraum aufliegenden elastischen Balkens. 
Bauingenieur Bd. 19 (1938), S. 76 u. f., insbes. S. 79 u.80; Ohde, Berechnung der Sohldruckverteilung unter Gründungs- 
körpern. Bauingenieur Bd. 23 (1942), S. 99 u. f., insbes. S. 102 und 105. 

4) Bauingenieur Bd. 7 (1926), S. 931 u. f., insbes. S. 949. 
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Nach Schleicher ist die mittlere Einsenkung des Halbraums 


(9), 
worin P die gesamte Belastung, 
F=b.l die Fläche des belasteten Rechtecks, 


m? 
:E wit der in der Elastizitätslehre üblichen Be- 


deutung der Buchstaben und » ein Zahlenwert ist, der sich mit dem Seitenverhältnis des 
Rechtecks ändert. Ist dieses Verhältnis 1:3, 5 und 10, so ist »=0,88, 0,82 und 0,71. 
Setzen wir w=#',„ nach (8) und (9), so erhalten wir 


Je kleiner man die Länge s wählt, desto mehr nähert man sich der Bettungsziffertheorie und 
für s=0 verwandelt sich (10) in jene Gleichung 


welche Schleicher zur Bestimmung der Bettungsziffer empfiehlt, um die Ergebnisse dieser 
Theorie mit der Halbraumtheorie zur möglichst guten Übereinstimmung zu bringen. 

Es liegt die Frage nahe, wie die Konstanten y und s ausgewählt werden müssen, damit 
Balken verschiedener Größe auf derselben nachgiebigen Unterlage möglichst zutreffend ver- 
glichen werden können. Eine eindeutige Antwort ist nur dann möglich, wenn die Unterlage 
der Bettungsziffertheorie gehorcht, nicht aber, wenn sie ein gleichmäßig federnder Halbraum 
ist, weil die Vorschrift (10) hierfür nicht eindeutig ist. Unter den verschiedenen Möglichkeiten 
erscheint mir die folgende Festlegung den Vorzug einer gewissen Natürlichkeit zu haben: 
Man bringe die Konstanten y und s in eine feste Beziehung zur Balkenbreite b, und zwar 


a 


und s=Pßb. 


Will man sich den Verhältnissen des gleichmäßig federnden Halbraums annähern, so muß 
man ß so wählen, daß «a für denselben Halbraum konstant bleibt. Durch Einführung von 


s—=ßb in Gl. (10) ergibt sich 
) b 


und für /=1,5 ist tatsächlich « annähernd konstant für Seitenverhältnisse 1:5=3 bis 10, 
also für den praktisch bei Balken vorkommenden Bereich. Wie man sich durch Einführung 
der &-Werte von Schleicher für diesen Bereich überzeugt, ist nämlich (mit einer Genauig- 
keit von 2°/,) 


Man kann also sagen, daß die Halbraumverhältnisse annähernd durch folgende Festsetzung 
für die Flüssigkeit wiedergegeben werden: 


C - 


Auf die zweckmäßige Auswahl von A und s bei nachgiebigem Baugrund hoffe ich an 
anderer Stelle eingehen zu können. Es kann hier nur der Schluß gezogen werden, daß s 
Werte unter 1,55 annehmen wird, da Baugrund hinsichtlich seiner elastischen Eigenschaften 
zwischen der Halbraum- und der Bettungsziffertheorie steht. Die Angaben von Schleicher‘) 
über C-Werte von Baugrund, welche auf beobachteten Senkungen von Bauwerken und Probe- 
belastungen beruhen, bewegen sich zwischen den Grenzen von 100 und 1000 kgj/em?. Daraus 
ergibt sich nach (11), daß bei einer Balkenbreite 
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b=0,50 m y=0,7 bis 7 kg/em? 
bei b=10 m v—=0,03 bis 0,3 kg/em® 
erwartet werden kann. 


6. Biegsamer Balken mit Gleichlast. 


Zur Berechnung des biegsamen, schwimmenden Balkens steht die Differentialgleichung (4) 
zur Verfügung. Am wichtigsten ist der Fall g=konstant, gleichmäßig verteilte Last, auf 
den sich auch Wieghardt beschränkt. An diesem Belastungsfall wird nämlich die Unzu- 
länglichkeit der Bettungsziffertheorie zur Beschreibung von nachgiebigem Baugrund besonders 
deutlich, da nach dieser Theorie ein gleichmäßiges Einsinken des Balkens olıne Verbiegung 
eintritt, während die Beobachtung von langen Bauwerken lehrt, daß die Mitte stärker einsinkt 
und somit Verbiegungen auftreten. 

Den Weg der direkten Lösung von (4) werde ich nicht beschreiten, obwohl er keine 
besonderen Schwierigkeiten bietet, sondern zeigen, daß es für jenen Bereich der Flüssigkeits- 
konstanten y und s, welcher beim Vergleich mit nachgiebigem Baugrund hauptsächlich vor- 
kommt, eine genügend genaue wesentlich einfachere Näherungsrechnung gibt. Das zweite, 
der Krümmung proportionale Glied von (4) ist nämlich meist im ganzen Verlauf des Balkens 
klein gegen die übrigen, so daß statt (4) die einfachere Differentialgleichung (2) der Bettungs- 
ziffertheorie verwendet werden kann. Nur an den Enden, wo die Krümmung infolge des 
Knickes in der Flüssigkeitsoberfläche unendlich groß ist, spielt das Krümmungsglied eine 
Rolle. Es macht also seinen Einfluß nur auf die Inte- 
grationskonstanten geltend. 

Damit ergibt sich folgende Rechenvorschrift für den 
biegsamen Balken: Außer dem Bettungsdruck p=yır 
wirkt an den Balkenenden je eine Kraft D, welche mit 
der Einsenkung und Verdrehung des Endquerschnittes 
Bild 3. Ende eines biegsamen Balkens. nach Bild 3 folgendermaßen zusammenhängt: 


w 
D=o:b oder mit o=ys? 


worin », die Einsenknng und ww,’ die Verdrehung des Endquerschnittes darstellt. Zum Beweis 

für die Zulässigkeit der Näherung kann man die Balkenbelastung 4 vergleichen mit dem auf 
2 

die Balkenunterseite wirkenden Druckanteil 4 p=-o0gn, der durch die Krümmung er- 

zeugt würde, wobei a über die Größenordnung der Fehler Aufschluß gibt. 

Als Beispiel zur Erläuterung und zum Beweis der genügenden Annäherung in jenem 
Bereich der Konstanten, welche bei Baugrund in Frage kommen, wählen wir den einseitig 
unendlich langen Balken. Die Bettungsziffertheorie liefert’) für eine am Ende senkrecht 
nach oben wirkende Kraft D: 


2D 


dw 


d’w 
dw? ybL? ” 


M=+DrLe sin & (Biegungsmoment) 


Darin ist L eine Vergleichslänge, auf welche die Abstände x bezogen werden: 


(Die Verwendung derselben Bezeichnung £ in (7) und (13) mit verschiedener Bedeutung ergab sich 
aus dem Bestreben, mit den angezogenen Quellen in möglichster Übereinstimmung zu bleiben.) 
Die Gleichlast 4 (auf die Längeneinheit des Balkens) hat nach der Bettungsziffertheorie eine 
gleichmäßige Einsenkung g/yb zur Folge. Demnach ist entsprechend Bild 3 


5) Nach Beyer: Die Statik im Eisenbetonbau, ?. Aufl., Berlin 1933, S. 142. 
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und nach (12) 


zu 
1+2 5; 


Damit ergeben sich Einsenkung und Biegemoment nach (13): 


Darin ist bei w, x nach außen, sonst & nach dem Balkeninneren positiv gezählt. 


Zur Abschätzung der Größenordnung des Fehlers berechnen wir 


d’w 
Ip=- oder mit 0o=ys’ 


esinE. 


1 auf, wo e sin &= 0,3224 ist. 


Ap 120(7,) 


q L 


Der Größtwert tritt bei = 


Wir nehmen den Balkenquerschnitt rechteckig mit der Breite b und der Höhe h an und 
berechnen 


Die Konstanten y und s führen wir nach (11) ein, das bedeutet Annäherung an den 
gleichmäßig federnden Halbraum, und erhalten 


4 


Die Genauigkeit unserer Näherung hängt in entscheidender Weise vom Seitenverhältnis des 
Balkenquerschnittes ab und man erhält das Ergebnis (siehe Zahlentafel 1), daß die Größen- 
ordnung der Fehler bei eigentlichen „Balken“ (b/h < 3) in erträglichen Grenzen bleibt. Die 
Abschätzung der Größenordnung der Fehler in der Zahlentafel gründet sich auf die Annahme 
s=1,b. Ist s kleiner, was bei Baugrund wahrscheinlich ist, so werden die Fehler kleiner. 
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Zahlentafel 1 


C = 1000 kg/em? C = 100 kg/cm? 


7 
1723 L max max — 


! 0.002 0, 0,0001 
Einsenkung 0,28 0,009 0.002 
0,031 : 0.008 
0.097 ‚37 0.027 
0.172 0,055 


Bild 4 (nebenstehend). Einsenkung und Biegemoment eines einseitig 
unendlich langen Balkens mit Gleichlast q auf der Längeneinheit. 


In Bild 4 ist Einsenkung und Biegemoment eines einseitig unendlich langen Balkens 
nach den Gln. (15) dargestellt. Man sieht, daß das Biegemoment in der Nähe des Balkenendes 
seinen Größtwert hat und nach dem Inneren abklingt. Dieses Ergebnis scheint nicht selbst- 
verständlich zu sein, denn für den biegsamen Balken mit Gleichlast ist®) z. B. ein Nähe- 
rungsverfahren vorgeschlagen worden, welches auch Eingang in die Baupraxis fand, nach 
welchem die Biegemomente mit dem Quadrat der Balkenlänge zunehmen, was nur bei kurzen 
Balken zutrifft. 

Ganz entsprechend wie beim einseitig unendlichen Balken läßt sich die Berechnung des 
endlichen Balkens aus der Bettungsziffertheorie mit 
Hinzufügung zweier Kräfte D an den Balkenenden auf- 
bauen. Ich werde an anderer Stelle den Rechnungsgang 
ausführlich behandeln und bringe hier nur zur Ver- 
anschaulichung das Ergebnis (siehe Bild 5) für einige 
Se ee ee gleichmäßig belastete Balken, welche sich nur durch ihre 
verschiedener Länge mit Gleichlast. Länge unterscheiden. 


7. Die Wieghardtschen Ergebnisse. 


Der „Bodendruck* einer durch den Wieghardtschen Ansatz gekennzeichneten Unter- 
lage wird an den Balkenenden nicht nur unendlich, wie beim gleichmäßig federnden Halb- 
raum, sondern er artet zu zwei konzentrierten Einzelkräften aus. Wieghardt ist 
der Meinung, daß die konzentrierten Drücke beim elastischen Balken nicht auftreten dürften, 
sofern man die technische Biegelehre durch eine genauere elastizitätstheoretische Rechnung er- 
setzt. Er findet tatsächlich, daß statt des konzentrierten Druckes ein Unendlichwerden des Boden- 
druckes erfolgt, was aus einem in unendlich kleiner Umgebung erfolgenden Wegdrücken der 
Balkenecke durch die Flüssigkeitsoberfläche verständlich wird. Merkwürdigerweise ergibt 
aber sein Zahlenbeispiel nicht nur für die Balkenenden, sondern auch für den übrigen Verlauf 
merkliche Abweichungen in den Bodendrücken bei genauer elastizitätstheoretischer Berechnung 
gegenüber der technischen Biegelehre. 

Diese Tatsache könnte man nur verstehen, wenn die Unterlage ungewöhnlich hart ist. 
Durch Vergleich unserer Gl. (7) mit den Wieghardtschen Gln. (8) und (9) ergeben sich die 
Konstanten, welche seinem Beispiel zugrunde liegen: 


s=+=32cm 


1 
= Balkenbreite=25 cm, c= 100000 


100000 
2.35.32 


—=62,5 kg/em? 


Vergleichen wir diese Konstanten mit unserer Annahme (11), so ergibt sich zunächst, daß s 
sich von dem für den gleichmässig federnden Halbraum geltenden Wert 1,55 wenig unter- 
scheidet. Benützen wir (10) zur Abschätzung des C-Wertes der Unterlage, so ergibt sich mit 
o=0,71 für ein Seitenverhältnis 10 des Balkens 


625 (14255 AUGE 256 5300 kg/em?. 


6) Kögler-Scheidig: Baugrund und Bauwerk, 2. Aufl., Berlin 1939, S. 213 bis 215. 
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Berücksichtigt man noch, daß Wieghardt mit E=1W06000 kg/em? (z. B. Holz) rechnet, 
während Beton der doppelten Wert hat, so ergibt sich, daß Wieghardt sein Zahlenbeispiel 
mit einer Unterlage berechnet hat, welche im Verhältnis zu seinem Balken rund zehnmal so 
hart ist wie der Boden nach den höchsten Angaben von Schleicher‘) im Verhältnis zu 
Fundamentbeton. 

Über die Integralgleichung (7) und deren Kern bemerkt Nemenyi’) folgendes: „Die 
Auffassung Wieghardts, daß die Unstimmigkeiten, die aus seinem Aufsatz resultieren, auf 
Mängel der Balkenbiegungs-Differentialgleichung zurückzuführen sind, ist nicht stichhaltig. 
Die Unstimmigkeiten rührten zweifellos wesentlich von der Wahl des Kerns her.“ Trotz 
seines (ebenda) erwähnten „sorgfältigen Studiums der Integralgleichungen“ und obwohl er die 
Tragwerke auf elastisch nachgiebiger Unterlage zum Gegenstand einer eigenen Vorlesung ‘) 
gemacht hat, gibt aber Nemenyi weder eine Erklärung für die Wieghardtschen „Un- 
stimmigkeiten“, noch aber hat er die außerordentliche Fruchtbarkeit des Wieghardtschen 
Ansatzes zur Berechnung von Balken auf nachgiebiger Unterlage erkannt. 

Durch meine Darlegungen hoffe ich zur weiteren Anwendungsmöglichkeit der Bettungs- 
ziffertheorie, welche „in neuerer Zeit heftig angegriffen wurde* (Borowicka, a.a.().) beizu- 
tragen, indem ihre Anpassungsfähigkeit an verschiedene Bodeneigenschaften durch zwei 
wählbare Parameter y und s gesteigert wird. 398 


7) Z. angew. Math. Mech. Bd. 11 (1931), S. 461, Fußnote 31. 
8) T. H., Berlin, S. S. 1931. 


Über den Spannungszustand im Halbraum 
bei halbkugelförmiger Druckverteilung. 


Von Nicola Usunoff in Sofia. 


Erster Teil. 


Verfasser berechnet die räumliche Verteilung der senkrechten Spannungs- 
komponente o, und der Schubspannung ı in einem Halbraum, auf dessen 
Oberfläche eine halbkugelförmig verteilte Belastung wirkt. Nach Aufstellung 
und Lösung des Grundintegrals werden die Spannungskurven für eine 
Reihe von Werten der o, und t konstruiert. Merkbare Übereinstimmung 
dieser o,.Kurren bzw. -Flächen mit denjenigen, die von einer gleichwertigen 
Einzelkraft herrühren, ist bei Werten kleiner als 0,03 p, vorhanden. 


Ein elastischer, unendlich großer Körper, der einseitig von einer Ebene begrenzt ist 
(Halbraum), wird durch eine auf einer Kreisscheibe dieser Ebene halbkugelförmig verteilte 
Belastung beansprucht. Es sollen die Spannungskomponenten für jeden im 
Körper beliebig gewählten Aufpunkt ermittelt werden. 

Für die Spannungsberechnung in Punkten, die in genügender Entfernung von der 
belasteten Fläche liegen, kann man die Last durch eine in Mittelpunkt 0 wirkende Einzelkraft 
ersetzen (Bild 1) und die Lösung wird dann nach der einfachen Theorie von Boussinesq 
erhalten. Dagegen muß für Punkte in der Nähe der Druckfläche eine die Grenzbedingungen 
erfüllende Lösung der elastischen Grundgleichungen der allgemeineren Theorie von Hertz 
gesucht werden. Für Aufpunkte in der allernächsten Nachbarschaft der Druckfläche oder auf 
dieser selbst hat Professor L. Föppl, ausgehend von den Boussinesqschen Formeln, 
mittels eines Grenzüberganges gezeigt, daß die Normalspannung gleich dem an der betrachteten 
Stelle herrschenden Druckanteil ist’). Hier soll nun von diesen zwei Extremgebieten ab- 
gesehen und nur das Zwischengebiet behandelt werden, d. h. das Gebiet von einer 
passend gewählten Nachbarschaft der Druckfläche bis zu Entfernungen, bei denen die ver- 
teilte Belastung schon praktisch als Einzellast betrachtet werden darf. 

Bei der kreissymmetrischen Druckverteilung werden die Spannungen, rings um die Lot- 
rechte durch den Mittelpunkt, symmetrisch verlaufen. Daher führt man Zylinderkoordinaten 
ein, wobei die Symmetrieachse mit x, der Radius mit r und der Winkel mit «a bezeichnet seien. 
Dann ist der Spannungszustand im Körper durch folgende vier Komponenten dargestellt: 
o, in der Richtung der Symmetrieachse, o; in der Richtung der Tangente, o, in der des Radius 
und 7 die einzige Komponente der Schubspannung°). Wegen der achsialen Symmetrie genügt 


EN Im Laufe der hochinteressanten Vorträge der Ingenieurkurse 1941 bei der Tee hnischen Hochschule in München 
gab mir Herr Professor L. Föpp!| die Anregung zu dieser Arbeit. Es ist mir eine angenehme Pflicht, ihm an dieser 
Stelle meinen herzlichen Dank auszusprechen. 

2) L. Föppl: Elastische Beanspruchung des Erdbodens unter Fundamenten. Forsch. Ing.-Wes. Bd. 12 (1941), S. 36/37. 


3) Z.B. A. Föppl: Die wichtigsten Lehren der höheren Elastizitätstheorie. Band 5 der Vorlesungen über 
Technische Mechanik. Leipzig 1918, S. 207/210. 
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es, wenn man den Verlauf dieser Spannungen nur in der Ebene eines Meridianschnittes be- 
stimmt, so daß die Frage in zwei Dimensionen zu behandeln ist. In dem vorliegenden 
ersten Teil der Arbeit wird sie nur für o, und r erledigt, die in einem einfachen Ver- 
hältnis zueinander stehen. 


I. Aufstellung des Doppelintegrals für x. 


Wenn eine einzige Last P in einem Punkt auf der Oberfläche des Körpers lotrecht 
wirkt (Bild 1), sind alle Spannungskomponenten im Aufpunkt A durch die Formeln von 
Boussinesq gegeben): 


m u (u + %) u? 
2 | 1 3P 


w 


wobei m die Poissonsche Konstante 
bedeutet. 


Ein beliebiges Flächenelement der 
Druckfläche rdrda (Bild 2) überträgt 
die Belastung prdrda. Sein Beitrag 
zur Druckspannung o, im Aufpunkte A 
ist nach der ersten Formel (1) 


Bild 2 (rechts). 


Bild 1 (unten). 


Sp-rdrda, 


wo von jetzt an anstatt = die Bezeich- 
nung Ah (Höhe des Aufpunktes) ge- 
braucht wird. 


Re Die Entfernung « des Aufpunktes 
vom Lastelement ist 


Wegen der halbkugelförmigen Druckverteilung ist der Druck p an einer beliebigen Stelle 
der Druckfläche durch die Formel 


(4) 


gegeben, wobei p, sein Höchstwert ist. 


Durch Integration über die ganze Druckfläche erhält man den Ausdruck 


d 
00 


(R+r?+h?— 2 Rr cos a)’: 


für die von der ganzen Belastung in dem beliebig gewählten Aufpunkt A hervorgerufene 
Druckspannung. Dies Doppelintegral ergibt o, als Funktion von R und h, d. h. als Orts- 
funktion in der Ebene des Meridianschnittes. Das Minuszeichen für o, in (1) ist, als Zeichen 
für „Druck“ der Einfachheit halber fortgelassen. 


II. Ausführung der ersten Integration. 


Wir führen zuerst die Integration nach a aus, r bleibt konstant und zur Abkürzung 
sei gesetzt: 


2Rr=C...... 


4) Z. B. A. und L. Föppl: Drang und Zwang. Bd. 2, München 1920, S. 226; auch J. Love: A Treatise on the 
mathematical Theory of Elasticity. Cambridge 1920, S. 189. 
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Abgesehen von dem gemeinsamen Faktor, nimmt das Integral über a allein die Form 


A-da da 


(B— cos 


(B — C eosa): 
an. Setzt man weiter 
1 2C B—C 


führt den Halbwinkel = = ein und dividiert Zähler und Nenner des Integranden durch cos° 5, 


so geht das Integral nach einfacher Umrechnung über in 
dB _ 4a 


0% \(a tig: (at 


. 0 
wobei noch 2=tg ß gesetzt worden ist. 


Eine teilweise Integration ergibt sofort ’°) 
|3a(a+ x)": a(a + ) a’(a + 
0 


a 


0 
v 


und durch Teilintegration des ersten Integrals erhält man daraus 


\ \a+tx "a 


0 
0 

(9). 
(1 + 


d 


_ 
a’ 
0 


Das bleibende elliptische Integral ist auf eine Normalform zu bringen. Sein Radikand hat 
die imaginären Wurzeln 


hat also keine mehrfache Nullstelle, ist infolgedessen unreduzierbar. 


) 
Mit a war bereits der Bruch = bezeichnet, dessen Variabilitätsbereich nach (6) 


0<a<l ist. Danach ist der Modul des Integrals 


Nach Ausführung aller Reduktionsprozesse®), nach Einsetzen des Endergebnisses in (9) 
und nach allen Umformungen bleibt auf der rechten Seite eine Summe, bestehend aus den 
elliptischen Normalintegrale erster und zweiter Gattung und dem störenden Glied 


3a’ Vatr 
Dieses verschwindet für =0, während für e=® die Wurzel den Grenzwert Eins, der 
zweite Faktor den Grenzwert Null hat. Infolgedessen hat das ganze Glied den Grenzwert 
Null, so daß Formel (9) jetzt lautet: 
_a+2 
(a+ 
0 
In den Reihendarstellungen der Normalintegrale Y, und Y, verschwinden alle Glieder bis 


0 


(11). 


auf das erste, welches im vorliegenden Falle den Wert end erhält. Die Lösungen dieser 
Integrale sind dann: 


5) Lösungsformeln für die dabei auftretenden einfacheren Integrale findet man z.B. in H. Egerer, Repetitorium 
der höheren Mathematik, München 1908, und E. Brahy, Exereises meöthodiques de Caleul Integral. Paris 1922, 

6) Siehe z. B. Serret-Scheffers: Lehrbuch der Differential- und Integralreehnung. Bd. 2. Leipzig 1911, oder 
auch Jahnke-Emde: Funktionentafeln. Leipzig und Berlin 1933. 
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a2). 


7 1 „2 


Nach Einsetzen in (11) und Zusammenaddieren der Koeffizienten, entsteht aus der rechten 
Seite die Reihe 


a+ta+l 


die verhältnismäßig schnell konvergiert, wenn der Modul &°? nicht nahe bei 1 liegt, das heißt 
nach (10), wenn h groß gegenüber R und r ist. In dem schlimmsten Fall k=1 ist die Reihe 
noch gut brauchbar, wenn es sich nur um Näherungslösungen handelt, denn der Zahlenwert 
des zweiten Gliedes beträgt höchstens ein Achtel von dem des ersten, und der des dritten 
Gliedes nur ein Zwanzigstel. Es würde dann in unserem Fall genügen, wenn wir die Reihe 
schon nach dem zweiten Glied abbrechen. Nach Einsetzen von k? aus (10) und Durch- 
dividieren geht die Lösungsformel (11) in 


über. Damit ist die Integration nach a von O bis 27 erledigt. 


III. Ausführung der zweiten Integration. 


Multipliziert man die rechte Seite von (13) mit dem Faktor (B om aus (8) und setzt 


für A seinen Ausdruck aus (6) ein, so entsteht der Integrand für die Integration nach r 
von O bis d. Nach allen zweckmäßigen Umformungen bekommt das zweite Hauptintegral 
die Gestalt 


d d 


d 


0 


y®—r.rdr OQyd@—r- 
(B—-CO%-yB+C 
0 


wobei nach (6) B und C als gegebene Funktionen der Variable r mit den Parametern R und h 
zu denken sind. 


Wegen r<.d ist die Konvergenz der Binomialreihe 
genügend schnell, so daß man nur die ersten drei Glieder zu nehmen braucht. Das erste 
Integral von (14) zerfällt dann in die Summe 
d 


J 


d d 
rdr 1 | rdr 
(B+OR 8a \(B+ 
0 v0 


Da das vierte Integral noch in der Form 


= "ya: 


geschrieben werden kann, erhält man für die Differenz 3.J, — J,, wie sie in (14) steht, den 
Ausdruck 


d 


_ 
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Setzt man zur Abkürzung: d+R=L, R’+h?’= M und beachtet, d& B+C=r+2Rr+M, 
B-C=r’—2Rr+NM ist, so kann man die obenstehende Integrale durch Substitutionen auf 
bekannte Typen zurückführen und vollständig integrieren’). Man erhält dann die Schlußformel 


Mi 1 1 


31 


(15), 


1 
wobei: 

+6 RW" R’d? 16 Rh? d? — 16 

ht — 25 
R? d* — 64 d? — 16 R’ — 240 R?h* — 4 Rh’ d? — 54 

604° d® + 200. R’ h* — 50 h? R’+ 40 Rh? d? — 56h? R’d? —- 48 Rh? d* 


-W+4d 14 RW de + 35 


1 3R 
—2d’+R’), 


Die bleibenden zwei Integrale J, und J, in (14) sind viel zu kompliziert und würden sicher 
zu einer noch längeren und schwerer brauchbaren Schlußformel als (15) führen. Im folgenden 
werden diese Integrale, für eine Anzahl Werte der Parameter R und h nach der Simpson- 
schen Regel berechnet. Die entsprechenden Werte für (3J, — J,) werden aus (15) ermittelt 
und nach (8) und (5) werden schließlich Zahlenwerte für o. gewonnen, mit welchen man die 
„Gleichspannungskurven“ konstruieren kann. 


IV. Ergebnisse. 

Zur Vereinfachung der Reehnung setzen wir für die Grenze d die Einheit ein, und 
messen r, Rund h nach dieser Einheit ab. Das Intervall von Null bis Eins sei dann in vier 
Teile zerlegt (n—=4), so daß man die Integranden für die Werte '/,, '/,, °/, von r auszurechnen 
hat, da sie für r=0 und r=1 verschwinden. Die Ergebnisse der ganzen Rechnung sind in 
Zahlentafel 1 zusammengefaßt. 


Zahlentafel }t. Für o,=Zahlenwert X ps: 


0 0.089 0,028 
% 061 00 00 0,087 0,06 
% 037 | 0883 0,153 0082 0,05 0.035 0,035 
1 0,2 0,214 | 0,135 | 002 007 0083 0,024 
1% 0.06 0,005 007 0058 0,04 0,029 0,022 
2 0.0053 0,024 | 0048 0.083 0,0336 0.025 0,0192 
3 0,0038 0,0126 0,019 0,0837 0,02 0,016 


Auf Bild 3 sind die entsprechenden Kurven aufgezeichnet, wobei R als Parameter 
gewählt ist. In Bild 4 dagegen ist h als Parameter gewählt. Auf Grund dieser Diagramme 
sind die Gleichspannungskurven für o, im Bild 5 konstruiert worden. Man hat dabei auch 
das unberechnete Gebiet zwischen hR=0 und h=!'/, teilweise ausgefüllt, indem man sich 
darauf gestützt hat, daß jede Kurve in den gleichwertigen Kraftvektor münden muß. Damit 
ist natürlich nur eine grobe Annäherung erreicht und deshalb sind die Kurven dort punktiert 
ausgezogen. 


?) Man benutze dazu etwa die unter 5) zitierten Werke. 


rn 
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T 
| 
| 


Bild 4 


(rechts). 


Bild 6 


(rechts). 


Bild 5 (unten). 


Bild 7 


(unten). 


| 
7 0 


Jeder Kurve in dem Meridianschnitt entspricht im Halbraum eine Rotationsfläche der 
Punkte gleicher Spannung. Man sieht jetzt, daß mit wachsender Entfernung von 0 die 
Spannungsflächen allmählich die Gestalt einer ellipsoidartigen Fläche 
annehmen. 

Wir wollen noch zum Vergleich eine Spannungsfläche angeben, wenn anstatt der ver- 
teilten Belastung eine gleichwertige Einzellast P im Punkte 0 wirkt. So hat man z. B. für 
%=0,03p, nach (1) die Gleichung 


3P h’ 


9 
ıf RA | | 
\ | 
03 
| 
ER | \ | | 
| | | | | 
Bild 3 
| | 
7% h=05 | | | | 
| | 
2 N 1075 AN | | | 
It FINN % 
R 
\\ 
03, 
3 2 2 3 
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Nach Bild 2 und mit Rücksicht auf (4) erhält man zwischen P und p, die Beziehung 


27 1 


22 1 
P=| lag \viorrar= gan. 
00 


0 0 


Durch den Einsatz h=u-sino, R=u:c0os« sei das Polarkoordinatensystem mit Anfang in 0 
und Winkel » zwischen R und « eingeführt, so daß obige Gleichung die Form 


annimmt. 


Die entsprechende o,-Kurve ist im Bild 5 strichweise aufgezeichnet, indem man einige 
Punkte von ihr aus (16) errechnet hat, nämlich 


0° 15° | 30° | | 
|o 20 | 346 | 47 


Bei dem gewählten Wert 0,03p, verlaufen beide Kurven schon ziemlich nahe beiein- 

ander, ausgenommen in der Nähe der Grenzfläche. Bei noch kleinerem Spannungswerte wird 
„der Unterschied zwischen ihnen weiter abnehmen bis bei großer Entfernung die Kurven 

zusammenfallen. 

Von der Fläche 0, = 0,03 p, nach außen würde es sich also kaum lohnen, den Spannungs- 
verlauf auf Grund der verteilten Belastung zu berechnen, denn man würde fast auf die 
gleichen Ergebnisse kommen, die man nach der einfachen Formel von Boussinesq erhält, 
wenn die Belastung als Einzellast in 0 wirkt. 

Zum Schluß wollen wir noch den räumlichen Verlauf der Schubspannung 7 bestimmen. 
Dank dem sehr einfachen Verhältnis, in welchem beide Spannungen o, und r zueinander 
stehen, ist diese Aufgabe durch die bisherigen Rechnungen zum größten Teil schon mitgelöst 
worden. Aus der ersten und vierten Formel (1) folgt direkt (Minuszeichen fortgelassen) 


und danach wurden die r-Werte aus den bereits bekannten #,-Werte errechnet, Zahlentafel 2. 


Zahlentafel 2. Für r= Zahlenwert X ps. 


05 05.05 1 1,5 2 


07 06 0582 0,18, 0,021 
028 | 047 | 0,483 35 | 0,16 
034 0385| 03 | 015 | 004 
021 023 0214 0,131 0,05 | 0,0114 
0,083 0,11 0.08 0,053 0,02 
0.023, 0,042 0,058 30,056 0,05 0,019 
0,01 0.016 0,021 0,024 0.03 0,029 0,019 


0,12 


Die Kurven mit h als Parameter sind im Bild 6 aufgetragen. Daraus entnimmt man die 
Koordinaten R und h der Punkte gleicher Schubspannung und konstruiert die Gleichspannungs- 
kurven für r im Bild 7. 


Eine einfache Überlegung führt zu dem Schluß, daß die r-Kurven bei ihrem weiteren 
Verlauf in der Nachbarschaft der Grenzfläche sich schließen müssen. Mit wachsendem r-Werte 
nehmen die Umfangslängen der Kurven ständig ab, bis an der Stelle der größten Schub- 
spannung die entsprechende Kurve sich auf einen Punkt zusammenzieht. 


Dieser liegt offenbar zwischen dem obersten Kurvenstück (r =0,6 p,) und der Druckfläche, 
links und rechts im Bild 7. Dem entspricht im Halbraum der Kreis der maximalen 
Schubspannung. Um diesen Kreis herum verlaufen alle, der Kurven entsprechenden 
r-Flächen in der Gestalt von Torusartigen Flächen, indem sie sich oberhalb des Kreises 
sehr stark zu der Druckfläche verdichten. Dagegen, sie dehnen sich weiter auseinander bei 
abnehmendem r-Werte unterhalb des Kreises. 


15° 90° 
u 55 | 58 
0,5 0.047 
0.6 
0,75 
1 0 
1.5 
3 
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Wie oben für ox, wollen wir auch hier zum Vergleich einige r-Fiächen ermitteln, die 
von der gleichwertigen Einzellast in 0 herrühren. Z. B. für =0,02p, erhält man nach (1) 


sP #:R 
u? 


0,02 pP, = 


oder in Polarkoordinaten 


und daraus die Zahlenwerte 


15° | 30° | 5" 


0 196 328 42 432 | 347 0 


Wiederholt man die Rechnung noch für 7 = 0,03 p,, so erhält man für « die entsprechenden 
Werte: 


0, 1,6, 2,68, 3,2, 3,53, 2,84, 0. 


Die Kurven sind gestrichelt im Bild 7 aufgetragen. Man sieht, daß sie mit den gleichwertigen 
Kurven von der verteilten Belastung nicht so gut übereinstimmen wie in dem Fall für o,. 
Eine Übereinstimmung im gleichen Maße ist hier bei noch kleinerem Spannungswerte zu 
erwarten. 400 


Harmonische Analyse diskreter Zahlenreihen. 
Von Dr.-Ing. Erich Rud. Berger z. Zt. im Felde. 


Es wird ein Verfahren angegeben, mit Hilfe der bekannten Abminderungs- 
faktoren die Koeffizienten der Trigonometrischen Interpolation instrumentell 
zu bestimmen. Für die Koeffizienten der Interpolation rom Grade k wird 
eine Beziehung aufgestellt, die auf einfache Weise den Abminderungsfaktor t" 
zu berechnen gestattet. Schließlich werden die Formeln angegeben, mit 
deren Hilfe man den Harmonischen Analysator auch zur Harmonischen 
Synthese verwenden kann. 


Wenn man von einer periodischen Funktion nur m einzelne äquidistante Ordinaten 
herausgreift, kann man diese diskreten Funktionswerte in verschiedener Weise interpolieren. 
Man kann sie z. B. durch Stücke von Kurven k-ter Ordnung verbinden, die mit den ersten k — 1 
Ableitungen stetig aneinanderstoßen. Die Koeffizienten der daraus sich ergebenden Fourier- 
Reihe bezeichnen wir. mit a“ und b%. Dieses Verfahren heiße Interpolation k-ter Ordnung. 

Man kann aber auch an Stelle der Integrations-Orthogonalität (in der Fourier-Entwick- 
lung) die Summations-Orthogonalität verwenden und die m Punkte durch einen m-gliedrigen 
trigonometrischen Ausdruck interpolieren, dessen Koeffizienten sich nach den Besselschen 
Formeln berechnen. Quade und Collatz') haben gezeigt, daß diese Trigonometrische Inter- 
polation den Grenzwert der oben beschriebenen Verfahren bildet und daher als „Interpolation 
von der Ordnung unendlich“ zu betrachten ist. Wir wollen diese Koeffizienten daher mit 
a, und b7 bezeichnen. 

Die Trigonometrische Interpolation ist im allgemeinen die vorteilhafteste Methode, da sie 
erstens als einzige einen endlichen Ausdruck liefert (die anderen Methoden liefern unendliche 
Reihen), zweitens die gegebenen Funktionswerte genau wiedergibt (die Reihen konvergieren 
nur dagegen), und da sie drittens die glatteste Kurve liefert; denn die höheren Harmonischen 
fehlen darin vollkommen. Dietsch und Rotzeig?), die dem Geradenzug den Vorzug geben 
wollen, gehen von falschen Voraussetzungen aus. Nach Quade und Collatz') ist jedoch 
dann die Interpolation n-ter Ordnung günstiger, wenn die ersten n—1 Ableitungen der 
gegebenen Funktion stetig sind, die n-te aber unstetig ist, wozu allerdings zu ergänzen wäre: 


wenn die Unstetigkeitsstellen bekannt sind und es möglich ist, sie sämtlich auf Gitterpunkte 
zu legen. 


1) Quade und Collatz: Zur Interpolationstheorie der reellen periodischen Funktionen; Sitzg.-Ber. d. Preuß. 
Akad. d. Wiss., Math.-Naturw. Kl. XXX (1938), S. 383 bis 429. 


. 2) Dietseh und Rotzeig: Eine neue Methode zur exakten Berechnung der Fourier-Koeffizienten: Gerlands 
Beitr. z Geoph. Bd. 38 (1933), S. 276 bis 281. 
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Kennt man die Koeffizienten für eine Art der Interpolation, kann man sie für jede 
andere Art aus diesen berechnen, denn alle diese Koeffizienten sind untereinander proportional: 
die Proportionalitätskonstanten sind die sogenannten Abminderungsfaktoren r: 


RM 


Diese Formeln gelten nur für n< 5 Die höheren Harmonischen hängen mit den niedrigeren 


zusammen durch die Beziehung 

wobei g eine beliebige ganze Zahl ist. Ebenso ist für die sin- Koeffizienten 


Die bekannten rechnerischen Verfahren der harmonischen Analyse (Schema-Verfahren) 
liefern unmittelbar die Koeffizienten a/ und 57. In den oben erwähnten (nur selten vor- 
kommenden) Sonderfällen kann man die Koeffizienten a‘ und b* mit Hilfe der Abminderungs- 
faktoren aus a) und b, berechnen. 

Bei den immer mehr in Aufnahme kommenden instrumentellen Verfahren der 


harmonischen Analyse schlägt man den umgekehrten Weg ein. Man verbindet die gegebenen 
Punkte durch einen gebrochenen Streckenzug und umfährt ihn mit dem harmonischen 


Analysator; das Instrument liefert dann die Koeffizienten a) und b}, aus denen man die 
günstigeren Koeffizienten a7 und b) nach der Formel j 


1 


berechnen kann, wobei r=- gesetzt ist. Es ist aber zu beachten, daß für m=2z und n =z 


denn es ist zwar 


22 22 
1 1 


weil die zweite Summe verschwindet, aber 


22 22 
1 


In analoger Weise erklärt sich das Verschwinden von b?. 
Auf diese Unregelmäßigkeit ist in keiner der zitierten Arbeiten hingewiesen. Sie ergibt 
sich aber ohne weiteres auch aus den Besselschen Formeln. Bei ungeradem m tritt sie 


nicht auf. 


Für n>5 sind alle Koeffizienten gleich Null . . 2 2 2 2 2 2 22... 


Es drängt sich nun die Frage auf, warum diese Koeffizienten verschwinden, die ja im 
Streckenzug alle noch vorhanden waren. Das erklärt sich daraus, daß eine Komponente von 


der Ordnung m — v (0 <v< ”) nicht zu unterscheiden ist von einer gleich starken Kompo- 
nente der Ordnung » mit negativer Phase, weil die Differenz 
cos [(m — + y] — cos (vr y)= — 2sin”, sin 


2 
an allen Gitterpunkten ar verschwindet. Entsprechendes gilt für die höheren Fre- 


quenzen m +» usw. 


; 


A h. . 
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Die verschiedenen Interpolationen unterscheiden sich nun dadurch, daß die höheren 
Frequenzen in immer stärkerem Maße durch niedrigere ersetzt werden, bis zum Schluß die 
Interpolation von der Ordnung unendlich, die Trigonometrische Interpolation, nur mehr die 
niedrigsten Komponenten enthält. Es muß daher die Beziehung gelten 


das ist ein Sonderfall von «a7: 


das ist evident, da a —a# ist, und alle weiteren Glieder versehwinden; für die 
und 
gilt keine derartige Beziehung, da jede dieser Harmonischen an sämtlichen Gitterpunkten 


verschwindet. 


Zusammen mit (1) und (2) ergibt sich daraus 


1 1 1 1 


k k+1 
a, Tu N | 


n 1 
daher mit r = — (r< 3) 
m 


1 1 1 1 


Aus (4b) folgt dieselbe Gleichung. 


2 
Frk=1, d= ergibt das 


1 


Dieselbe Gleichung erhält man auch durch Differenzieren der bekannten Partialbruch-Zerlegung 
für acetgra. Allgemein ist 


1 _ de 


Für die Stufenkurve (k—=0) gilt diese Formel nicht. Hier liegen nämlich die Sprungstellen 
nicht auf den Gitterpunkten, sondern in der Mitte dazwischen. Interpoliert man in dieser 
Art, was für beliebig hohe Grade möglich ist, so ergeben sich Abminderungsfaktoren 7 nach 
der Formel 


dk 


Der Abminderungsfaktor stellt eine Beziehung her zwischen einer Reihensumme (af) 
und einem Integral (a%). Diese Beziehung gilt nicht nur für ganzzahlige », sondern für 
beliebige Werte, da man sie ohne Verwendung der Periodizität durch partielle Integration 
ableiten kann (Sprungstellenmethode®)). Man kann daher Reihensummen bestimmen, indem 


3) Köhler und Walter: Harmonische Analyse von Funktionen mit Sprüngen, Eeken und ähnlichen Besonder- 
heiten. Arch. El. Techn. Bd. 25 (1931), S. 747. 
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man das ihnen entsprechende Integral instrumentell auswertet. Eine praktische Anwendung 
davon wurde soeben beschrieben. Eine zweite ist die Verwendung des Harmonischen Ana- 
Iysators zur Harmonischen Synthese. Es ist nämlich 


= — yesnt=|- 
"m J sınrn) " 
rn 2 
yınnz=i-——) bi. 
"m sinrn) " 


Dabei berechnet sich m aus den Schritten, mit denen « springt: 


_n-4x 
2n 


Man kenne nun von einer Funktion F(n) die Fourier-Koeffizienten A: und Bz. Die 
ersten k von den A: werden über £ aufgetragen (dabei ist 1£—=1), durch einen Streckenzug 
verbunden, der bei A, beginnt und bei A;,+,=0 endet. Dabei ist für A, der doppelte Wert 
einzusetzen, so daß das konstante Glied von F(n) gleich A,/2 ist. Analysiert man nun 
diesen Streckenzug und ergibt sich dabei der cos-Koeffizient «„, so ist nach (9a) der gerade 


- Teil von Fin): 
2 


Ebenso ist nach (9b), wenn man die B: aufträgt und bei der Analyse die sin-Koeffizienten /, 
erhält, der ungerade Teil von F(n): 


Vn)= 


... 


Um dabei eine dichte Wertefolge in » zu erhalten, muß man ein Instrument verwenden, das 
gestattet, die Periodenlänge stetig zu verändern. Von den handelsüblichen Instrumenten 
erfüllt nur der Analysator Mader- Ott diese Bedingung. 


Die Periodenlänge, die man am Instrument einzustellen hat, ist bei Verwendung des 


N-ten Zahnrades P="" N, im £-Maßstab gemessen. Man kann daher nicht in der Um- 
gebung von n=0 arbeiten. Man beschränkt sich zweckmäßig auf das Gebiet 5= nn. 


Für die Abszissen O<n <> formt man um n=7—n, wobei sich die Koeffizienten folgender- 


maßen verändern: 
A:=(—) A: Be=(—)+'Bs 


Man arbeitet dann in der gleichen Weise wie oben, im Abschnitt z <sSn<sa,d.h. 72 zenZV. 
Den Rest der Periode untersucht man im Gebiet - a <n=0; hier gilt die Relation 


Damit sind sämtliche Werte - a<£n=-+7 untersucht. 


Die kürzeste bei diesem Verfahren auftretende Periode ist P\n=2N. Um k Glieder 
zu summieren, muß man den Bereich &=0 bis &=k-+1 analysieren; dazu benötigt man ein 


Zahnrad N> + - 


Die längste vorkommende Periode ist Pya=4N. Der &-Ma&stab ist so zu wählen, 
daß man damit den größten am Instrument einstellbaren Wert nicht überschreitet. 388 


m=— 1=— 
m 
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Über die praktische Integration ebener Vektoren. 
Von 8. v. Borbely in Kolozsvär. 


Es wird ein Verfahren zur zeichnerischen Integration eines von einer 
reellen Veränderlichen abhängigen Integrales einer komplexen Funktion 
gegeben, das die Durchführung der Konstruktion direkt in der komplexen 
Ebene gestattet, ohne dabei durch Komponentenzerlegung auf zwei reelle 
Integrationen zurückgreifen zu müssen. 


Bei manchen Aufgaben der praktischen Mathematik kommt als selbständige Frage- 
stellung — Integration der Hodographenkurve, praktische Auswertung komplexer Linien- 
integrale, Bestimmung der Fourierkoeffizienten usw. — oder als Zwischenkonstruktion — 
z. B. bei der Lösung von Differentialgleichungen in Polarkoordinaten — die Aufgabe vor 
eine, von einer skalaren Veränderlichen f abhängige komplexe Funktion (N) = 
zu integrieren. Dabei ist es von Interesse, die Integration tatsächlich in der komplexen 
Zeichenebene der Funktion w(t) auszuführen, ohne daß man dabei durch Komponenten- 
zerlegung auf zwei reelle Integrationen zurückgreifen müßte. Das konstruktive Gegenstück 
des im Reellen — wegen seiner Einfachheit, und nur durch die Fehler der zeichnerischen 
Durchführung eingeschränkten Genauigkeit — so vielseitig verwendeten Rungeschen Inte- 
grationsverfahrens fehlt im Komplexen. Dies ist an und für sich nicht verwunderlich, 

t 
handelt es sich doch bei der konstruktiven Bestimmung des Integrals RW) ’W—= \w(t)dt 
um die Lösung der folgenden beiden Differentialgleichungen 


arctg (r - 5) =9+y=pll), 


am: (ds: 
(won, 


aus denen die unbekannten Funktionen R(f), O(f) bei praktisch gegebener rechten Seite, ab- 
gesehen von Sonderfällen, exakt nur durch schrittweise Annäherung bestimmt werden können. 

Im Gegensatz zum komplizierten Bau obiger Gleichungen ist ihre komplexe Integral- 
kurve durch sehr einfache Tatsachen gekennzeichnet. 

Die erste stellt eine Richtungsbedingung der Integralkurve dar, daß nämlich ihre Tan- 
gentenrichtung in jedem ihrer Punkte mit der Richtung des entsprechenden Integranden- 
vektors übereinstimmt; die zweite gibt die Parameterbezifferung der Integralkurve in der 
Form an, daß ihre Bogenlänge S gleich dem Integral des Absolutbetrages des Integranden 
ist. Demnach erhalten wir eine erste Näherung der Integralkurve, indem wir obige Bedin- 
gungen längs eines Polygons zu erfüllen trachten. Zu diesegn Zweck teilen wir den Inte- 

t 173 
granden durch die Teilpunkte und bestimmen nach Runge aus S= w| \\w,dt, 
5 to 
dann zeichnen wir das Polygon T,A,T,A,...4; T; A, +1, wobei 4, T, A, + ı| und 
+ gesetzt wird. 

Dieses Polygon erfüllt an den diskreten Stellen T, exakt beide Bedingungen der Integral- 
kurve, die mit den Berührungspunkten T; in sie einbeschriebene Kurve dagegen dort exakt 
nur die Richtungsbedingung, nicht aber die Bedingung der vorgeschriebenen Bogenlänge. 
Wollten wir diese erste Näherung verbessern, so müßten wir die Lage der Punkte T, so ab- 


ändern, daß für die neue Kurve unter Einhaltung der Richtungsbedingung T,n= NS; sein 
i=1 

soll, wobei außer dem Umstand, daß die Abtragung der Bogenlänge praktisch eine sehr un- 

bequeme Zwischenkonstruktion darstellt, auch der Umstand wesentlich ungünstig ist, daß wir 

über die vorzunehmenden Verschiebungen der Berührungspunkte 7; keine Anhaltspunkte be- 

sitzen, bei der Verbesserung also auf ein mehr oder minder geschiektes Probieren an- 

gewiesen sind. 

Selbst wenn wir aus dem Integral durch Reihenentwicklung die zu der richtigen 
Integralkurve zugehörigen Abschnitte des Tangentenpolygons berechnen, erhalten wir nur 
dann konstruktiv einfach zu verwertende Ergebnisse, wenn wir uns mit einer Näherung 
zweiter Ordnung begnügen. 

Wir wollen uns im folgenden immer auf den Winkel g als unabhängige Integrations- 


variable beziehen, und schreiben \ w(d)dt= \ r(g)e'"dy, wobei = gesetzt ist. 
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Wie leicht nachzuprüfen, liefert der Ansatz 


| | 


2 1 


A 4 4A 


Da und 1° von der Größenordnung ist, können wir mit einer 
Genauigkeit bis auf Glieder zweiter Ordnung in 1 schreiben 
kp) +2 1 
Die in den eckigen Klammern stehenden Ausdrücke sind leicht aus PR Kurve des 


Integranden abzugreifen, mit ihnen wird das Tangentenpolygon schrittweise gezeichnet. Die 
Zahlenwerte der Koordinaten der Punkte 7; erhalten wir nach Multiplikation der 
Millimeterkoordinaten dieser Punkte mit dem Maßstabsfaktor #5 3 (1455): 

Obwohl diese Konstruktion schon ziemlich genaue Ergebnisse liefert, befriedigt sie in 
zweifacher Weise nicht. 

Erstens kann es bei dieser Näherung sehr wohl vorkommen, daß die Ungenauigkeit der 
Grundformel größer ausfällt, als es die zeichnerisch erreichbare Genauigkeitsgrenze zuläßt. 

Zweitens sind im allgemeinen die Werte r (,) gegenüber der Korrekturgröße Ar(g;) 
sehr groß. Ist die Korrekturgröße noch bequem abgreifbar, so wächst die Zeichnung schon 
nach einigen Schritten über die Abmessungen unseres Zeichenbogens. Wollen wir das durch 
entsprechende Wahl von x verhindern, so erhalten wir die Korrekturgrößen als kaum auf- 
tragbare Bruchteile unserer Einheit, die jedoch nicht zu vernachlässigen sind, da es auf sie, 
wegen der schrittweisen Fehlersummation, wesentlich ankonmt. 

Aus diesen Gründen versuchen wir die Aufgabe von einer anderen Seite aus, einer be- 
friedigenden Lösung zuzuführen. 

Betrachten wir die Krümmung K der Integralkurve, so ergibt sich unter Berücksichtigung 


r (4 


R 
obiger Differentialgleichungen ünd der Gleichung eotg 


K= 
Die Evolute der Integralkurve ist also dadurch gekennzeichnet, daß sich ihre Banibias als 
Differenz der entsprechenden Absolutbeträge des Integranden darstellen läßt, wobei alle ihre 
Tangenten auf die entsprechenden Vektoren des Integranden senkrecht stehen. 

Ist insbesondere der Integrand ein Kreis um den Nullpunkt des Koordinatensystems, so 
schrumpft die Evolute des Integrals auf einen Punkt zusammen, die Integralkurve wird 
wiederum ein Kreis mit demselben Radius. 

Könnten wir also die Integration der Funktion r(g)ei% abschnittsweise auf die Inte- 
gration von r(®)ei# zurückführen, wobei P=const ein in dem betrachteten Abschnitt 
passend bestimmter Mittelwert der Variablen bedeutet, so hätten wir damit das komplexe 
Gegenstück der Rungeschen Methode gefunden. Es ist leicht einzusehen, daß im Komplexen 
eine solche einfache Abgleichung nicht möglich ist. Betrachten wir nämlich die Abweichung 
ee. abgleichenden Kreises von der gegebenen Funktion, so müßte bei exaktem Abgleichen 

\ @)—r(P)]jerdp=0 gelten. Deuten wir das Integral als Summe der ebenen Zentral- 
kräfte vom absoluten Betrag [r(p) —- r(®)], so kann ® so gewählt werden, daß entweder die 
auf einen passenden mittleren Radius bezogene Summe der Kraftkomponenten, oder die auf 
eine mittlere Tangente bezogene Komponentensumme verschwindet, nicht aber (im allgemeinen) 
gleichzeitig beide zusammen. 
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Da meistens die Komponentensumme in Richtung des mittleren Radius überwiegt, werden 
wir die Abgleichung nach diesem Gesichtspunkt durchführen, es bleibt aber dann als zweite 
Konstruktion noch immer die Erfassung der unberücksichtigten Querkraftkomponenten bestehen. 

Diese beiden Konstruktionen so zu vereinigen, daß sie in einem Arbeitsgang mit prak- 
tischer Genauigkeit durchgeführt werden können, ist unser Ziel. Setzen wir die Reihen- 
entwickung des Integrals in der Form an 

J J 


\ r(g)etdg = \ 
ER 


so ergeben sich für die ersten beiden Koeffizienten 
(0) r’(W)], (0) — 2r’(0) +3ir” (0)]. 


Die Querkraftkomponenten gehen erst in die Glieder dritter und höherer Ordnung der 
Entwicklung durch das Auftreten komplexer Konstanten ein. Die Berücksichtigung dieser 
Glieder ist praktisch notwendig. Sie werden im allgemeinen bei den mir bekannten Kon- 
struktionen zu Unrecht vernachlässigt. Dadurch werden zwar z. B. die Methoden zur Lösung 
von Differentialgleichungen in Polarkoordinaten wesentlich vereinfacht, die so gewonnenen 
Ergebnisse können jedoch numerisch nur als rohe Annäherungen gelten, bei denen auch eine 
evtl. durchgeführte Verbesserung durch wiederholte Annäherungen, wegen dieses Fehlers der 
Grundkonstruktion, nur eine scheinbare ist. 

Eine Darstellung der Korrekturgröße bis einschließlich der Glieder vierter Ordnung in 1 
erhalten wir leicht nach dem Ansatz mit Gewichten 


Ir 
+m, \ e falls m, = 
Es gilt also nach einer einfachen Umreehnung 
Ss 


4? 


6 gesetzt wird. 


J 


+4r(3)+r0) \erap+ \ et dg 


J 


die wir der zeichnerischen Integration im Komplexen zugrunde legen wollen. 

Wir teilen den Integranden durch die Winkel 9,=g,+44. Für jedes dieser Inter- 
valle stellt unsere Formel zwei Integrationen über mittlere Kreissegmente dar, wobei der 
erste Teil die Integration über die Radialkomponenten, der zweite in erster Näherung die 
Korrektur der Querkomponenten liefert. Beschränken wir uns auf Intervallbreiten I=0,1...0,2, 
so liegt die Genauigkeit der Formel selbst bei wiederholter Anwendung noch über der erreich- 
baren Zeichengenauigkeit. 

In Bild 1 sind zwei Schritte der Konstruktion dargestellt. 7, sei der Ausgangspunkt 
der Konstruktion. Wir ziehen durch 7; die Normale N, zu 9, und tragen auf ihr den mitt- 


2 1 
leren Radius R,=r, des Intervalles @,...9;+, ab, wobei die Vor- 
zeichen der Differenzen sorgfältig zu beachten sind. Der so gewonnene Punkt M, ist der 
Krümmungsmittelpunkt der ersten Integralkurve. Durch M, wird die Normale N,+1, zu 9,4, 


gelegt und die Integralkurve dieses Intervalls bis zur Normalen N;+, (Kreisbogen T}; T} + |) 
gezeichnet. Von Tj+, tragen wir auf der Normalen N,;;+, den mittleren Radius R,;, des 
Intervalls ab, durch den Punkt wird gelegt, danach der 


Kreisbogen T} +, T}+» gezeichnet usw. Genau so erfolgt die Konstruktion des zweiten 
(Querkraft-)Integrals, nur liegen dabei die Normalen »;|| zu g,, und die entsprechenden 
Krümmungsradien der einzelnen Integrale sind Die Vektoren T,K;+1, 
geben die Richtung der Verschiebung der Punkte Tj +1, Tj+» infolge Berücksichtigung der 
Querkraftkomponenten an, der Absolutbetrag dieser Verschiebung ergibt sich aus | 7, Ky+ı| 


nach Multiplikation mit dem Maßstabsfaktor z.B: Indem wir an die Punkte Tj +1, diese 
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Bild 1. 


Zahlentafel 1. 


[2 50° 60° 


Näherung [4° 52,70 2 87,84 105,45 
zeichnerisch | 44 52,7 872 | 1046 122,2 


Tabellarisch 


5236 69,81 | 8727 | 10472 122,17 


Fehler +034 +0,45 | +0,57 


Näherung | 4° 14.45 26,74 44,48 52 
zeichnerisch | 44 143 444 69,0 107,0 


Tabellarisch 14,38 | 
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Ar, 
as 77 
Kar N 7, 
A / 
NM. 
/ Wir 
N ; 
/ 
> MIN In 
DAN / N 
% 
15 y 
\ Bild 2. 
ho 
/7 
\ N / Zo 4 
Bild 3. 
\1 77, 
XI" 
; Kos 
70° 
TEE 
| 20 
y 26,65 44,19 69,31 107,29 
+- 0,47 
—0,31 0,29 
XUM 


4.38 Nr. > Okt. 194 Buchholz. Eine Reihentransformation bei einer Fourierschen Reihe 


Korrektur anbringen, erhalten wir die bis auf Glieder vierter Ordnung genauen Integral- 
punkte 7,;+,. Im allgemeinen wird die Abweichung der Punkte Tj}, und 7};, klein sein, 
falls jedoch nötig, legen wir durch T;+, die Tangenten der Integralkurve || zu g,;, und 
zeichnen mit Hilfe des so bestimmten Tangentenpolygons die gesuchte Integralkurve. 

Es ist zweckmäßig auf N, erst |T, @,|=r; abzutragen, die Teilung der entsprechenden 
Differenzen durch Anlegen eines auf Transparentpapier gezeichneten und die Teilung '/s 
resp. °/s tragenden Dreistrahlenwinkels durchzuführen, gleichzeitig die Teilgrößen in den 
Stechzirkel zu summieren und den Punkt M, durch Abtragung dieser Korrektur auf N; fest- 
zulegen. Falls während der Integration es zweckmäßig erscheinen sollte, von der Intervall- 
breite | zu der Intervallbreite 4,— «1 überzugehen, so steht dem bei der Konstruktion der 


ersten Integralkurve nichts im Wege. Zur Wahrung des Maßstabsfaktors =D der zweiten 
Integration müssen wir jedoch von da ab statt x4,,r; «x A, r; als Krümmungsradius abtragen. 


ei 
dpy=x2-+iy mit der Intervallbreite 
4= 0,1745 (10°) für =0,0°...70,0° und der Einheit E=1mm durchgeführt. Gleichzeitig 
wurde diese Aufgabe nach der quadratischen Tangentennäherung in ca. 5fachem Maßstab 
gelöst. Es ergaben sich die in Zahlentafel 1 mitgeteilten Werte. 


In Bild 2 ist die Konstruktion für 1 \ 


iy 
In Bild 3 wurde \ ar mit 1=0,1 im Intervall 9=05...2,0 und 1=0,2 im Inter- 


05 

vll 9=2,0...3,0 mit der Einheit E=200 mm ermittelt. Dabei ergab sich ein mittlerer 
Fehler von ca. 

Zum Schluß sei noch darauf hingewiesen, daß über die instrumentelle Erzeugung der 

Integralkurve auf Grund ihrer Evoluten- und Krümmungseigenschaften demnächst an anderer 

Stelle berichtet wird. 389 


Eine einfache Reihentransformation 


bei einer sehr allgemeinen Fourierschen Reihe. 
Von Herbert Buchholz in Berlin. 


(Mitteilung aus dem Zentrallaboratorium für Fernmeldetechnik der AEG.) 


Für eine Fouriersche Reihe mit dem allgemeinen Glied .cospg-(p a 
wird in der Arbeit eine Reihentransformation durchgeführt, die die Be- 
‚rechnung der Reihensumme besonders dann erleichtert, wenn « nur wenig 
kleiner als 1 oder gleich 1 ist und k eine reelle Zahl, die nicht wesentlich 
qrößer als 1 ist. Die Reihentransformation besteht im wesentlichen in der 
Verwandlung der vorgelegten Reihe in eine Potenzreihe nach x, die in ge- 
wissen Fällen noch ein logarithmenbehaftetes Glied enthalten kann. Die 
Koeffizienten der Reihe bilden in den einfacheren Fällen die Polynome von 
Bernoulli und Euler. Die allgemeine Transformationsformel schließt 
eine große Zahl bekannter Sonderfälle ein. 


1. Allgemeines. 
Bei der mathematischen Behandlung technischer und physikalischer Aufgaben wird man 
häufig auf Reihen der Form (I) oder (II) geführt, in denen die vier Parameter a, r, r und k 
durchweg reelle Größen 


p=1 
mit den Variabilitätsbereichen 


0<jal<1, für k>0 sind. 
Trotzdem es sich wegen cos (7 +.) p=(—)’-cos(px) zunächst zu erübrigen scheint, in den 
obigen Gleichungen nach den beiden Fällen (+r)? und (—r)”’ zu unterscheiden, wird es sich 
doch im Hinblick auf die späteren Ergebnisse empfehlen, hier an dieser Unterscheidung fest- 
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zuhalten. An den Angaben über die Schwankungsbereiche ist bei der Größe a nur wesent- 
lich, daß es nicht gleich einer positiven ganzen Zahl sein darf. Für die Beweisführung könnte 
es sogar komplexe Werte haben. Wir werden jedoch der Einfachheit halber und in Rücksicht 
auf die praktischen Erfordernisse im folgenden bei der Annahme reeller Werte für a bleiben. 

Ist in den beiden Reihen (I) oder (Il) k eine genügend große Zahl und r nicht gerade 
in der Nähe von 1 gelegen, so macht die direkte Berechnung der Reihensumme nach Vor- 
schrift dieser Gleichungen keine besondere Schwierigkeit. Liegt jedoch k etwa unterhalb 
von 5 und ist zudem r nur wenig kleiner oder gar gleich 1, so wird die Berechnung des 
Summenwertes der beiden Reihen schon recht mühselig, und sie wird es um so mehr, je 
kleiner % ausfällt. In solchen Fällen kann die hier aufzustellende Reihentransformation recht 
gute Dienste leisten. Außerdem gestattet sie für eine große Zahl bereits bekannter und auch 
neuer Formeln eine gemeinsame Herleitung. Bevor wir uns jedoch dem Beweise dieser 
Reihentransformation zuwenden, wird es ganz nützlich sein, erst noch auf die besonderen 
Fälle hinzuweisen, in denen man für die obigen Reihen eine geschlossene Summenformel 
schon kennt. 

Rein äußerlich unterscheiden sich diese Summenformeln nach dem Merkmal, daß in der 
einen Gruppe der Ausdruck für die Summe die gleichen periodischen Eigenschaften zeigt wie 
die Reihenglieder selbst, während bei der zweiten Gruppe der Ausdruck für die Reihensumme 
nur für eine bestimmte Periode Anspruch auf Gültigkeit machen kann. An Formeln der ersten 
Art kennt man nur die beiden Formelpaare, die dem Fall k=0 und dem Fall k=1 und a=0 
IFr.cose 


(x, 


+r-sinz 
= r-+ (l.la), (x, r;a)= (1.1b) 


lF2r-cosc+r 

entsprechen. Infolge des periodischen Charakters der rechtsstehenden Ausdrücke gelten diese 
1 

yir2r-cosc-+r” 


+) (+) r- sin. 
(#,r;0)=In 12a), (&,r;0)= + arctg . (1.2b 
( ) ( ) ( ) "(1Fr.cos.r) ( ) 
Formeln für jedes positive und negative reelle x. 
Von der zweiten Gruppe von Summenformeln sind die nachstehenden vier gleichfalls 
seit langem bekannt [1]'). Sie gelten jeweils nur innerhalb des hinter jeder Formel angegebenen 
Periodenintervalls, das in allen Fällen von der Länge 2x ist. Durch fortgesetzte Integration 


SP |< x<27). . .(13a), ST (n<zxz<-+n). . .(1.3b), 


> 


, 


2 .2 


nach x lassen sich aus diesen Bezieliungen beliebig viele weitere Summenformeln herstellen, 
die sich entweder auf die Funktion se, oder auf die Funktion CF beziehen. Alle diese 
Formeln sind auch noch an den Grenzen des Intervalls gültig, wenn die zugehörige Fourier- 
sche Reihe im ganzen abgeschlossenen Intervall gleichmäßig konvergiert, und das ist be- 
kanntlich immer dann der Fall, wenn k>1 ist. Damit hängt es auch zusammen, daß die 
auf der rechten Seite der letzten vier Gleichungen stehenden Ausdrücke für die Randwerte 
des Intervalls bei den beiden Gl. (1.3) verschiedene Werte ergeben, während sie bei den 
Gl. (1.4) zu gleichen Werten führen. Aus der im folgenden hergeleiteten Reihentransformation 
werden sich die Gl. (1.3) und (1.4) als einfache Sonderfälle ergeben. 


2. Die Reihentransformation für die Funktionen und 


Wir gehen zum Zwecke der Durchführung dieser Reihentransformation von der dureh 
Gl. (2.1) definierten Reihe aus. In ihr soll wie vordem 0=Z/al<1, z«=0 und k >0 sein. 


Der Parameter 5 in (2.1) sei vorerst stets positiv und von Null verschieden. Die Reihe (2.1), 


p=1 

die in der angeschriebenen Form auch als eine Dirichletsche Reihe aufgefaßt werden kann, 
ist dann sogar schon für jedes k=0 absolut und gleichmäßig konvergent. Funktionen der 
Art wie E\* ' wurden zum erstenmal von Lerch [2] untersucht, jedoch allein mit dem Ziel, 
für eine solche Funktion hinsichtlich des Parameters k eine Ergänzungsbeziehung zu finden, 
wie sie vordem von Riemann für die speziellere {-Funktion angegeben worden war. 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttum-Verzeichnis am Schluß der Arbeit. 
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Nach dem Beispiel von H. Mellin [3] und [6] ersetzen wir in (2.1) die Exponential- 
funktion durch die Integraldarstellung (2.2), die gültig ist für jedes z mit are 2’ </2, also 


e’= 


sicherlich auch für ein z2=ß —ix mit 5>0. Nach dem Einsetzen von (2.2) in (2.1) ist es 
dann unter den hier vorliegenden Voraussetzungen gewiß erlaubt [4], die Reihenfolge von 
Summation und Integration zu vertauschen, sofern die Integrationsabszisse e in (2.2) so ge- 
wählt wird und so gewählt werden kann, daß Re(k —s)>1 oder e>1--k ausfällt. Das ist 
aber nach (2.2) jederzeit möglich und stellt überhaupt nur für ein k<1 eine über die sowieso 
schon bestehende Vorschrift e>0 hinausgehende Forderung dar. Im Hinblick auf die De- 


» = 1 N j > 


finitionsgleichung (23) der von H. Mellin eingehend betrachteten verallgemeinerten [- 
Funktion [5] ist damit für die Funktion Er von Gl. (2.1) die Integraldarstellung (2.4) ent- 
standen. 
Nach Gl. (2.3) geht für ein a=0 die Funktion Z(—s,1—a) in die gewöhnliche Rie- 


1 & e>— für \ 
\ I (p—ix)‘-ds c>1—k für O<k<if (2. 


mannsche £-Funktion £(—s,1) = über. . So wie diese besitzt auch die Funktion a) 
im endlichen Teil der s-Ebene nur die einzige singuläre Stelle s=1, in deren Umgebung sie 
1 


sich gemäß der Gl. (2.5a) verhält. In den sehr weit entfernten Teilen der s-Ebene regelt 
sich das Verhalten der beiden Funktionen £(s) und £(s,a) überhaupt in merklich der gleichen 
Weise und wird dort durch die Gl. (2.5b) beschrieben, in der s=o-+i-r ist. Die Werte, 


In |! nur dann vorhanden, wenn oder beliebig klein 


die £(s,a) in den Punkten s=+n mit n=0,1,2,... annimmt, sind aus den Gl. (25e...h) 
zu ersehen. Hierin sind die B,„; die Bernoullischen Zahlen, deren Werte in dieser Be- 
zeichnung gleich null zu setzen sind für ungeradzahlige n, die p„(a) die Bernoullischen 
Polynome in der Definition von P. E. Böhmer [7] und [8] und B,(a) dieselben Polynome in 


n 
Pula) +)? . . 
mit 9,@)=2—1; für n22, für 


für jedes s . . . (2.5h) 


der Definition des Tafelwerks von H. T. Davis [9] und = I” (#)/I'(@). Der Zusammen- 
hang zwischen den Polynomen Bernoullis und den Bernoullischen Zahlen geht aus 
Gl. (2.5c) hervor. 


3 
| 
| 
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Nach dieser notwendigen Abschweifung in die Theorie der verallgemeinerten {-Funktion 
untersuchen wir zunächst mit Hilfe der über sie mitgeteilten Formeln das Verhalten des Inte- 
eranden von (2.4) in den sehr weit entfernten Teilen des Integrationsweges sowie der rechts 
davon gelegenen s-Halbebene. Wir schreiben zu diesem Zweck s vorübergehend in der 
Doppelform 

s=R-ei* s—=o-ir 


und haben demgemäß den Näherungswert des Integranden von (2.4) für R>x und für 
'z.=nrl2 zu bestimmen. Dabei wollen wir uns noch der ständigen Abkürzung bedienen 


In Rücksicht auf die aus den Ergänzungssätzen der /. und [-Funktion folgenden Beziehung 


T 
läßt sich dann mit diesen neuen Zeichen nach Gl. (2.5b) dem asymptotischen Näherungs- 
ausdruck für den Integranden von (2.4) die folgende Form geben: 
XeXp R- sin y| (3 sanz) - R-cosy:-In | 
Ebenso wie das Integral (2.2) ist demnach auch das Integral (2.4) auf seinem Integrations- 
weg stets sogar absolut konvergent, solange | < 7/2 und also $>0 ist. Unter der zusätz- 
lichen Bedingung (=?+P°)/2<27 verschwindet aber der Integrand von (2.4) auch auf dem 
unendlich fernen Halbkreis der rechten s-Halbebene. Wird daher der geradlinige, parallel 
zur imaginären Achse verlaufende Integrationsweg von (2.4) durch einen solchen Halbkreis 
zu einem geschlossenen Weg ergänzt, so liefert dieser Teil des Integrationsweges zum Wert 
des Integrals keinen Beitrag. Nun liegt aber rechts vom Integrationsweg von (2.4) einmal 
die unendliche Kette von einfachen Polen s=0,1,2,..., die alle von der /:Funktion her- 
rühren. Sie haben das Residuum 


(-rrin!- 


Nach Gl. (2.5a) liegt aber daselbst auch noch der von der Z- Funktion herstammende einfache 
Pol s=k -1 mit dem Residuum 


— (BP 
Für ein mit n=0,1,2,3, besteht daher für die Funktion E, die unter 
der Voraussetzung (2.8a) absolut konvergente Reihenentwicklung (2.8). Die für die Konvergenz 
9 . . (286) 


des Integrals (2.4) und für seinen Umformungsprozeß noch wesentliche Bedingung (2.Sb) hat 
wegen 7=0 in (2.7) für die Reihe selbst keine Bedeutung mehr. 


(k ==1,2,3,...) 


Liegt der praktisch besonders wichtige Fall vor, daß k einer positiven ganzen Zahl 
gleich ist, so werden auf der rechten Seite von (2.8) sowohl das alleinstehende Glied als 
auch in der Reihe selbst das Glied für n=k- 1 unendlich groß, und zwar mit entgegen- 
gesetztem Vorzeichen. Bestimmt man für k>n-+1 den Grenzwert dieser Differenz auf die 
gewöhnliche Weise oder berechnet man statt dessen direkt das Residuum des Integranden für 
den Doppelpol u=k 1, so entsteht im Falle eines positiv ganzzahligen Wertes von k 
für Er '(#,ß,a) die wiederum nur unter der Bedingung (2.8a) absolut konvergente Ent- 
wieklung (2.9). Der neben dem Strich am Summenzeichen dieser Formel stehende Zeiger 


(k=1,2,3,...) 


| 
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soll dabei darauf hinweisen, daß bei der Bildung der Reihensumme das durch den Zeiger 
gekennzeichnete Reihenglied fortzulassen ist. 

Zerlegt man in der Gl. (2.9) nach Multiplikation mit e-«@#-ix%) beide Seiten in ihre 
reellen und imaginären Bestandteile und setzt sie einander gleich, so kommen unmittelbar 
die gewünschten Transformationsgleichungen zustande. Für ein positiv ganzzahliges k 
nehmen sie die folgende Gestalt an: 


cospx 
(p-a)k 


ea’. e-P? e'? 


n 
‚ eos(ax— nd). (+ -a) 
(2.10), 


k-1 


PA —a) yar+ —9- sin (ax —(k -1)9)) 


9 (k 


sinp« 
(p—a)k 


e-?r— 


sin (ax 


nd). 


k—1 


sin(ax d(k-—1)) a)+In yar +? 
Die Gültigkeit dieser beiden Gleichungen ist in Rücksicht auf die Konvergenz der rechts- 
stehenden Reihen nach wie vor an die Bedingung (2.1la) gebunden, während der durch die 


Gl. (2.11b) definierte Hilfswinkel #9, der in der Integraldarstellung (2.4) seinem Absolutwert 
(B=0) 


nach kleiner als =/2 sein mußte, aus einem schon oben angegebenen Grunde in den obigen 
Reihen für 5>0 auch gleich 1/2 sein darf. Dabei muß nur gehörig auf das Vorzeichen 
von ® geachtet werden. Nach Gl. (2.11b) stimmt dieses Vorzeichen stets mit dem von x 
überein. Daraus geht einmal hervor, daß die Gl. (2.10) und (2.11) auf beiden Seiten in der- 
selben Weise auf einen Vorzeichenwechsel von & reagieren. Diese Gleichungen sind also für 
jedes x gültig, das der Bedingung (2.11a) entspricht, sei es nun positiv oder negativ. 


. (2.1la), 9=aretg (.r|P) . (2.11b) 


Die Gl. (2.11b) zeigt aber auch zugleich, wie im Grenzfall 3>0 mit & und 9 zu ver- 
fahren ist. Offenbar sind deren beide Grenzwerte dann nach den Gl. (2.12) zu bestimmen. 
Dadurch kommt man für #>0 von (2.10) und (2.11) her zu Gleichungen der folgenden Gestalt: 


cospyxX 
Cat; a)— 3 


n=V 


2 


[P(k) — Y(1— a) — In — 


k=1,2,3, 


: 
- 
1839). 
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(x,1;a)= (p—a)* k=1,8%... 


2)". Z(k—n,1—a) 


=sgn#: 


Für k=2 gehört auch, wie man leicht erkennt, der Punkt «=0 zum Gültigkeitsbereich dieser 
Gleichungen. 

Als Beispiel möge für k —=1 unter der wesentlichen Voraussetzung eines a=#=1 mittels 

der Gl. (2.14) der Ausdruck (2.15) gebildet werden, um einmal auch auf diese Weise für ein 


p=1 
dem Intervall O...27 angehöriges x den Wert der in (2.15) auftretenden Summe zu be- 
stimmen. Es gibt zwar, um sie zu ermitteln, einen einfacheren Weg als den hier beschriebenen, 
aber gerade darum bietet sich uns in dieser Aufgabe eine höchst erwünschte Möglichkeit der 
Prüfung. Nach einigen einfachen Zwischenrechnungen ergibt sich zunächst für die Summe von 
(2.15) die Darstellung 


(2.15) 


sin’, 


—sina.«- cos "lea —n,1-a)— 


Nach Gl. (2.5g) hat aber hierin die angeschriebene unendliche Reihe den Summenwert 
— (sina.x)/a, während andererseits die Beziehung gilt: 


... 0.0.0.0. .(2.16). 


Im ganzen kommt somit die Formel (2.17) zustande, die in der Tat zutrifft. 
p-sinpx (a — x) 
2. > 


Ein anderer wichtiger Sonderfall der Gl. (2.10) und (2.11) liegt vor für a=0 und p=V%. 
In der unendlichen Reihe der Gl. (2.13) fallen dann alle in » ungeraden Glieder fort und in 
der Gl. (2.14) alle in » geraden Glieder. Wird dementsprechend n=2m oder 2m +1 ge- 
setzt, so kommen wegen [(s,1)={(s) die beiden folgenden Formeln zustande: 


| | 


-1) 


(2.18), 


(#,1;0) = 


—)” 
k-1 


— \ 
n! 
n—1 
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Beachtet man, daß nach Gl. (2.5c) wegen die Funktion S(—n) für allen=2,4,6,... 
verschwindet, so kann es vorkommen, daß die beiden unendlichen Reihen, die in den obigen 
Gleichungen auftreten, bereits nach einer endlichen Zahl von Gliedern abbrechen, und zwar 
tritt das ein 

bei der Gl. (2.15) für ein geradzahliges k, 

bei der Gl. (2.19) für ein ungeradzahliges k. 
Da gleichzeitig für die so beschaffenen k-Werte auch das mit In |.! behaftete Glied zu Null 
wird, so reduziert sich in den angegebenen Fällen die rechtsstehende Gleichungsseite auf ein 
Polynom. Wir haben damit die Fälle vor uns, die schon im Abschnitt 1 bei den Gl. (1.3) 
und (1.4) erwähnt wurden. In der Tat, setzt man in (2.18) k=2 und in (2.19) k=1, so ent- 
stehen sofort die schon damals angegebenen Beziehungen für C5* (x,1,0) und (x, 1,0). 
Wir stellen diesen längst bekannten Formeln in den Gl. (2.20) und (2.21) die entsprechenden 


1: — Im- (2m)! In sin 2 (2.20), 


p=1 


> sin x 


p=l 


Entwicklungen für die Funktionen (#,1,0) und gegenüber, von denen die 
erste in ihrem ganz rechtsstehenden Summenwert aus “an (1.2a) erschlossen werden kann. 

Schließlich empfiehlt es sich im Hinblick auf die Anwendungen auch noch den Fall 
eines verschwindenden Wertes von x hervorzuheben. Die dafür geltende Beziehung leitet 
man am besten aus Gl. (2.10) her. Sie hat die Form von Gl. (2.22), und für nahe an Null 
gelegene Werte von p, für die die linksstehende Reihe bei kleinen Werten von k sehr schlecht 


— n 


n—0 


konvergiert, genügen für die numerische Berechnung mittels der rechtsstehenden Reihe nur 
wenige Glieder. Für k=0 entsteht aus Gl. (2.22) die bekannte Entwicklung für expie Phe?—1). 

Auch für ein von Null verschiedenes x kann aus den Gl. (2.10) und (2.11) wenigstens 
solange der gleiche Nutzen gezogen werden, als gemäß der Gl. (2.11a) (“? +?) klein bleibt 
gegenüber 27. Selbst für ein B=0 wird dann aber die Konvergenz der eben aufgestellten 
Reihen schon recht mäßig, wenn « in der Nähe von x liegt. Da sich nun aber der bei 
solchen Werten von x ratsame Übergang von x zu z— .r in den Gl. (I) und (II) lediglich 
durch einen Vorzeichenwechsel der in p ungeraden Glieder kundtut, so liegt es nahe, um 
auch für Werte von . in der Nähe von 7 den Vorteil der Reihentransformation wahrnehmen 


zu können, die gleichen Untersuchungen wie an der Funktion A auch an der Funktion u” 
anzustellen. Dabei werden wir uns wesentlich kürzer fassen können. 


3. Die Reihentransformation für die Funktionen ' und 


Wir gehen an die Untersuchung der Funktion Ei? von Gl. (3.1) unter denselben 
Annahmen über den Variabilitätsbereich der vier Parameter x, 5, a und k heran wie bei der 


B—-ix)(n— a) 


(3.1) 


Funktion Et ’ von Gl. (2.1) und ersetzen auch hierin wieder die Exponentialfunktion durch 
die Gl. (2.2). Die Substitution des Integrals in die gegebene Reihe führt dann diesmal nach 
der Vertauschung der Reihenfolge von Summation und Integration auf die durch Gl. (3.2) 
definierte Variante der verallgemeinerten Z-Funktion, die sich in der darin angebenen Art 


| . 22) 
XUM 
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durch die Differenz von zwei der früher benutzten Funktionen mit zwei um den Wert 1/2 von- 
einander verschiedenen Parametern wiedergeben läßt. Die Schreibweise {* für die hier zu 
verwendende [-Funktion gestattet ohne weiteres die rein formale Übertragung der schon im 
Abschnitt 2 gefundenen Beziehungen, so daß wir z, B. sofort für die Funktion Ei? von 
Gl. (3.1) die Integraldarstellung (3.3) anschreiben können. 


Aus der Erklärungsgleichung (3.2) der Funktion {* ergibt sich zugleich, daß die 
Funktion Z*(s,1--«a) in der ganzen endlichen s-Ebene regulär und analytisch oder m. a. W. 
eine ganze transzendente Funktion ist. Sie hat demnach im Gegensatz zur Funktion {(s, 1—«) 
an der Stelle s=1 keinen Pol. Den früheren Gl. (2.5e...h) stellen sich für sie die 
Gl. (34a...e) an die Seite. Hierin hat 9,„(«a) dieselbe Bedeutung wie früher, während E,„_, (a) 


n -&*(l—n,a) =2""' 


1 
E.n-,()=— E,„-, (1) =()"- (2?*—1)- 


*(0,a) =- . (3.4b), . . .(84bN), 


2 


(5 von . . &4e), 


(n=0,1,2,3,..) . . (3.4d), 
*ls,a) s beliebig . . . . . (34e) 


das Eulersche Polynom bedeutet in der Definition des schon erwähnten Tafelwerkes von 
H. T. Davis [9]. 

In den weit entfernten Teilen der s-Ebene zeigt die [*-Funktion nach Gl. (3.2) im wesent- 
lichen das gleiche Verhalten wie die Z-Funktion mit dem einzigen Unterschied, daß wegen 
des Faktors 2° in (3.2) in der entsprechenden asymptotischen Näherungsgleichung (2.7) unter 
dem Zeichen In im Zähler nur der Faktor x steht. Mit dieser einzigen Änderung können 
die früheren Angaben über das Konvergenzverhalten des Integranden von (2.4) unverändert 
auch für den Integranden von (3.3) übernommen werden. Die daher auch hier wieder zu- 
lässige Anwendung des Residuensatzes auf das Integral (3.3) gestaltet sich jedoch im vor- 
liegenden Fall trotz aller dieser Ähnlichkeiten wesentlich einfacher, da es diesmal in jedem 
Falle, es sei nun k ganzzahlig oder nicht, in der rechten s-Halbebene nur einfache Pole gibt. 
Es gilt daher im jetzigen Falle für jedes # die Entwicklung: 


n! 


Nach Trennung der reellen und imaginären Bestandteile dieser Gleichung entstehen daraus 
die folgenden endgültigen Beziehungen: 


cos 


(3.6), 


(3.7). 


- sin (ae 


| | 
n 
= 
1 n! 
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Dabei muß dann im vorliegenden Falle in Rücksicht auf die Konvergenz der rechtsstehenden 
Reihen die Bedingung (3.8a) erfüllt sein, während für #9 wie vordem die Definitions- 
gleichung (3.8b) besteht. 


.... (88a), . . 8b). 
Von den verschiedenen Sonderfällen, die in den Gl. (3.6) und (3.7) enthalten sind, wollen 
wir hier nur auf den Fall eingehen, wo $=0 ist. In der Vorzeichenfrage sind dabei die 


gleichen Vorsichtsmaßnahmen zu beachten wie früher. In einer für positive und negative x 
in gleicher Weise gültigen Schreibart ergibt sich dann aus (3.6) und (3.7) für = 0: 


Sinpx 
Sala) = 2 


. (3.10). 


Die besonderen Formeln, die hieraus für a—=0 entsteben, brauchen wohl nicht angeführt zu 
werden. Wir heben von ihnen nur die beiden Gleichungen heraus, die unter der weiteren 


vereinfachenden Annahme eines k—1 gelten. Sie lauten im Hinblick auf die Gl. (3.4a) und 
(3.4b’) wie folgt: 


> = „+ By 92m er «| 


ST 


81). 

Der zuletzt angeschriebene Summenwert von Reihe (3.11) geht aus Gl. (1.2a) für r—1 
hervor. Seine vorherige Kenntnis würde natürlich in diesem Einzelfall die Herleitung einer 
transformierten Reihe entbehrlich machen. Im übrigen zeigt sich, daß in dieser transformierten 
Reihe hier nichts anderes zu erblicken ist als die im Intervall 03 .r< 7 absolut konvergente 
Taylorsche Reihe für In (cos x/2). Für kleine Werte von x führt sie tatsächlich sehr schnell 
zu einem gut brauchbaren Näherungswert. Für die Reihe (3.12) hat sich durch die Reihen- 
transformation unmittelbar der schon früher angegebene Summenwert von Gl. (1.3b) ergeben. 

Für die numerische Berechnung von Reihen nach Art der Gl. (I) und (Il) bestehen da- 
mit im Falle kleiner oder verschwindender 5-Werte die folgenden beiden Möglichkeiten: für 
Werte von x in der Nähe von 0 bis über 7/2 hinaus erfolgt sie am besten mittels der Gl. (2.10) 
und (2.11) und ihren Abwandlungen. Für Werte von x in der Nähe von 7 werden dagegen 
mit größerem Vorteil die Reihen (3.5) und (3.6) benutzt. Eine praktische Nutzanwendung der 
hier mitgeteilten Formeln wird u. a. in einer Arbeit des Verfassers gemacht, die inzwischen 
im Septemberheft der Zeitschrift für technische Physik erschienen ist. 


4. Schlußbemerkungen. 


Zu der Frage der numerischen Berechnung von Reihen nach Art der Gl. (I) und (II) 
mit Hilfe der im vorstehenden hergeleiteten transformierten Reihen möge noch ergänzend be- 
merkt werden, daß man dabei mit großem Nutzen von dem im Text schon mehrfach erwähnten 
Tafelwerk von H. T. Davis Gebrauch machen kann. In seinem zweiten Bande enthält es 
u. a. auch Tafeln von den Bernoullischen und Eulerschen Polynomen B,„(x) und E,„(x) 
bis zur Ordnung »—=8 für alle x des Intervalls O0...1, und zwar in Schritten von 0,01 mit 
einer Genauigkeit von 10 Dezimalstellen. Außerdem finden sich in dem gleichen Bande auch 


(p-a) 
== = ( —n,1 - a): 
| 
| | 
XUM 
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Tafeln von den Funktionen 7 (»), die z. B. in den Gl. (2.13) und (2.14) vorkommen. Sie 
sind in dem genannten Tafelwerk bis zur Ordnung n—=4 enthalten, und zwar für denselben 
Bereich und für einen teilweise noch feineren Tafelschritt. Schließlich erwähnen wir noch 
die von Gram berechneten und von E. Nörlund herausgegebenen Tafeln der Riemann- 
schen {-Funktion [10], die im Falle eines verschwindenden a benötigt werden und vor allem 
bei nicht ganzzahligen Werten von k von Vorteil sind. 
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Verfahren zur Beurteilung einer sehr geringen Korrelation 
zwischen zwei statistischen Merkmalsreihen. 
Von Hans Gebelein in Eßlingen. 


Während die bekannten Korrelationsmaße in mehr oder weniger befriedigen- 
der Weise dazu dienen können, einen verhältnismäßig engen Zusammen- 
hang zwischen zwei statistischen Merkmalsreihen zu beschreiben, versagen 
sie, wenn es sich darum handelt, den Grad der gegenseitigen Abhängigkeit 
zweier noch eben nicht ganz unabhängiger statistischer Merkmalsreihen zu 
beurteilen, wie dies’namentlich bei der Behandlung biologischer und medi- 
zinischer Fragen erforderlich ist. Hier Klarheit zu schaffen und brauch- 
bare Untersuchungsrerfahren zu gewinnen ist ein Problem, das eng mit 
der Lexisschen Dispersionstheorie zusammenhängt. Im folgenden wird 
jene Theorie so erweitert, daß es gelingt, innerhalb zweidimensionaler 
statistischer Gesamtheiten mit vielfach gegliederten Merkmalen geringfügige 
Korrelationen nachzuweisen. Insbesondere zeigt sich, daß ein bestimmtes 
Vielfaches des Pearsonschen Kontingenzmaßes der dem Problem an besten 
angemessene Parameter ist, indem diese Größe einerseits in den beiden ein- 
fachsten Sonderfällen sich auf die Lexissche Zahl L bzw. auf das Pearson- 
sche 4? reduziert, während anderseits für sie ein universelles, durch Tafel- 
werke der Rechnung zugängliches Verteilungsgesetz von der Art einer wiel- 
dimensionalen Gaufschen Normalverteilung gilt. 


1. Problemstellung. 


Zur Erläuterung des Fragenkreises, der hier behandelt werden soll, betrachten wir zu- 
nächst ein Beispiel aus der immunbiologischen Praxis'). Bei einer Untersuchuug in den 
Städten Breslau und Gottesberg, Kreis Waldenburg, verteilten sich die Kleinkinder von 1 
bis 5 Jahren in dem auf die Impfung gegen Diphtherie unmittelbar folgenden Jahr nach den 
Merkmalen „an Diphtherie erkrankt bzw. nicht erkrankt“ und „geimpft bzw. nicht geimpft“ 
entsprechend den Angaben der Tabellen 1a und ib. Der fragliche Impferfolg muß in einer 
erwünschten Beeinflussung der Verteilung für die Erkrankungen durch die Impfung bestehen; 
daher dürfen die beiden genannten Merkmalsreihen nicht stochastisch unabhängig voneinander 


1) Aus H. v. Schelling: Über die exakte Behandlung des Zusammenhanges zwischen biologischen Merkmals- 
reihen. Arbeit a. d. Staatl, Inst. f. exper. Therapie, Heft 39 (1940), S. 35 bis 71. 
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sein. Wie hoch ist aber diese Abhängigkeit, die auch nicht angenähert eine vollständige ist, 
in jedem der beiden Fälle einzuschätzen, und wie lassen sich weiter die beiden Impferfolge 
miteinander vergleichen ? 
Beispiel 1: Erkrankung von Kleinkindern an Diphtherie im ersten Jahre nach der 
Impfung, beobachtet in Breslau und Gottesberg. 


| nicht nicht 

ı erkrankt erkrankt erkrankt erkrankt 
geimpft. . . . . .| 32264 113 
nicht geimpft... 14813 166 nicht geimpft . 


2 
> 
Tabelle Ia. In Breslau beobachtet. Tabelle Ib. In Gottesberg beobachtet. 


Es liegt nahe, zur Beantwortung dieser Fragen den üblichen Korrelationsk effizienten r 
heranzuziehen, der im vorliegenden einfachsten Fall einer 2X 2-Tafel bekanntlich auch mit 
den beiden Pearsonschen Korrelationsverhältnissen und dem Kontingenzmaß (Mean square 
Contingeney) nach Pearson übereinstimmt. Die Rechnung liefert r—0,0461 für die Bres- 
lauer und r—=0,1115 für die Gottesberger Beobachtungsreihe. Beide Werte sind überraschend 
niedrig, womit zwar nur die Tatsache zum Ausdruck kommt, daß auch olıne Impfung die 
meisten Kinder nicht erkranken, was aber zur Folge hat, daß diese Größen sich als Beweis- 
stücke für den Impferfolg nicht gut eignen. Bedenklicher jedoch ist, daß der Vergleich beider 
Beobachtungsreihen mittels des Korrelationskoeffizienten mißlingt, denn wie sich später heraus- 
stellen wird, liefert trotz des kleineren r die Breslauer Beobachtungsreihe infolge ihres größeren 
Umfangs gegenüber der Gottesberger Beobachtung den stärkeren Beweis für den Erfolg des 
Impfverfahrens. 

Erforderlich ist für die statistische Bearbeitung derartiger Fragen ein Urteil darüber, 
ob die Besetzung der Felder einer Korrelationstabelle, die im allgemeinen Fall k Spalten 
und I! Zeilen aufweisen möge, den Zufallsgesetzen entspricht, oder ob sich in ihr Besonder- 
heiten zeigen, die sich durch Zufall nicht begründen lassen. Als Grundlage für die folgenden 
Untersuchungen wird zunächst mittels der Kombinatorik die Aufteilungswahrscheinlichkeit 
für Matrizen vorgeschriebener Zeilen- und Spaltensummen errechnet. Diese Verteilungsfunktion 
ermöglicht es dann, Kriterien für die Beurteilung auch sehr geringer Korrelationszusammen- 
hänge zu entwickeln. Außerdem treten von dem so gewonnenen Standpunkt aus bemerkens- 
werte Zusammenhänge mit einer Reihe bekannter Verfahren der mathematischen Statistik 
zutage. 


2. Kombinatorische Aufteilungswahrscheinlichkeit für Matrizen. 


Zur Vorbereitung auf das zu lösende kombinatorisch - statistische Problem für Matrizen 
betrachten wir zunächst folgende einfachere Aufgabe °): 

Es liege eine Menge von m Elementen vor, unter denen k Klassen mit den Besetzungs- 
zahlen m,,m,,... m; zu unterscheiden sind, derart, daß 


m, +m,+...+m.=m ist (m; >0, k<m). 


Mit welcher Wahrscheinlichkeit wird bei der Entnahme von Teilmengen, bestehend aus » 
Elementen, eine solche mit den bestimmten Besetzungszahlen »,, n,,... (mit O<£n;<m,, 
erhalten? 

Die Lösung geschieht am einfachsten durch folgende Überlegung: Aus der gegebenen 
m-Menge sind im ganzen n Elemente herauszuheben, was auf (7) Weisen geschehen kann. 
Wir betrachten nun eine bestimmte dieser Teilmengen, z. B. diejenige mit den Besetzungs- 
zahlen n,...n;, und fragen, auf wie viele Weisen sie erhalten werden kann; denn je häufiger 
sie zustande kommt, um so mehr ist mit ihrem Auftreten zu rechnen. Zu ihrem Zustande- 
kommen erforderlich ist erstens die Entnahme von n, Elementen der ersten Klasse in der m- 


m 
Menge, was auf (7) Weisen geschehen kann, ebenso die Entnahme von », Elementen der 
1/ 


zweiten Klasse, wofür es (7) Möglichkeiten gibt usw., bis zur Entnahme von », Elementen 
2 


der k-ten Klasse, was auf wu Arten möglich ist. Daher sind unter den | Pr Teilmengen 


mit n Elementen genau 
n, N, N% 


2) Vgl. H.Gebelein: Zahl und Wirklichkeit — Grundzüge einer mathematischen Statistik. Leipzig 1942, S. 250 ff. 
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solche mit der angenommenen Klassenaufteilung enthalten. Die gesuchte Wahrscheinlichkeit 
ist der Quotient dieser beiden Anzahlen und wir erhalten daher das Frgebnis 


022222... 
(1) 


(m, n+..+nm=n, 0sn,<m; füri=1,2,...k). 


Es ist dies eine arithmetische Wahrscheinlichkeitsverteilung mit (k -—- 1) voneinander 
unabhängigen, ganzzahligen Veränderlichen. Man überzeugt sich leicht, daß die Summe der 
W(n,...nx), erstreckt über alle positiv ganzzahligen Wertesysteme n;, welche die Bedingung 
n,+...+n,=n erfüllen, wie es sein muß, den Wert Eins ergibt. Bemerkenswert ist, daß 
man für die Probefunktionen 


n, mln; —1)(n;—2) usw. 


die Erwartungswerte auf Grund der Verteilung (1) sofort angeben kann. Daher kann man 


anch ohne Schwierigkeit für die relativen Einzelhäufigkeiten ni innerhalb der Teilmengen 


Erwartungswert und Streuung berechnen. Diese wichtigen Ergebnisse lauten 


n m n n(m —1) m m 

Wir wenden uns nunmehr der Frage nach der Aufteilungswahrscheinlichkeit für Matrizen 
zu. Diese Aufgabe lautet: 

Gegeben sei eine Menge von m Elementen, darunter m, solchen der ersten, m, solchen 
der zweiten usw., endlich m; solchen der k-ten Art. Diese Gesamtmenge sei die Zusammen- 
fassung von 1! Teilmengen mit je »; Elementen, in welchen die k Klassen mit bestimmten 
Besetzungszahlen n,;, vertreten sind. Die n;; bilden eine Korrelations- 
tabelle in Form einer Matrix mit k Spalten und ! Zeilen. Die Werte der Spaltensummen 
sind die k Klassenumfänge m; in der Gesamtmenge; die Werte der Zeilensummen die Um- 
fänge der ! Teilmengen, aus denen sich die Gesamtmenge zusammensetzt. Gefragt ist nach 
der Wahrscheinlichkeit für das Auftreten einer bestimmten Matrix (n;;), wenn die Systeme 
der Anzahlen m,... m; und n,...n; vorgeschrieben sind. 


Korrelationstabelle. 


| | k 


) - | 
k Klassen i Sn; 
| | i=1 
@ 
| 
® 
un; m m, ... Mk m 


Nach den getroffenen Annalımen ist für alle @ und j 


k k k 

iwl j=1 i=1 j=1 


Die Frage nach der Aufteilungswahrscheinlichkeit W (n; ;) wird nun gelöst durch wiederholte 
Anwendung des obigen Ergebnisses (1). Indem aus der Gesamtmenge mit m Elementen und 
den Besetzungszahlen m,... m; zunächst eine Teilmenge mit », Elementen herausgehoben 
wird, gilt für die Besetzung der ersten Zeile die Verteilung 


| 


— 
4 
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Für die Entnahme der zweiten Zeile mit n, Elementen steht nun noch eine Menge mit (m —n,) 
Elementen und den Besetzungszahlen (m, —n,,), --- (m; —n;,) zur Verfügung. Daher lautet 
für sie die Verteilung ; 

\ /\ Na N 


Auf diese Weise wird fortgefahren bis zur (I—1)-ten Zeile. Für deren Gewinnung enthält 
die Restmenge nur noch Elemente und zwar je + 
solche der einzelnen Klassen. Daher ist die Verteilungsfunktion für die vorletzte Zeile 


N ,1-ı NK,1-ı 


+ 


Die letzte Zeile endlich ergibt sich zwangsläufig, denn ihre formale Verteilungsfunktion nimmt 
den Wert Eins an. Sie lautet nämlich 


Die Verteilungsfunktion für die Matrix ist nun das Produkt aller dieser Verteilungsfunktionen 


für die einzelnen Zeilen. Beachtet man, daß die an gleicher Stelle stehenden Faktoren sich 
zusammenfassen lassen, indem z. B. 


N, N, 
m! (m—n,)! 


=1, 


ist, so lautet das Ergebnis 


(4). 


Die letztere Schreibweise läßt erkennen, daß eine Vertauschung der Spalten und Zeilen diese 
Aufteilungsfunktion unverändert läßt. 

Für die folgenden Rechnungen sollen nun zunächst einige Erwartungswerte bereitgestellt 
‚werden, nämlich diejenigen der Prüffunktionen 


Das Verfahren ist in all diesen Fällen das gleiche. So folgt z. B. aus der Identität 


(m; —1)!- m; m;;! 


n,!...;— 


wegen der Tatsache, daß rechts bis auf den konstanten Faktor a wieder eine Wahr- 


scheinlichkeitsverteilung nach (4) steht, die über alle Wertesysteme der n;; summiert Eins 
ergibt, der Erwartungswert 
m;n; 
Ein )=- 


und ganz entsprechend lassen sich auch die Erwartungswerte der übrigen oben aufgeführten 
Größen gewinnen. Die Ergebnisse lauten: 


min; 


E (n; =- = — > 


Ein;(m;—D)= 


m; (m; —1) 


m(m—1) ?’ 


m(m—1) 


m! 
k 
m,! m,! m;! 
m!IIn;;! 
N,:...N: 
| 
19 
XUM 
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Es ist nicht schwierig, mittels dieser Formeln Erwartungswert und Streuung für die 


Anzahlen n;; oder auch für die relativen Einzelhäufigkeiten = innerhalb der betreffenden 
er 


Zeile auszurechnen. Es ist entsprechend den Gl. (2) 


oder mit Verwendung der im folgenden häufig gebrauchten Abkürzungen 


3. Die E. Webersche Differenzenprobe. 
Greifen wir zurück auf den Befund der Breslauer Kleinkinder von Beispiel 1. Ohne 


Zweifel kommt der Impferfolg dadurch zum Ausdruck, daß nur in —=0,349 Jo der geimpften 


166 
Kinder erkrankten, während diese Häufigkeit bei den ungeimpften Kindern 14979 = 1,108 %/u 


betrug. Die Differenz 0,759°/o beider Häufigkeiten drückt die Verbesserung aus; um ihre 
Bedeutung einzuschätzen, ist es aber erforderlich, Erwartungswert und Streuung dieser Größe 
zu kennen. Der Vergleich der Abweichung dieser Größe von ihrem Erwartungswert und ihrer 
mittleren quadratischen Abweichung, das ist die Quadratwurzel aus der Streuung, liefert die 
sehr brauchbare Differenzenprobe. 

Wir betrachten den fraglichen Ausschnitt aus der Korrelationstabelle. Die Prüfung er- 
streckt sich nur auf zwei Elemente der Matrix, überdies aus der 
gleichen Spalte. Die Prüffunktion lautet: 


N;jr N;j N; j 5» N; 


d= 


Nach Gl. (6) ist ihr Erwartungswert ERROR 
= m m m; m 
daher kommt es nur auf die Streuung an. Für diese gilt 
) 2 
Str (d= m; _ (m 
nr m N; m 
n» m n; m n» m/\n; m 
wegen Gl. (5) und (6) 
n;)' 
Damit lautet das Ergebnis: 
d= ny E(d=0, Str(d= m(m—]) (8). 


Dieser Ausdruck für die Streuung von d ist bis auf den für die praktische Rechnung 
belanglosen Faktor - = | identisch mit einer Formel, die von Erna Weber für die Diffe- 


renzenprobe im Falle einer 2X 2-Tafel empfohlen worden ist’), und die in biologischen 
" Erna Weber: Einführung in die Variations- und Erblichkeitsstatistik. München 1955, S. 180f. 


. 
- 
| 
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Kreisen viel Beachtung gefunden hat. H. v. Schelling macht in der obengenannten Schrift*®) 
darauf aufmerksam, daß die Anwendung der bekannten Bernoullischen Formel zur Be- 
rechnung der Streuungen von Minuend und Subtrahend und nachfolgende Addition hier nicht 
zum richtigen Ergebnis führt; er setzt sich für die von E. Weber selbst nicht befriedigend 
begründete Formel ein und liefert für sie einen Beweis, dessen Grundgedanke unten im An- 
schluß an Gl. (9) zur Sprache kommen wird. 


Nun kann die in Beispiel la auftretende Differenz von 0,759°/ beurteilt werden. Mit 
41356, m; = 279, 33377, 113, 14979 und n; jr = 166 wird 
279. 47077 
47356 - 47355 \ 14979 
Daher ist die mittlere quadratische Abweichung 
y Str (d) = 7,563 — 0,07563 
und für die hierauf bezogene Differenz ergibt sich 


d 0,7592 
0,0703 


Es kann als ganz ausgeschlossen gelten, daß dieser Wert zufällig zustande gekommen sein 
könnte, denn bekanntlich stellen die Werte +3 bereits ungefähr die Grenze der praktisch in 
Betracht zu ziehenden Fälle dar. Die vorliegende Beobachtungsreihe gibt demnach einen 
starken Hinweis auf den Impferfolg. 

Die Anwendbarkeit der durch Gl. (8) gegebenen Differenzenprobe ist jedoch keineswegs 
auf den einfachen Fall einer 2X 2-Tafel beschränkt, sondern dieses Verfahren kann bei 
irgendwelchen Korrelationstabellen zur Kontrolle zweier etwa bei der Durchsicht auffallender 
Besetzungszahlen aus der gleichen Spalte dienen. 


Str — 57,19. 10°°. 


32377 


Beispiel 2: Differenzen der Augenzahlen bei fünf Wurfserien mit zwei Würfeln. Es 
wurde mit zwei Spielwürfeln gewürfelt und bei jedem Wurf die Differenz der Augenzahlen 
aufgeschrieben. Fünf Wurfserien, die von verschiedenen Personen ausgeführt wurden, lieferten 
die in Zahlentafel 2 zusammengestellten Ergebnisse. 


Zahlentafel 2. 


10 


Es handelt sich hier nicht darum, ob etwa die Würfel „richtige“ sind, sondern die an- 
zustellenden Proben laufen darauf hinaus, eine etwaige systematische Abhängigkeit der Wurf- 
ergebnisse von den ausführenden Spielern festzustellen, die allerdings nicht zu vermuten ist.. 


Besonders stark aus der Reihe fallen z. B. der zweite und der vierte Wert der ersten 
Spalte. Für diese liefert die Differenzenprobe mit »; 7 =10, n,;=11l, n»=55, n;=%, m; =57 
und m =342 nach den Gl. (8) 

10 11 57.25 (1, 1 
d= 55 93 —=0,06225, Str (d) = 32.31 + 5) —=0,04047 
d __0,06225 
Str(d) 0,06361 
Im Unterschied zum vorangehenden Beispiel ist hier die Größe der Differenz demnach durch- 
aus nicht ungewöhnlich; die Werte der beiden Beobachtungen lassen sich zwanglos als zu- 
fällig erklären. 


Str(d)=0,06361 und damit —=0,IT86 . 


4) H. v. Schelling: a.a.0.S. 4f. 


Wie 0 1 2 3 4 Ei 
Serie | 
I 11 u. m 13 6 70 
u 10 15 13 ) 6 55 
u 7 16 9 11 8 55 
IV 11 28 24 22 10 | 92 
18 16 9 | 70 
19* 
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4. Die Schellingsche Trefferprobe. 


Während die Differenzenprobe für die Beurteilung einzelner Felder einer Korrelations- 
tabelle wertvoll ist, besteht auch ein Bedürfnis nach Prüfverfahren, welche auf einmal ein 
Urteil über die ganze Aufnahme zu gewinnen gestatten. Ein bemerkenswertes Verfahren zu 
diesem Zwecke wurde im Jahre 1935 von H. v. Schelling angegeben’). 


H. v. Schelling beschränkt sich dabei auf quadratische Korrelationstabellen, also auf 
den für die Korrelationstheorie besonders wichtigen Fall I=k. Er geht davon aus, daß im 
Grenzfall vollständiger Abhängigkeit der beiden Merkmale nach entsprechender Anordnung 
der Merkmalsreihen nur die Hauptdiagonale von Null verschiedene Besetzungszahlen aufweist; 
und er wählt daher als Prüffunktion die Spur der quadratischen Matrix 


k 


k 
Für deren Erwartungswert folgt nach (Ga) E(s)— m - Zvi 4, für die Streuung aber wird 


Str(s) = mp; 


=> 


Nun ist, ebenfalls nach (6a), 


m’, m? 
und nach (5) 
= Em ) + m’p; pr 


3 2 


Damit wird endlich 


k 


2 k 2 
m 


i=1 
Fassen wir die Ergebnisse zusammen: 
k 
s=- In, Ei=m: 


i i=1 


(9). 


m? 


Die Schellingsche Kennziffer für die Treffsicherheit von Prognosen, von ihm als 7 (Treffer) 
bezeichnet, ist 

s—-E(s) 
VStri(s) ' 


Werte für T bis etwa 3 können als zufallsbedingt gelten, größere Werte geben einen Hin- 
weis auf systematische Zusammenhänge und zwar um so mehr, je größer T ausfällt. 


T— (9a); 


, H. v. Schelling: Die Beurteilung der Treffsicherheit von Prognosen. Dtsch. Statist. Zentralblatt Bd. 37 
(1935), S. 141 bis 148. 
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Beispiel 3: Impfreaktion und Schlußtiter. Bei einer Versuchsreihe von ©. Paschlau 
ergab sich folgender Zusammenhang zwischen der Impfreaktion und dem durch die Impfung 
erreichten Schlußtiter (bei Impfung mit Di-Toxoid Asid)®). 


Zahlentafel 3. 


Impf- Schlußtiter 
reaktion | < 0,005 | bis 0,03. bis 1,0 >10A.E. 


0 0 0 


2 
+. 4 13 13 
2 


++ 
4+++| » 0 12 


m; 20 36 


Die für die Berechnung von T erforderlichen Größen sind in Zahlentafel 4 zusammen- 
gestellt. 


Zahlentafel 4. 


02723 | 011864) 0.613653 0.01902 
0,09091 0,50000 0,04546 0.409090 | .0,01859 
0,22727 0,25000 0,05682 0,52273 0.02970 
0,40909 0,4545 0.025836 


— 0,18905 — 0,09267 


s=3l, E(s)=88- 0,18905 — 16,64, 
2 


Str (s) - (0,18905? + 0,09267) = 11,43, Y'Str(s) 3,381 


und hiermit ergibt sich 


m_31—1664 
T = vn 4,249 . 


Zum Abschluß sei noch erwähnt, daß im Falle k=2 die Maßzahl 7 sich als identisch 
erweist mit dem Ausdruck 781, (d) der Differenzenprobe, dessen Bestandteile nach (8) dann 
mittels der Beziehungen (3) sich ebenfalls noch vereinfachen lassen. Auf dieser merkwürdigen 
Tatsache beruht der Beweis, den v. Schelling neben der unmittelbaren Ausrechnung für 
die E. Webersche Formel gibt”). Der tiefere Grund für diese Übereinstimmung ist der, daß 
es sich im Falle einer 2X 2-Tafel um ein Problem mit nur einem Freiheitsgrad handelt, da 
wegen (3) bei Vorgabe nur einer Feldgröße, z. B. n,, die übrigen Besetzungszahlen sich zwangs- 
läufig ergeben. Daher ist beidemal wegen der Linearität der Prüffunktionen der Summand 
mit », deren allein wesentlicher Bestandteil, und man erhält, wie eine einfache Rechnung 
zeigt, in der Tat die fragliche Formel 


Na 


T=ym-—1:- 
Vn,n,m,m, 


(9b) 


auch für 


Der Faktor von Ym — 1 in der Formel für T ist übrigens genau der Korrelationskoeffizient r. 
Es ist daher nun möglich, auch die zu Anfang für Beispiel la und 1b errechneten Werte 
r=0,0461 bzw. r=0,1115 richtig einzuschätzen. Dem Wert T=1 entspräche ein Korre- 
lationskoeffizient der Größe Va ‚ also gleich 0,00460 im ersteren und 0,0284 im zweiten 
Fall. So gering die Korrelation auch ist, sie ist doch im ersten Fall rund 10-mal und im 


6) Aus H. v. Schelling: Über die exakte Behandlung ... S. 36. 
a.a.0.S. 


nj 
10 
44 
22 | 
| 12 
88 
3 
4 
| 
n.,— E(n,,) 
Bin, 
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zweiten Fall rund 4-mal größer als jene Werte und damit beidemal so hoch, daß kaum zu 
befürchten ist, sie könnte durch Zufallsergebnisse vorgetäuscht sein. 


5. Das universelle Verteilungsgesetz. 


Bevor wir fortfahren, weitere Prüfverfahren zu betrachten, greifen wir zunächst noch- 
mals auf die Verteilungsfunktion Win; nach Gl. (4) zurück, um für sie eine für große m 
asymptotisch geltende Formel zu gewinnen. (Die folgende Umformung entspricht in der Dar- 
stellung und den Bezeichnungen der lehrbue hmäßigen Wiedergabe des Laplaceschen Grenz- 
übergangs in meinem Buche „Zahl und Wirklichkeit“.) 

Der entscheidende Kunstgriff bei der Behandluug der Gl. (4) ist die Einführung neuer 
Veränderlicher an Stelle der Größen n;;. Beachtet man, daß nennenswerte Wahrscheinlich- 
keiten nur für solehe Werte der n;; bestehen, die innerhalb einer gewissen Umgebung von 
deren Erwartungswerten mp;1q; liegen, so ist es naheliegend, als neue Parameter 2;; die 
Unterschiede zwischen »;; und mp; einzuführen. Durch die Substitutionen 


wird mit m; = mp; und n;=mgj 


m;! 


k 
m! n;;! m! TI (mpg+t 2)! 


Wegen der Gl. (3) gelten für die z;; die Beziehungen 


k k 
i=1 jel 


Vor allem ist aber wichtig, daß die interessierenden z;; klein sind im Vergleich zu den 
Größen mp;qj. Aus den Gl. (6a) folgt nämlich, daß die mittleren quadratischen Abweichungen 
wie Ym zunehmen, so daß sie und damit auch die 2;; für einigermaßen großes m gegen- 
über den Größen mp;q; mehr und mehr zurücktreten. Dies hat zur Folge, daß bei hin- 
reichend großem m für alle i,j 
<mpig; 
ist, so daß asymptotisch in m die Formel 
) > 
(mpg 


gilt‘). Hiermit wird nun 


1,7 


k,l 


= 7 ‚mpi q)!-e 

denn wegen (11) ist 
k 
ij 1 ij 


Für die Verteilung W(n;;) lautet das bisherige Ergebnis 


(m pi)! u (m - 


Es ist nun noch die Konstante € zu bestimmen. 
asymptotischen Formeln °) 


Vgl. Zahl und Wirklichkeit, S. 270 


Für deren drei Bestandteile gelten die 


ı Kl zj2 
XUM 
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i=1 


m! ; 


m! 


mi, y2am Pit (pi 


k 


i=1 


Damit wird nach einigen Kürzungen 
— 


(k— 1) 
V Pı- 


Das Endergebnis für die Verteilung W(n;;) aber lautet: 


1 

Dies ist eine Gaußsche Verteilung in einem Raum von (k — 1)(1--1) Dimensionen. Sie 
stellt demnach das universelle Grenzgesetz der Aufteilungswahrscheinlichkeiten für Matrizen 
dar. Bekanntlich stehen für die Arbeit mit dieser Verteilung, namentlich durch die Berech- 
nungen Karl Pearsons ausgedehnte Tafelwerke zur Verfügung. Besonders bemerkenswert 
aber ist, daß die Abhängigkeit von allen Größen »;; mittels einer einzigen Veränderlichen 
erfaßt wird, die somit den dem Problem am besten angemessenen Parameter darstellt. Dieser 
N;j 


Parameter lautet sowie mit Verwendung der Abkürzung h;;= 


k,l 


Er stimmt demnach bis auf einen DES mit dem Pearsonschen Kontingenzmaß 


(hi; "Pi 1 n}j 


Sonderfälle: 
a) Für eine 2X2-Tafel wird f?=r? und daher nach (12) mit k=1=2 


m 


1 1 


Das in diesem einfachsten Fall einparametrige Wahrscheinlichkeitsgesetz der n;; geht dem- 
nach im Grenzfall m>x in eine gewöhnliche Gaußsche Normalverteilung für den Korre- 
lationskoeffizienten r über. Wie man sofort abliest, ist 


(r)= . (12a). 


E(r)=0 und Str (nach (9b) genauer 


b) Bei der bekannten Pearsonschen Verteilung w (2) =(:e ? ” handelt es sich um 
den Schluß von einer sehr großen, nach k Merkmalen gegliederten Ausgangsmenge des Um- 
fangs m auf eine ebensolche, unbegrenzt viel weniger mächtige, aber immer noch sehr große 
Teilmenge vom Umfang n=4g: m. Es liegen also nur !=2 Teilmengen vor mit den relativen 
Umfängen q,=qg und ,=1-—g>-1. Daher vereinfacht sich die Konstante des Verteilungs- 
gesetzes (12) zu 


C=- 


1 1 
das Argument aber wird mit n,,=n,, n,=m; —n, mq,=n und m, =m—.n 


1 
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k 


k 
Pi4: 


(m; — (mn) pi)’ 


(m n)pi 


Dies ist genau die aus der Pearsonschen Theorie bekannte Kennzahl. Es ist demnach (2? 
die Verallgemeinerung der Pearsonschen Prüffunktion y’ auf den Fall einer Matrix mit k 
Spalten und I Zeilen. Die naheliegende Anwendung des Pearsonschen Verfahrens auf 
Korrelationstabellen durch Bildung or Größe 


N (ni; _ 
E(n,;) 


ist hiermit als zulässig erwiesen; zu beachten ist lediglich, daß wegen der Nebenbedingungen (3) 
die Dimensionenzahl der Gaußschen Verteilung nicht gleich der um Eins verringerten An- 
zahl der Veränderlichen ist, sondern (k — 1) (1-1) beträgt”). 

Mit Rücksicht auf eine spätere Anwendung vermerken wir zum Schluß, daß, nachdem 
bekanntlich Erwartungswert und Streuung von z? 

und Str ()=2(k —1) 
ist, diese Größen für (2? die Werte 


besitzen. Bemerkenswert ist endlich noch die für große x, geltende asymptotische Formel für 
die Größe P; (z,), die angibt, welche Wahrscheinlichkeit dafür besteht, daß im bestimmten 
Fall „= z,ist. Sie lautet in der Pearsonschen Theorie '") 


\ 


für sehr großes %,, 


und daher für 2 im Falle einer Matrix mit k Spalten und ? Zeilen 


= 
P;,ı(2,) > — für sehr großes 2, . . (14a). 


1. 
0 


6. Prüfung der Korrelationstabelle mittels der verallgemeinerten Lexisschen Zahl. 


Die Untersuchungen des letzten Abschnitts haben erkennen lassen, daß unter allen mög- 
lichen Prüffunktionen für Kontrollrechnungen an Korrelationstabellen die Größe Q? eine Aus- 
nahmestellung einnimmt. Wir legen daher sie einem Prüfverfahren zugrunde und errechnen 
zu diesem Zwecke nochmals genau für endliches m den Erwartungswert. Es ist nach (13) 


und (6a) 


i=1 i=1 


Dr Vergl. auch R.A. Fisher, Statistical Methods for Research Workers. 7. Aufl. Edinburgh-London 1938, S. 93. 
10) Zahl und Wirklichkeit. S. 282 Gl. (21). 


m 
1 1 
1 
1 
| XUM 
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Dieses Ergebnis 
m 
E (2?) = 


läßt sich zur Normierung der Prüffunktion (2? verwenden. So entsteht mit 
EEE... m—1 Er 
eine statistische Kennzahl, die den Erwartungswert E(L)=1 besitzt. Ihre Streuung exakt 
auszurechnen wäre wegen der vorkommenden Glieder vierten Grades in den n;; recht um- 
ständlich, jedoch folgt für große m aus Gl. (14) hierfür sofort 
.)= Sa) _ 2 
Die Ergebnisse lauten zusammengestellt 
k,l 


Str(L) = Str 


Die Größe L ist eine Verallgemeinerung der bekannten Kennzahl aus der Lexisschen 
Dispersionstheorie, der sog. Lexisschen Zahl. In jener Theorie handelt es sich um unseren 
Fall k=2 bei beliebigem I. In diesem Sonderfall wird mit p, =p, p,=1-—p, n,;j=n,; und 
NN; 


pn; Ad—p)n; 


I 
j=1 


und hieraus durch Division mit E (2?) = 


m—1 
J= 
Dies ist aber die Lexissche Zahl''), womit erwiesen ist, daß unsere Größe L deren Ver- 
allgemeinerung auf Matrizen darstellt. 

Man spricht von übernormaler Dispersion, wenn /, größer als Eins, von unternormaler 
Dispersion, wenn Z, kleiner als Eins ausfällt, wobei aus der Streuung der Bereich folgt, 
innerhalb dessen Abweichungen von Eins auch von Zufallsergebnissen herrühren mögen. Besteht 
zwischen den beiden Merkmalsreihen ein Zusammenhang, so äußert sich dieser durch das 
Auftreten einer übernormalen Dispersion. 

Das Höchstmaß übernormaler Dispersion tritt ein, wenn die einzelnen Teilmengen mit 
je n; Elementen homogen sind. Der- Korrelationszusammenhang wird dann so eng als mög- 
lich, was sich darin äußert, daß im Falle k =! das Kontingenzmaß f? den Wert Eins annimmt, 
während im Falle k==1 zwar f? den Wert Eins nicht erreicht, jedoch die Maximalkorrelation 
K? Eins ergibt'’). Es werde angenommen, die Teilmengen 1 bis l, enthalten nur Elemente 
der ersten Art, die Teilmengen l, +1 bis !, nur solche der zweiten Art usw. Dann sind die 
Teilmengen der Nummern I;_,—+1 bis !; homogen aus Elementen der i-ten Klasse. Daher 
ist für diese n;;—=n; und demzufolge 


k I; I; 
min; min; 
i=1 j=h-ı+tl 
11) Vgl. Zahl und Wirklichkeit. S. 313 Gl. (36). 


12) H. Gebelein: Das statistische Problem der Korrelation als Variations- und Eigenwertproblem und sein Zu- 
sammenhang mit der Ausgleichsrechnung. Z. angew. Math. Mech. Bd. 21 (1941), S. 364 bis 379. 


I; 
m;=1. 
j=h-ı +1 
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Die ganze in (15) vorkommende Doppelsumme aber wird 


—1 m—1 


ij=1 


Naturgemäß ist im angenommenen Fall k <!. Ist umgekehrt k > 1, so vertauschen sich ledig- 
lich die Anzahlen der Zeilen und Spalten. Für die Maximalwerte von Z lautet daher das 
Ergebnis 


Für eine 2X 2- Tafel stimmt, wie auf Grund der Beziehungen (12a) und (9b) leicht einzusehen 
ist, die Kennzahl Z mit dem Quadrat der obigen Kennziffer T von (9a) für die Treffsicherheit 
überein. Dagegen erhält man im allgemeinen verschiedene Ergebnisse bei der Prüfung um- 
fangreicherer Korrelationstabellen. Das L-Verfahren dürfte der Schellingschen T-Probe aus 
zwei Gründen überlegen sein. Erstens beschränkt sich seine Anwendbarkeit nicht auf qua- 
dratische Matrizen; zweitens ist es invariant gegenüber einer Umstellung der Merkmale inner- 
halb der einzelnen Merkmalsreihen. Auch die Tatsache, daß für großes m nicht nur Er- 
wartungswert und Streuung der Größe L bekannt sind, sondern die durch Tabellen erfaßbare 
Wahrscheinlichkeitsverteilung selbst, kann als ein Vorzug gelten. Die mit der L-Probe ver- 
bundene Rechenarbeit ist, wenn man die Formeln (15) benützt, nicht erheblich. 


Zahlenbeispiele: Die beiden statistischen Aufnahmen von Beispiel 1 liefern ohne 
weitere Rechnung wegen L= T? im Falle der Breslauer Versuchsreihe L = 100,8, im Falle 
der Gottesberger Versuchsreihe aber nur 4A= 15,36. Die Größe P nach Gl. (14a) wird daher 
im ersten Fall rund 10°°°, im zweiten etwa 10°*. Man ersieht, wie außerordentlich viel 
beweiskräftiger das Breslauer Ergebnis ist im Vergleich zur anderen Beobachtung, und dies 
trotz des absolut geringeren Korrelationskoeffizienten r. Im Falle des Beispiels 2 mit den 
beiden Würfeln ag die Rechnung mit A —=6 und !=5 und den Zahlenwerten der Zahlentafel 2 


—= 1,0493 und daher mit m = 342 
m; 


i=1 
341 
(1,0193 0,841; YStr(L) =0,316. 


Die sich zeigende geringe unternormale Dispersion kann zu keinen Folgerungen Anlaß geben, 
denn der Wert für Z fällt noch innerhalb der einfachen Streuungsgrenzen 1 + 0,316. Es er- 
gibt sich keinerlei Hinweis auf das Bestehen eines Zusammenhangs zwischen den beiden 
Merkmalsreihen. ae 3 endlich liefert mit k=1=4 und Zahlentafel 3 


Mi —1,474 und daher mit m =88 
‚mn; 


87 
dazu (2) = 0,472. 


Der Wert Z liegt etwas höher als die oben errechnete Größe T—=4,25. Da er weit außer- 
halb der Streuungsgrenzen 1 + 0,472 liegt, besteht an der Tatsache eines Zusammenhangs 
zwischen den beiden Merkmalsreihen kein Zweifel. Hingegen gibt f? als Korrelationsmaß 
mit seinem Wert 


(1,474 1) — 0,158 


keinen anschaulichen Begriff von der Enge und vor allem keinen zutreffenden Begriff von 
der Bedeutung dieses Korrelationszusammenhangs. 


Während bei den üblichen Korrelationsmaßen ein Vergleich mit dem nichtstatistischen 
Grenzfall des vollständigen funktionsmäßigen Zusammenhangs zwischen den beiden Merkmals- 
reihen angestellt wird, um so ein Urteil über die Enge der Verknüpfung zu gewinnen, ge- 
schieht bei den hier beschriebenen Prüfverfahren ein Vergleich mit dem entgegengesetzten 
Grenzfall vollständiger Unabhängigkeit beider Merkmalsreihen, der selbst ein statistisches 
Problem darstellt. Das Ergebnis dieses Vergleichs ist ein Urteil über die Bedeutung des 
zutage tretenden Korrelationszusammenhangs. 406 


k,l 
| 
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BUCHBESPRECHUNGEN 


Dr.-Ing. WALTER TRAUPEL, Neue allge- 
meine Theorie der mehrstufizenaxi- 
alen Turbomaschine. 150 S. m. 32 Fig. u. 
i Tafeln. Zürich und Leipzig 1942. Verlag Gebr. 
Leemann & Co. Preis brosch. 7.20 M. 

Der Verfasser hat in dem wertvollen Büchlein 
mit Erfolg versucht. das weit verstreute und von 
sehr verschiedenen Gesichtspunkten ausgehende 
Schrifttum in einen einheitlichen Gedankengang 
zusammenzufassen. Darüber hinaus liefert er in 
der Arbeit für die Theorie der Strömungsmaschine 
verschiedene wichtige Beiträge. Das Buch enthält 


zunächst eine Ableitung und Zusammenstellung der 


für Turbine und Verdichter gemeinsam geltenden 
Grundgleichungen, die sich aus der Theorie des 
mittleren Stromfadens ergeben. Anschließend be- 
handelt es die bei den Strömungsversuchen an 
Schaufelgittern auftretenden grundsätzlichen Fra- 
gen, vor allem in bezug auf die Zulässigkeit der 
meistens gemachten Vernachlässigungen bei den 
Messungen. Die Ergebnisse der Gitterversuche 
werden in verschieden definierten Verlustkenn- 
zahlen zusammengefaßt und deren Brauchbarkeit 
bei der Berechnung von. einzelnen Stufen und 
zanzen Strömungsmaschinen auf Grund der Gitter- 
versuche nachgewiesen. Den Abschluß bildet eine 
Theorie -der räumlichen Strömung zähigrkeitsfreier 
Flüssirkeiten in der Turbomaschine. 

Die meisten Ableitungen und Betrachtungen (des 
Buches beziehen sich auf die unzusammendrück- 
bare Flüssigkeit. Es werden aber an mehreren 
Stellen Hinweise darauf gegeben. inwiefern sich 
die zusammendrückbare Flüssigkeit unterscheidet. 
Verschiedentlich bringt der Verfasser durchaus 
neuartige Gedanken und Darstellungsweisen, die 
sich für die Theorie als sehr fruchtbar erweisen. 
Immerhin dürften nicht alle Ausführungen in 
diesen Abschnitten vom Physiker oder Mathema- 
tiker ohne gewisse Bedenken hingenommen werden 
können. Das ist aber bei allen Arbeiten so. in 
denen neue Wege gezeigt werden. 

Es ist kaum anzunehmen. daß das Buch die 
zur Zeit im Gebrauch befindlichen Berechnungsver- 
fahren von Strömungsmaschinen wesentlich ab- 
ändern wird, trotzdem diese Verfahren dringend 
einer Erneuerung bedürfen. Auch die weitzehend 
verfeinerten in dem Buch dargestellten Theorien 
können nämlich nicht zu einer wirklich sicheren 
Vorausberechnung von Beschauflungen führen. 
Der Verfasser weist selbst mehrfach ausdrücklich 
auf die Schranken hin, die der Theorie und auch 
unseren Versuchsverfahren bis heute anhaften. 
Trotzdem sollte das Buch überall dort gelesen 
werden. wo man sich um das tiefere Eindringen 
in die Vorgänge innerhalb der Strömungsmaschine 
bemüht. Hier kann es fruchtbare Dienste leisten. 
die sich auch für den Strömungsmaschinenbau 
auswirken werden. 

Augsburg. E. Sörensen. 415 

Dr. phil. HERMANN ATHEN, Ballistik. 
(Hochschulwissen in Einzeldarstellungen.) 298 S. 
m. 40 Abb. Leipzig 1941, Verlag Quelle & Meyer. 
Preis geb. 8.60 M. 


Das Buch ist in 5 Hauptteile gegliedert: 1. Äußere 
Ballistik, 2. Anwendung der Wahrscheinlichkeits- 


lehre auf die B., 3. Schußtafel-B., 4. 
5. die B. des Bombenwurfs. 

Es bringt neben den alten bekannten Verfahren. 
(die für besondere Zwecke noch immer ihre Bedeu- 
tung in der Ballistik behauptet haben, die Ergeb- 
nisse der letzten Jahrzehnte (die in dem grund- 
legenden Werk’ von Cranz nicht mehr berücksich- 
tiet werden konnten). 


innere B.. 


Besonders hervorzuheben sind neuere Betrach- 
tungen über den Luftwiderstand. Sätze über ähn- 
liche Flugbahnen, Parameterentwicklungen für 
Flugbahnscharen, der Abschnitt über photogramme- 
tiische Flugbahnvermessung und die Behandlung 
(er Störungsrechnung. 

Einige Dinge sind dem Ref. bei der Durehsicht 
aufgefallen: Bei der Herleitung des Stabilitätsfak- 
tors von Geschossen darf die Querkraft nieht ver- 
nachlässiget werden. Ein „amerikanisches“ Luft- 
widerstandsgesetz gibt es nicht. Den Exterior 
Ballistie Tables liegt (bis auf einen Maßstabsfaktor) 
das Gesetz von Gävre zugrunde. — Das „Normal- 
Iuftgewieht” bei Heer und Marine ist verschieden. 
1.22 geren 1.245 ke/m?. — Bei den Betrachtungen 
über den Luftwiderstand sollten die Ergebnisse der 
Aerodynamik ausführlicher erwähnt werden. 

Die innere Ballistik tritt im Vergleich zur äuße- 
ren etwas zurück. 

Ein ausführlicher Schrifttumnachweis und eine 
Zusammenstellung wichtiger mathematischer Sätze 
und Formeln ergänzt das Buch. 

Das Werk vermittelt trotz teilweise gedrängter 
Fassunz in klarer Weise die theoretischen Grund- 
lagen der Ballistik. wobei jedoch vielfach Brücken 
nach der praktischen Seite hin geschlagen werden 
im Gegensatz etwa zu dem Werk von Popoff).. Der 
Verfasser hat sein Ziel, die wissenschaftlichen 
Grundlagen der Wehrtechnik zu vertiefen und zu 
ihrer Verbreitung beizutragen, voll erreicht. Da der 
Inhalt des Buches zum unentbehrlichen Rüstzeug 
jedes Ballistikers gehört, so wird sich das Werk 
neben den übrigen Darstellungen zweifellos einen 
guten Platz erobern. 

Essen. K. Stange. 410 

Dr.-Ing. RICHARD EMIL KUTTERER, Balli- 
stik. mit einem Vorwort von Geh. Reg.-Rat Dr. 
C,Cranz. (Die Wissenschaft. Einzeldarstellungen 
aus der Naturwissenschaft u. der Technik. Bd. 97, 
herausgegeben von Prof. Dr. Wilhelm Westphal.) 
VII + 208 S. Braunschweig 1942, Verlag Friedr. 
Viewer & Sohn. Preis geb. 13,50 M. 

Als 97. Band in der Reihe „Die Wissenschaft. 
Sammlung von Einzeldarstellungen aus den Ge- 
bieten der Naturwissenschaft und der Technik“ er- 
schien diese Einführung in die Ballistik von Re- 
gierungsrat Dr.-Ing. Richard Emil Kutterer, 
einem Schüler von ©. Cranz, der diesem Werk 
ein anerkennendes Vorwort schrieb. 

Das Buch betont vorwiegend die technisch- 
physikalische Seite der Ballistik. Mathematische 
Ableitungen werden durchweg möglichst elementar 
gehalten, z. T. auch ohne Beweis angegeben, je- 
doch meist so. daß der Anschluß an eingehendere 
theoretische Darstellungen möglich bleibt. Eine 
klare Darstellungsweise, ausgezeichnete Aufnahmen 
des fliegenden Geschosses und eine große Anzahl 
eindringlicher Zahlenbeispiele machen dieses Werk 
zu einer gut verständlichen Darstellung der Grund- 
züge der Ballistik. 

Aus dem ersten — verhältnismäßig kurzen — 
Abschnitt „Äußere Ballistik“ sind zu erwähnen die 
ausführliche Darstellung der experimentellen Me- 
thoden und die Einbeziehung aerodvnamischer 
Gesichtspunkte in einzelne Probleme der äuße- 
ren Ballistik. Die mathematischen Grundlagen 
zur Berechnung von Flugbahnen werden nur 
kurz gestreift. Der zweite Teil „Innere Ballistik“ 
ist ziemlich eingehend: auch er bringt neben der 
theoretischen Behandlung des innerballistischen 
Gleichungssystems hauptsächlich eine Darstellung 
experimenteller Methoden der Innenballistik, wo- 
bei der Verfasser umfangreiche eigene Erfahrungen 
auf diesem Teilgebiet mit Nutzen verwertet. Den 


XUM 


300 


Buchbesprechungen 


Z. angew. Math. Mech. 
Bd. 22 Nr.5 Okt. 1942 


Schluß des Buches bilden Tabellen- und Formel- 
auszüge sowie Schrifttums-, Namen- und Sach- 


verzeichnisse. 

Das Buch ist dem Neuling in der Ballistik eben- 
so zu empfehlen wie den erfahrenen Ballistikern, 
unter denen gerade die reinen Theoretiker dem 
Verfasser für diese klare Einführung in den tech- 
nisch-physikalischen Teil der Ballistik dankbar 
sein werden. 

Berlin. H. Athen. 392 


THEODOR VAHLEN, Ballistik. 2. Aufl. neu- 
bearb. u. herausgeg. unter Mitwirkung von Alfred 
Klose. X + 267 S. m. 65 Abb. Berlin 1942, Ver- 
lag Walter de Gruyter & Co. Preis geb. 18 M. 

Das Buch bringt in 18 Kapiteln eine umfassende 
Übersicht über den gesamten ballistischen Problem- 
kreis. Behiandelt werden u. a.: die in der Außen- 
ballistik wirkenden Kräfte, Grundlegung der Be- 


wegungsgleichungen und allgemeine Flugbahn- 
eigenschaften, Reihenentwicklungen und speziell 


Potenzreihen, Grenzfälle und allgemeine Lösungs- 


klassen, praktische Lösungsmethoden, Geschoß- 
bahnstörungen, konische Pendelung, Drallfragen, 


Übergangsballistik, Innenballistik, ballistische 
Wahrscheinlichkeitsrechnung, Endballistik. 

Die 2. Auflage unterscheidet- sich in ihrem Ge- 
samtcharakter nicht von der 1. Auflage (1922). In 
einzelnen Kapiteln sind jedoch wesentliche Neue- 
rungen aufgenommen worden, davon besonders er- 
wähnenswert: die veränderte Behandlung des Luft- 
widerstandes, die Grundlegung der außenballisti- 
schen Bewegungsgleichungen, das Kapitel über die 
praktischen Lösungsmethoden der Außenballistik, 
die Betrachtungen über Iterationslösungen in der 
„kosmischen“ Ballistik und das neue innenballisti- 
sche Lösungsverfahren von Vahlen. 

Das Werk- vermittelt durch schärfste Syste- 
matik und mathematische Sauberkeit einen hervor- 
rarenden Einblick in den wissenschaftlichen Grund- 
gedanken der theoretischen Ballistik. Insofern ist 
im Rahmen des vorliegenden Buches das von 
Vahlen erstrebte Ziel eines klaren und systemati- 
schen Aufbaues der Ballistik voll erreicht worden. 

Berlin. Dr. H. Athen. 411 


Dr.-Ing. HANS KNOBLOCH, Die BAllistik 
in der Luftwaffe, vom Bombenwerfen, 
Schießen und vom Treffen. (Luftfahrt- 
Lehrbücherei, Bd. 12.) 78 S. m. 72 Abb. Berlin 1942, 
Verlag Dr. M. Matthiesen & Co. Preis kart. 3,80 M. 

Das Buch gibt einen leichtverständlichen Ein- 
blick in die Probleme der Ballistik bei der Luft- 
waffe. Neben den allgemeinen Begriffen, wie 
Schwerkraft, Luftwiderstand, Drall und Pfeilflug, 
wird die Bombenballistik vom Standpunkt des 
Soldaten, des Meßtechnikers und des Ballistikers 
mit der Darstellung des Bombenwurfes aus dem 
Horizontaltlug, Sturzflug und Tiefflug, sowie des 
Wurfes von Brandbomben ausführlich behandelt. 
Ein kurzer Abschnitt hat die Ballistik für das 
Schießen vom Flugzeug aus zum Gegenstand; nach 
der Schilderung allgemeiner Voraussetzungen wird 
sie vom Standpunkt des ruhenden und des mit- 
bewegten Beobachters behandelt. 

Das Buch wendet sich an alle, die an den spezi- 
ellen Fragen der Ballistik der Flugzeugwaffen In- 
teresse haben: für den Leser mit mathematischen 
Elementarkenntnissen sind viele Beziehungen zu- 
sätzlich noch in Formeln gefaßt. Die Darstellungs- 
weise des Büchleins ist klar und einfach: durch 
zahlreiche Abbildungen wird das Verständnis ver- 
tieft. Dem Charakter des Buches entsprechend 


wäre ein Eingehen auf die Wirkung der Bomben 
und Geschosse im Ziel, sowie auf die besonderen 
Fragen der Treffwahrscheinlichkeit beim Schießen 
im Luftkampf erwünscht gewesen. 


Dr. H. Athen. 412 
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Dr. FRANZ BAUR, Prof. a. d. Universität Frank- 
furt/M., Einführung in die Großwetter- 
forschung. (Math.-Phys. Bibl. Reihe I: heraus- 
zegeben von W. Lietzmann und W. Kliefoth, 
Bd. 88.) 2. verb. Aufl. 53 S. m. 14 Abb. im Text. 
Leipzig-Berlin 1942, Verlag B. G. Teubner. Preis 
kart. 1.20 M. 


Die kleine, aber außerordentlich inhaltreiche 
Schrift von F. Baur „Einführung in die Groß- 
wetterforschung“ ist in zweiter, verbesserter Auf- 
lage erschienen. In ihr behandelt der Verfasser, 
der als Forscher auf dem Gebiet der Langfrist- 
vorhersage und als Leiter des Forschungsinstitutes 
für langfristige Vorhersagen des Reichswetter- 
dienstes allgemein bekannt geworden ist, nach- 
stehende Kapitel: „Großwetter und mittlere Luft- 
druckverteilung — Wärmehaushalt und allgemeiner 
Kreislauf der Atmosphäre — Witterungsperioden 
und Wetterrhythmen — Geophysikalische Eintlüsse 
auf das Großwetter — kosmische Einflüsse auf das 
Großwetter — Die langfristige Wettervorhersage.“ 
Gegenüber der ersten Auflage sind vor allem die 
Abschnitte über den Wärmehaushalt. über Luft- 
druckrhythmen und über die 10tägige Wetter- 
vorhersage verbessert und erweitert worden. Die 
Darstellung ist so anschaulich und fesselnd und 


die Bedeutung der behandelten Probleme — vor 
allem für wirtschaftliche Fragen — so groß, daß 


dieser Schrift eine recht 
wünschen ist. 


z. Zt. Großenhain. 


weite Verbreitung zu 
H. Goldschmidt. 407 
WERNER HEISENBERG, Wandlungen in 


den Grundlagen der Naturwissen- 
schaft. 6 Vorträge. 3. erw. Aufl. 95 S. Leip- 


zie 1942, Verlag S. Hirzel. Preis kart. 3,50 M. 


Zu den Aufsätzen der beiden ersten Auflagen 
(„Wandlungen der Grundlagen der exakten Natur- 
wissenschaft in jüngster Zeit“ aus dem Jahre 34; 
„Zur Geschichte der physikalischen Naturerklä- 
rung“, 32: „Prinzipielle Fragen der modernen Phy- 
sik“, 35) sind drei weitere hinzugekommen: „Ge- 
danken der antiken Naturphilosophie in der mo- 
dernen Physik“ (Antike, Band XID: „Die Goethe- 
sche und die Newtonsche Farbenlehre im Lichte 
der modernen Physik“ (41) und „Die Einheit des 
naturwissenschaftlichen Weltbildes“ (41). Der hohe 
Genuß, mit dem man jeden einzelnen dieser Auf- 
sätze liest. beruht einmal darauf, daß hier ein For- 
scher zu uns spricht. der entscheidende Wendungen 
in der Entwicklung der Physik unserer Tage selbst 
mit herbeigeführt hat, dessen Führung man sich 
also wie kaum einem anderen anvertranen kann. 
Dazu kommt aber, daß die Gedankengänge in einer 
Klarheit und Durchsichtigkeit ablaufen, daß die 
Darstellung durch eine Schönheit und Treffsicher- 
heit des Ausdrucks gekennzeichnet ist, die man 
leider nicht allzu häufig in den für ein weiteres 
Publikum bestimmten Aufsätzen wissenschaftlichen 
Inhalts findet. Unwillkürlich wird man an (die 
glänzenden Helmholtzschen populärwissenschaft- 
lichen Vorträge erinnert. Interessant ist dabei, 
die Wandlung in der Grundeinstellung zur Wissen- 
schaft festzustellen, die in den Heisenbergschen 
Aufsätzen zum Ausdruck kommt: trotz zuversicht- 
licher Stimmung, trotz des Vertrauens in den Fort- 
schritt der Wissenschaft finder wir doch darin eine 
zurückhaltende Bescheidenheit, eine, man könnte 
vielleicht sagen agnostische Haltung, die der 
Wissenschaft des ausgehenden 19. Jahrhunderts 
fast unbekannt war. Heisenberg charakterisiert 
diese Einstellung in seinem letzten Aufsatz mit fol- 
genden Worten: „Wir sind uns mehr als die frühere 
Naturwissenschaft dessen bewußt, daß es keinen 
sicheren Ausgangspunkt gibt, von dem aus Wege 
in alle Gebiete der Erkenntnis führen, sondern daß 
alle Erkenntnis gewissermaßen über einer grund- 
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losen Tiefe schweben muß, daß wir stets irgendwo 
in der Mitte anfangen müssen, über die Wirklich- 
keit zu sprechen mit Begriffen, die erst durch ihre 
Anwendung allmählich einen schärferen Sinn er- 
halten. und daß selbst die schärfsten, allen Anfor- 
derungen an logischer und mathematischer Präzi- 
sion genügenden Begriffssysteme nur tastende Ver- 
suche sind, uns in begrenzten Bereichen der Wirk- 
lichkeit zurechtzufinden.“ 


Berlin. E. Mosch. 408 


Dr. phil. nat. habil. MAX STECK, Dozent der 
Math. a. d. Techn. Hochschule München, Das 
Hauptproblem der Mathematik. XVI 
+232 S. Berlin 1942, Dr. Georg Lüttke Verlag. 
Preis geb. 10.60 M. 

Verf. sieht das Hauptproblem der Mathematik in 
der Untersuchung ihrer denkerischen Grundlagen, 
aus denen durch logische Operationen die Erkennt- 
nisse gewonnen werden, „die von apodiktischer 
Sicherheit, die einfach die Wahrheit selbst sind“. 
Dieses Hauptproblem stellt sich ihm unter zwei 
Aspekten dar: Der Problemkreis 4 befaßt sich mit 
der Aufstellung eines Axiomsystems im Sinne des 
Formalismus; der Problemkreis B hat die Gewin- 
nung gewisser ideeller Gebilde zum Ziel, hinter 
denen denkerisch reine Ideen des Mathematischen 
stehen, die ebenfalls axiomatisch erfaßt werden 
müssen, und die dann erst die wahren Grundsätze 
der Mathematik liefern. Die mathematische For- 
schung unserer Zeit befindet sich nach Meinung des 
Verf. auf falschem Wege. Sie hat sich dem For- 
malismus und Logizismus verschrieben. Die Denk- 
weise der von ihm so genannten „deutschen Linie“, 
die in Plato ihren Ursprung hat und über die grie- 
chischen Mathematiker und den Neuplatonismus zur 
mittelalterlichen deutschen Mystik und weiter über 
Nicolaus Cusanus, Dürer, Copernieus und Kepler 
zu Leibniz und Kant führt und noch in den nach- 
kantischen Idealismus hineinreicht, mußte infolge 
des Einflusses namentlich jüdischer Forscher der 
„westlichen Linie“ weichen, die bei den Sophisten 
beginnt und über Aristoteles schließlich zu Des- 
cartes (auch Galilei ist davon nicht unberührt ge- 
blieben) und dann über die Aufklärung zum Posi- 
tivismus des 19. Jahrhunderts führt. Die Mathe- 
matiker der deutschen Linie, denen eine „urbildlich- 
noötische“ Denkweise eigen ist, waren sich noch 
der unzertrennlichen Verbundenheit der beiden 
Problemkreise A und B bewußt, während die west- 
liche Linie mit ihrer „abstrakt-nominalistischen“ 
Denkweise dem Problemkreis B verständnislos ge- 
genübersteht und nicht die reine Wahrheit 
anstrebt, sondern diesen Begriff durch den der 
Widerspruchslosigkeit des Axiomen- 
systems 4 ersetzt. Mathematische Hauptrichtun- 
zen dieser westlichen Linie sind der Empirismus 
(Pasch), der Formalismus (Hilbert) und die Logistik 
(Russel - Whitehead - Scholz). Verf. bemüht sich, 
namentlich auch mit Hilfe zahlreicher Äußerungen 
mathematischer Forscher aus allen Lagern — die 
Belegstellen nehmen allein etwa 40% des Buches 
ein —, die Unrichtigkeit und Schädlichkeit jener 
Richtungen unter Beweis zu stellen, wobei er be- 
sondern heftig den Hilbertschen Standpunkt an- 
greift. Der Ton seiner Ausführungen ist oft reich- 
lich scharf: seine Erregung macht sich in 
zwischen die gegnerischen Darlegungen einge- 
streuten ironischen Bemerkungen und Ausrufungs- 
zeichen Luft, die, ebenso wie etwa der persönliche 
Aneriff gegen Scholz (S. 149, Anm.), kein sach- 
liches Material zur Widerlegung der Gegner bei- 
steuern. — Verf. sieht in der Platonischen Ideen- 
lehre und dem Kantschen erkenntnistheoretischen 
Standpunkt die Ausgangspunkte zur Lösung des 
Hauptproblems, einer Lösung, die er in einer Syn- 
these von Formalismus und Intuitionismus erblickt. 


‚werden (dürfte. 


— Da es sich beim mathematischen Grundlagen- 
streit doch wohl um gegensätzliche geistige Grund- 
einstellungen erkenntnistheoretischer, ja metaphv- 
sischer Art handelt, so dürften auch die tempera- 
mentvollen Ausführungen des Verf. kaum imstande 
sein, Andersgläubige zu bekehren. die dann aller- 
dings das herbe Urteil des Verf. auf sich nehmen 
müssen (S. 213), daß sie das Wesen des Mathe- 
matischen überhaupt nieht erkannt haben, ein Ur- 
teil. durch das die Bedeutung eines Mathematikers 
vom Range Hilberts — der doch wohl auch jenem 
Vorwurf ausgesetzt wäre — kaum geschmälert 
Verf. schließt seine Arbeit mit der 
Mahnung, auch in der Mathematik das Erbe, das 
die großen Deutschen und die arischen Forscher 
denkerisch errungen haben, getreulich zu verwalten 
und zu mehren. 

Berlin. . Mosch. 409 

Prof. Dr. PAUL RIEBESELL, Mathematik 
des täerlichen Lebens. 52 S. Hamburg 
1942, Hanseatische Verlagsanstalt. Preis kart. 2M. 

Plaudereien über Fragen des täglichen Lebens, 
die mit Hilfe des Zahlenrechnens, der Mathematik 
oder der Statistik behandelt werden können. Das 
Heft wird hoffentlich viele anregen, über die ange- 
schnittenen Probleme nachzudenken und sich ein- 
gehender mit ihnen zu beschäftigen. Allerdings 
werden zum Teil Kenntnisse und vor allem eine 
Denkfähigkeit vorausgesetzt, wie wir sie heute z.B. 
nur bei sehr wenigen der Abiturienten unserer 
höheren Schulen finden. 


Dresden. Willers. 391 

Prof. Dr.-Ing. P. WERKMEISTER, em. 
Prof. a. d. Techn. Hochschule Dresden, Ver- 
messungskunde IL Stückmessung und 
Nivellieren. 7. Aufl. (Samml. Göschen. Bd. 
468). 168 S. mit 145 Fig. Berlin und Leipzig 
1942, Verlag Walter de Gruyter u. Co. Preis geb. 
1,62 M. 

Nach der im Jahre 1938 vom Verlag de Gruyter 
herausgegebenen sechsten Auflage des Bandes 468 
der Sammlung Göschen, der aus dem Gebiete der 
Vermessungskunde die Kapitel „Stückmessung 
und Nivellieren“ behandelt, erscheint im gleichen 
Verlag schon: wieder eine neue, die siebente Auf- 
lage des genannten Bändchens, ein Beweis für die 
berechtigte Beliebtheit auch dieser Arbeit von 
P. Werkmeister. Inhaltlich hat sich an den 
vorzüglichen Ausführungen nichts geändert, und 
auch in der Ausstattung hat die Neuauflage trotz 
der kriegsbedingten Schwierigkeiten in keiner 
Weise gelitten. 

Freiberg-Sa. 


ord. 


H. Müller. 394 


Ferner sind bei der Schriftleitung folgende 


Bücher eingegangen 
bleibt vorbehalten): 


Dr.-Ing. FRIEDRICH KLINGER, Einführung 
in die n-dimensionale algebraische 
Geometrie, mit besonderer Berück- 
sichtigunge der Physik für Teeh- 
niker und Physiker. Herausgeg. mit Unter- 
stützung der Akademie der Wissenschaften in 
Wien. VIE + 152 S. m. 1 Tafel. Wien 1942, Sprin- 
ger-Verlag. Preis brosch. 12 M. 


Dr. KARL GEY, Oberstudiendirektor in Leipzig. 
Vierstellige Tafel der Logarithmen 
und Antilogarithmen, Praktisches 
Hilfsmittel für Überschlags-. Kon- 
troll- und abgekürzte Rechnungen, 
Anhang: Die trigonometrischen Funktionen für die 
ganzen Grade. 8. S. Bonn und Berlin 1942, Ferd. 
Dümmlers Verlag. Preis brosch. 0,90 M. 
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Rudolf Rothe zum Gedächtnis. 

Wenige Tage nach Vollendung seines 69. Lebens- 
jahres starb am 26. Oktober 1942 in Berlin der ord. 
Professor Dr. Rudolf Rothe, Emeritus und Ehren- 
senator der Technischen Hochschule Berlin. Mit 
ihm verliert die Deutsche Mathematikerschaft einen 
namhaften Gelehrten und die Berliner Technische 
Hochschule einen ihrer angesehensten Lehrer. (ler 
weit über die Grenzen seines Faches hinaus be- 
kannt und beliebt war. Wir 
Glück hatten, Rudolf Rothe näher kennenzulernen, 
die wir als seine Schüler von ihm lernten, als 
Assistenten ihm zur Seite standen, als Kollegen 
mit ihm gemeinsam arbeiteten. wir sehen in Rudolf 
Rothe mehr als nur einen bedeutenden Wissen- 
schaftler und hervorragenden Pädagogen, wir 
trauern auch um den Menschen: (die Lauterkeit 
seines Wesens, (len Ernst seiner Lebensauffassung, 
seine gelassene Heiterkeit, die auch den anstren- 
eendsten Arbeitstag überdauerte, die vorbildliche 
Ptlichttreue, die Nachsicht. mit der er die Schwä- 
chen der Mitmenschen ertrug, die Strenge gegen 
sich selbst — dies alles liebten wir an Rudolf 
Rothe, dem feinsinnigen musischen Menschen, dem 
gütigen Mentor, dem aufrechten deutschen Mann, 
der am 27. Januar 1919 seine Vorlesungen aus- 
fallen ließ, weil der Kaiser Geburtstag hatte, dem 
ganzen Kerl. der noch mit fast 60 Jahren das 
Matterhorn bestiee! Rudolf Rothe wird uns allen 
immer unvergeßlich bleiben. 

Rudolf Ernst Rothe wurde am 15. Oktober 1873 
in Berlin geboren. Durch den frühen Tod des Vaters. 
der seine Familie in engen wirtschaftlichen Ver- 
hältnissen zurückließ,. lernten Rudolf und seine 
drei jüngeren Schwestern schon frühzeitig (den 
Ernst des Lebens kennen. Rudolf besuchte das 
Sophiengymnasium,. das er 1892 mit dem Reife- 
zeugnis verließ. Er studierte dann an der Univer- 
sität Berlin, wo damals H. A. Schwarz Weierstraß’ 
geistiges Erbe betreute und mehrte, Mathematik 
und Naturwissenschaften. Neben H. A. Schwarz 
und L. Fuchs gewann vor allen der Flä-hentheore- 
tiker Johannes Knoblauch Einfluß auf den jungen 
Rothe, mit dem er persönlich eng befresndet war. 
Die frohe Kameradschaft im Mathematischen Ver- 
ein an der Berliner Universität, dem Rothe auch 
noch später als Alter Herr die Treue hielt, sorgte 
dafür, daß die strenge Wissenschaft dem jungen 
Studenten nicht über den Kopf wuchs. Übrigens 
zeigte Rothe schon damals Verständnis für 
die Frage (der mathematischen Ingenieurausbil- 
(dung, indem er in den Übungen von Lampe an der 
Teehnischen Hochschule Charlottenburg assistierte, 
sehr zum Mißverenügen von H. A. Schwarz. der 
es gar nicht verstehen konnte. was ein reiner Ma- 
thematiker in Charlottenburg zu suchen hatte, „wo 
die Leute mit dem Reißbrett über die Straße 
eehen“. Rothe promovierte 1897 mit einer differen- 
tialgeometrischen Arbeit (Untersuchungen über (die 
allgemeine Theorie der isothermen Flächen), an 
die sich eine ganze Reihe schöner Beiträge zur 
Flächentheorie anschloß. Nach Beendigung des 
Studiums trat Rothe 1897 als wissenschaftlicher 
Hilfsarbeiter bei der Physikalisch - Technischen 
Reichsanstalt ein, in der damals Warburg regierte. 
In den folgenden Jahren veröffentlichte Rothe 
neben seinen differential-geometrischen Arbeiten 
verschiedene physikalische Untersuchungen, die 
erkennen lassen, mit welchem Eifer er sich 
seinen Aufgaben bei der Reichsanstalt wimete, 
bei der’ er inzwischen etatsmäßiger Mitarbeiter 
geworden war (1901). Noch nach langen 
Jahren dachte Rothe gern an diese Zeit zurück, 
wobei er stets versicherte, daß er dort bei der 
Physikalisch-Technischen Reichsanstalt erst riehtig 


NACHRICHTEN 


aber. die wir das, 


Physik gelernt hätte. 1905 habilitierte sich Rothe 
an der Technischen Hochschule Charlottenbure und 
wurde drei Jahre später als Nachfolger von Horn 
nach Ulausthal berufen, wo er neben der Hö- 
heren Mathematik auch die Mechanik vortrug. 1914 
folgte Rothe einem Ruf nach Hannover und kehrte 
schon im folgenden Jahr nunmehr endgültig an 
die Technische Hochschule Charlottenburg auf 
den Lehrstuhl von Hettner zurück. Neben 
seiner Lehrtätiekeit an der Technischen Hoch- 
schule hielt Rothe im Weltkriex zeitweise auch an 
(der Berliner Universität Vorlesungen über Diffe- 
rential- und Integralrechnune und Differential- 
geometrie. Vor allen Dingen aber führte Rothe da- 
mals gemeinsam mit Geheimrat Uranz, dem Nestor 
(ler deutschen Ballistik, wichtige Untersuchungen 
durch, die mit Rücksicht auf die erforderliche Ge- 
heimhaltung gar nicht oder erst sehr viel später 
veröffentlicht werden «durften. So z. B. bestimmte 
Rothe auf zeichnerischem Wege Flugbahn und 
Geschwindigkeit der Granaten des Parisgeschützes. 
Auf den Zusammenbruch vom Jahre 1918, unter 
dem Rothe. ein Mann von reinster vaterländischer 
Gesinnung, unsagbar gelitten hat, folgten lange 
Jahre rastloser, wissenschaftlicher. pädagogischer 
und auch administrativer Tätigkeit, war doch 
Rothe in der schwierigen Intlationszeit Rektor der 
Technischen Hochschule und später noch mehrere 
Jahre Dekan der allgemein-wissenschaftlichen Fa- 
kultät. Bei allen diesen Pflichten fand Rothe noch 
Zeit und Kraft zu vielseitiger literarischer Arbeit: 
in den ‚Jahren 1925 bis 35 entstanden die drei 
Bände seiner „Höheren Mathematik“. 1927 er- 
schienen die von Rothe bearbeiteten Vorlesungen 
über Variationsreehnune von Karl Weierstraß 
(1915 hatte Rothe schon den Band „Anwendungen 
der elliptischen Funktionen“ herausgegeben). 1931 
schrieb er (gemeinsam mit Öllendorf und Pohl- 


hausen) die angewandte Funktionentheorie, 1937 
die  Differentialgeometrie. Daneben bearbeitete 


Rothe in mehreren Auflagen des Ingenieurtaschen- 
buches „Hütte“ den mathematischen Abschnitt. 
hielt im Telegraphentechnischen Reichsamt und im 
Außeninstitut der Technischen Hochsehule Vor- 
träge für in der Praxis stehende Ingenieure und 
leitete ein schulmathematisches Kollogainm, in dem 


er mit den Lehrern der Höheren Schule daran 
arbeitete, die Beziehungen zwischen Schule und 
Hochschule enger und fruchtbarer zu gestalten. 


Schließlich widmete Rothe auch noch der studen- 
tischen Sozialfürsorge als Mitglied der verschiede- 
nen Ausschüsse Zeit und Kraft. Als nach Wieder- 
herstellung der Wehrhoheit in Deutschland an der 
Techmischen Hochschule eine wehrtechnische Fa- 
kultät eingerichtet wurde, übernahm Rothe zu den 
alten Pflichten noch das Ordinariat für Ballistik. 
nachdem er schon vorher regelmäßig über dieses 
Thema Vorlesungen gehalten hatte. Seine eigenen 
wissenschaftlichen Veröffentlichungen im einzelnen 


aufzuführen ist hier nicht der Ort. Zu den phy- 
sikalischen Arbeiten der Frühzeit und zu den 
zahlreichen differential - geometrischen Veröffent- 


lichungen gesellten sich im Laufe der Jahre noch 
viele verschiedensten Inhalts: über Fragen der 
praktischen Analysis, über den Mittelwertsatz unıd 
(die Taylorsche Formel, über wiederholte Quadra- 
turen. über die Abelsche Integralgleichung, über 
schulmathematische und  ballistische Probleme. 
Überragend war Rothes pädagogische Kunst: wie 
vorzüglich verstand er es, seine Zuhörer zu fesseln! 
Das Auditorium Maximum der Technischen Hoch- 


schule war in den ‚Jahren des großen Andrangs zu 
(den Hochschulen bis auf den letzten der 650 Plätze 
und darüber hinaus gefüllt: auch auf den Stufen 
und Fensterbrettern saßen die Studenten und hör- 
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ten ihrem verehrten Professor Rothe andächtig zu. 
In seiner frischen humorvollen Art, mit kraftvoller 
Stimme und vorbildlich klarer Schrift an der Tafel 
machte Rothe seinen Hörern die ihnen von Hause 
aus kaum sympathische Mathematik so schmack- 
haft, daß sie geduldig aushielten, bis Rothe es für 
gut befand, die Vorlesung zu schließen. Die Stu- 
(lenten wußten aber auch, daß sie an Rothe den 
verständnisvollsten und gütigsten Freund hatten, 
bei dem sie jederzeit für ihre Nöte und Anliegen 
ein offenes Ohr fanden, und bewiesen ihm ihre Zu- 
neigung bei jeder Gelegenheit. Auch zwischen 
den Bergleuten an der Technischen Hochschule 
und Rothe bestand eine sicherlich nicht wissen- 
schaftlich — Rothe prüfte die Bergleute im Vor- 
examen —, sondern nur menschlich und  histo- 
risch (von seiner Clausthaler Zeit her) fundierte 
Freundschaft, die sich bei dem alljährlich wieder- 
kehrenden Barbarafeste aufs herzlichste äußerte. 
Rothes besondere Liebe galt natürlich den jungen 
Mathematikern und Physikern, die, soweit sie be- 
dürftig und würdig waren, aus den Erträgnissen 
der Rotheschen „‚JJohannastiftung” ein Stipendium 
erhalten konnten. 

Rothes- Arbeitskraft schien unverwüstlich. We- 
niee Urlaubswochen, (die er in seinen geliebten 
bayerischen oder tiroler Bergen verlebte, genügten 
völlig, um ihn für ein weiteres Jahr angestrengte- 
ster Arbeit zu stärken. Dabei waren diese Ferien 
keineswegs der Ruhe gewidmet. sondern mit be- 
achtlichen alpinistischen Leistungen ausgefüllt, von 
den nie abreißenden Korrekturarbeiten ganz zu 


schweigen, die Rothe auf jeder Reise begleiteten. 


ZUSCHRIFTEN AN 


H. Peineke, Fehlerbetrachtung zur Rinaldini- 
schen Regel. Z. angew. Math. Mech. Bd. 22 (1942), 
Ss. 116—117. 


1. Fehlerbetrachtungen zur Rinaldinischen Regel 
und auch zur Näherungskonstruktion des Herzogs 
Karl Bernhard zu Sachsen-Weimar-Eisenach für die 
Seite des regelmäßigen n-Eckes sind schon 1878 
von S. Günther angestellt worden (Zeitschrift 
für Realschulwesen, Wien, Bd. 3 (1878), S. 526 
bis 533 und S. 704): ihre Ergebnisse sind aus- 
zugsweise in dem bekannten Werk Konstruk- 
tionen und Approximationen von Th. Vahlen 
(Leipzig u. Berlin 1911) auf S. 302ff. wieder- 
eegeben. Die Güntherschen Fehlerkurven wurden 
nachgeprüft und teilweise berichtigt von H. J. 
Fischer und Kurt Schmeiser (Unter- 
suchungen zur angenäherten Kreisteilung, Sitzungs- 
berichte der Heidelberger Akad. d. Wiss., math.- 
naturw. Klasse, 1934, 18. Abhandlung). 

In der zuletzt genannten Arbeit von H. J. 
Fischer und K. Schmeiser wird auch eine 
im Jahre 1931 von dem Lehrer Jakob Klee an- 
gegebene sehr einfache Näherungskonstruktion zur 
Teilung des Viertelkreises in beliebig vorgegebe- 
nem Verhältnis auf ihre Genauigkeit untersucht. 
Das Fehlermaximum liegt bei ihr zwischen 86” und 
87”, d. h. der Abstand zwischen wahrem und kon- 
struiertem Teilpunkt auf der Kreislinie kann maxi- 
mal erst dann 1 mm groß werden, wenn der Durch- 
messer des geteilten Kreises 4.80 m lang ist! Be- 
nutzt man die Kleesche Konstruktion insbeson- 
dere zur Bestimmung des Mittelpunktwinkels des 
regelmäßigen n-Eckes, so liefert sie in mancherlei 
Hinsicht viel bessere Ergebnisse als die Methoden 
von Rinaldini und Karl Bernhard von Sach- 
sen-Weimar-Eisenach. Beide weisen ein weit größe- 
res Fehlermaximum von mehreren Minuten auf. Das 
Produkt ng, aus Seitenzahl n und Kkonstruiertem 
Mittelpunktwinkel 9, strebt, falls n gegen Unendlich 


Auch als Rothe im Jahre 1939 emeritiert wurde, 
bedeutete «das bei ihm natürlich nieht. daß er nun 
die Hände in den Schoß legte. In den letzten 
Jahren wurde der Unermüdliche allerdings durch 
ein seit 1938 immer wieder aufflackerndes Leiden 
zeitweise zur Untätigkeit gezwungen. Und nun hat 
der Tod diesem an Arbeit und Erfolg überreichen 
Leben ein Ziel gesetzt, zu früh für alle, die Ru- 
dolf Rothe kannten und ihn nie vergessen werden: 
zu früh, obwohl nach menschlichem Ermessen nie- 
mand mehr von seinem Leben erwarten darf. als 
Rothe so reichlich zuteil wurde. War doch auch 
sein Leben erfüllt im Sinne der trostreichen Worte 
von Wildenbruch: 

Sterben ist nur eines Tages Enden, 

Tod nur Schlaf der niemals wach Gewesenen 

Nie entschläft. wer einmal wach zelebt. 


Freiberg. Grüß. 431 


Persönliches. 


Ende September verstarb der Obmann des Aus- 
schusses für Rechentechnik im Reichsausschuß 
für wissenschaftliche Fertigung Stwudienrat H. 
Schwerdt. der durch seine Bücher über No- 
mographie bekannt geworden ist. 

Dem o. Prof. Geheimen Rat Dr. Dr. e. h. Se- 
bastian Finsterwalder, München, wurde aus 
Anlaß der Vollendung des achtzigsten Lebensjahres 
in Anerkennung seiner Verdienste um die mathe- 
matischen Wissenschaften die  Goethemedaille 
verliehen. 


DEN HERAUSGEBER 


geht, bei Rinaldini gegen 396 ° 57’24”, bei Karl 
Bernhard gegen 362° 22’ 12”, bei Klee gegen 
360 ° 32° 30°, d. h. beinahe gegen den Sollwert 
360 ° 00° 00”. Der relative Fehler ist für große n 
bei der Kleeschen Konstruktion also besonders 
klein. Sie verdient deshalb durchaus, neben den 
beiden älteren Konstruktionen genannt zu werden. 
Da die Kleesche Konstruktion nicht allgemein 
bekannt sein dürfte, sei sie hier kurz beschrieben: 
BAC sei der gegebene rechte Winkel, 


AD:DB=s:t 


das vorgeschriebene Streckenverhältnis. In B uni 
D werden Lote auf AB errichtet: DE wird gleich 


DB, der Winkel BDG gleich 8 gemacht. Dann 


“Bild 1. 


liefert der Kreis um G mit dem Radins EG den 
gesuchten Näherungsteilpunkt H, d. h. der Punkt H 
teilt den Viertelkreis BC mit der obengenannten 
Genauigkeit im Verhältnis s:t. 

Tübingen. Max Müller 
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2, Es liegt nahe, durch eine kleine Änderung der 
Rinaldinischen Regel dafür zu sorgen, dab 
die stark verschiedenen positiven und negativen 
Höchstwerte von- Ja gleich werden. Dies ist zu 
erreichen, wenn man die Abstände MO=MP 
statt = rtg60° nun—rtgy macht, wo y etwas 
kleiner als 60° ist. Der Fehler da berechnet sich 
dann nach der Formel 


tg 
tg’y+ 2° 


cos (a + 


Für 7—=59° 122° schwanken die Fehler Aa 
zwischen + 19%’ und 19%’. Legt man Wert 
darauf, daß die Konstruktion mit Zirkel und 
Lineal allein ausführbar bleibt, so kann man 
»=59° 89 wählen, wodurch die Entfernungen 


der Punkte und P von r werden, 


L 
20 
Bild 2. 


Die Fehler 4a liegen dann zwischen + 18%’ und 
— 20%’. Das Bild 2 zeigt den Verlauf der Fehler- 
kurven für Ja als Funktion von « für die beiden 
Fälle »=60° und 7=59° 122°. 

Sind also Fehler des Zentrierwinkels « bis zu 
Höchstwerten von +20’ zulässig, was bei kleinen 
n-Werten öfter der Fall sein wird, so könnte die 
in der angegebenen Weise veränderte Rinaldini- 
sche Regel ohne weiteres benutzt werden. 


Berlin. H. Manrer. 403 


Die Zuschriften haben Herrn Peineke vor- 


gelegen. 


Berichtigungen. 


D. Küchemann, Bemerkung über den Ein- 
fluß des Seitenverhältnisses auf den Zusammen- 
hang von Normalkraft und Tangentialkraft eines 
be 2 Z. angew. Math. Mech. Bd. 20 (1940) 
290 bis 


Herr Prof. Tollmien hat mich freundlicher- 
weise auf ein Versehen in der obigen Mitteilung 
aufmerksam gemacht, das ich hiermit berichtigen 
möchte. 


Ich habe bei der Ableitung stillschweigend vor- 
ansgesetzt, daß die Flügelsehne in der Nullauftriebs- 
richtung des Profils liegt und Tangential- und Nor- 
malkraft parallel bzw. senkrecht zu dieser Rich- 
tung definiert. Da jedoch in der Praxis der Profil- 
messungen die Sehne meist nach geometrischen Ge- 
sichtspunkten gewählt wird, trifft die gemachte 
Voraussetzung im allgemeinen nur bei symme- 
trischen Profilen zu. Es ist daher zweckmäßig, den 


Winkel a, zwischen geometrischer Sehne und Null- 
auftriebsrichtung einzuführen. In Gl. (1) ist statt 
a dann @«—a, zu setzen, während in den Gl. (2) 
und (3) auf der rechten Seite ein Glied — €„-@ 
hinzukommt. Die Gültigkeit des aufgestellten 
Satzes sowie die allgemeinen Überlegungen dazu 
werden also durch den neuen Zusatz nic ht berührt. 


| 
- 4 - 12- 
Jeitenverköltnis 
N 3 
:6 
:7 
En? — 
0-0.84. 00.389 
| 


Bild 1. 


— Bei dem in Bild 1 verwandten Profil Gö 389 be- 
trägt a, etwa —5°. Die Differenz zwischen den 
früher eingetragenen theoretischen und den ge- 
messenen Werten ist also nicht der Profilwider- 
stand. In-die dieser Berichtigung beigefügten Ab- 
bildung ist die vollständige theoretische Beziehung 
eingetragen, wobei die Konstanz des Profilwider- 
standes in einem großen Anstellwinkelbereich gut 
herauskommt. 
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P. und E. Brauer, Über unvollstän 
Webersche Funktionen. Z. angew. Math. 
Bd. 21 (1941), 8. 177. 


Herr Prof. Dr. Hugo Kasper, Brünn, wies uns 
freundlicherweise auf eine Unstimmigkeit hin, 
welche ihre Ursache in einer durch Schreibfehler 
entstellten Gl. (1) hat. Diese Ausgangsgleichung 
muß richtig lauten: 


dige Anger- 
Mech. 


Wilh. Müller, Beitrag zur Theorie der lang- 
samen Strömung zweier exzentrischer Kreiszylinder 
in der zähen en Z. angew. Math. Mech. 
Bd. 22 (1942), S. 177—189. 


Der Titel der Arbeit soll lauten: ER zur 
Theorie der langsamen Drehung zweier exzen- 
trischer Kreiszylinder in der zähen Flüssigkeit.“ 


Berlin NW 7. 


Für den Textteil verantwortlich: Prof. Dr. Fr. A. Willers, Dresden-A. 20, Dorotheenstr. 12. — Printed in Germany. — 
Copyright 1939 by VDI-Verlag G. m. b. H., 


— Druck von A. W. Ziekfeldt, Osterwieck am Harz. 


